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CHAPTER I 


HISTORICAL INTRODUCTION TO HADAMARD 
TRANSFORM SPECTROMETRY (HTS) 

The Idea of modulating or encoding the optical output of 
a spectrometer goes hack to the original work of Golay (1949) 
and Fellgett (1951). Its purpose is to allow many different 
wavelengths of radiation to fall on a detector simultaneously, 
and thereby to increase the signal -to-noise ratio (SNR) of the 
resulting spectrum. This improvement comes about because each 
element of the spectrum is effectively viewed a larger fraction 
of the total available observing time. One idea is to encode 
or modulate each spectral wavelength exiting the spectrometer 
output with an audio frequency that contains the optical wave- 
length information. The use of a conventional wave analyzer 
then allows recovery of the original optical spectrum. There 
are many variations of this technique. 

In 1968, Ibbett et^ al and Decker et al Independently 
suggested the use of sequentially stepped multiplex spectro- 
meters. In both systems radiation enters a dispersive instru- 
ment through a single slit and is analyzed at a number of exit 
slits. Decker et al pointed out that two constraints should 
be Imposed on the encoding scheme: 

(1) To obtain the optimum signal to noise ratio, each spectral 
element should be viewed during exactly half the step positions. 

(2) To impose the smallest dynamic range requirements on the 
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detector amplifier system, each step position should pass 
light from exactly half the spectral elements. They also 
worked out a scheme that satisfied the two constraints for 
masks having elements m=4n+2, where n is an arbitrary integer 
or zero. Ibbett et al introduced the Hadamard pattern for the 
mask. As discussed below, this is a pattern based on a set of 
binary orthogonal matrices first studied by the French mathe- 
matician Jacques Hadamard. Ibbett et_ al^ also described the. 
application of their scheme to a real time computer aided 
measurement . 


In 1969a Sloane et. al worked out a number of binary cyc- 
lic coding schemes for multiplex spectrometry and evaluated 
the performance of each scheme in terms of a linear, least 
mean square, unbiased estimate. These schemes include a Hada- 
mard matrix H, and various modified Hadamard matrices, which 
these authors refer to as G matrix and S. matrix (Fig. 1-1). 

A Hadamard matrix H of order N is an N x N matrix Hjj 
of +l‘s and -l r s which satisfies: 


where 1^ is an N x N unit matrix 

A modified Hadamard matrix G of order M is a partitioned 
matrix from the H matrix: 



(M = N-l) 


1 


H 


+ + + + + + + + 

+ +- **-+* + 
+ -- -+ - + + 
+ -- +- + +- 
+ “+ ta " + + w ” 
+ + “ + + — - 
+ - + +- -- + 


G 


+ + - - - + - 
+ *- — “+** + 
«w m *-|— 

-+- + +-- 
-■++-“- + 


0 0 1110 1 
0 1110 10 
1110 10 0 
110 10 0 1 
10 10 0 11 
0 10 0 111 
10 0 1110 


Figure 1-1. An 8 x 8 Hadamard matrix and two 7^7 cyclic 
matrices that can be derived from it. 



where the first row and first column of H are all +l’s. A 
feature of the G matrix is that it can be written in cyclic 
form — a factor which we will show to be of considerable 
practical importance. 

A modified Hadamard matrix S is a matrix obtained from 
G by replacing +l’s by 0’s and -l*s by +l f s. 

The properties of H, G, and S will be discussed in sec- 
tion 2-6. 

When we talk of encoding by means of a Hadamard matrix 
we have the following in mind. A mask is used to modulate - 
open or close - a series of entrance and exit slits in a spec- 
trometer. If a certain slit location is open, we can desig- 
nate it by a +1; if it is closed we can designate it by a 0; 
if it can be used to subtract from the signal incident on the 
detector, we designate it by -1. The sequence' of +l*s and 
-l’s characterizing a mask in a given modulating position 
corresponds to a row of a matrix. The whole set of mask posi- 
tions corresponds to the set of rows of the matrix. If the 
sequence of mask patterns corresponds to the rows of a Hada- 
mard matrix we say we are encoding with a Hadamard pattern. 

Sloane et al also introduced the idea of using a cyclic 
matrix for coding masks. This greatly decreases the experi- 
mental cost and facilitates operation, since any N slits of 
a single mask 2N-1 slits long can be used to provide one of 
the required mask patterns. 

The first single entrance Hadamard spectrometer (HTS) 
was .built by Decker and Harwit (1969). The spectrometer had 


a single entrance slit and 19 exit slits. The exit mask was 
stepped manually. The authors used this spectrometer to take 
the spectrum of the mercury vapor 1.7v band to demonstrate the 
Hadamard transformed spectrum’s fidelity and freedom from sys- 
tematic errors. 

With 19 exit slits, the HTS had a theoretical slgnal-to- 
noise advantage of 2.18 over the conventional spectrometer, 
which is rather hard to verify experimentally. Decker (1971) 
therefore proceded to build a 255-slit HTS. In this spectro- 
meter, the radiation, after being decoded by the exit mask, 
exits along the same path it comes in. This reverse pass de- 
disperses the beam and allows it to be brought to a focus at 
the entrance plane. Thus the dimensions of the focused image 
are roughly the same as the dimensions of the entrance aperture, 
and the detector size can be minimized. This is important since 
sufficiently large detectors sometimes do not exist, and if 
available tend to be noisy . Decker experimentally verified 
the theoretically predicted multiplex advantage of an HTS. 

DeGraauw and Yeltman (1970) were the first to use an HTS 
for astronomical work during the 1970 solar eclipse. Houck 
et al (1973) subsequently used an HTS to obtain near infrared 
spectra of Mars from airplane altitudes. 

Besides putting an encoding mask at the exit plane, one 
can also put another encoding mask at the entrance plane of 
a spectrometer. In that way the radiation is modulated at 
both the entrance and exit apertures. Harwit et al (1970) 
worked out this scheme of doubly multiplexed dispersive spec- 


trometry, The double multiplexing scheme allows one to in- 
crease the total amount of radiation that can be transmitted 
through a spectrometer, Furthermore, by a proper reduction 
of the data, one can also obtain a one dimensional picture 
of the source at the entrance plane. For a spectrometer of 
m entrance slits and n exit slits, one needs m x n data points 
to recover m spatial spectra, with each spatial spectrum con- 
taining n spectral elements. For a homogeneous source one. 
does not need the spatial information, so (n + m - 1) data 
points vrill be enough to recover the spectra. Harwit et al 
(197 i la)describe two schemes for recovering the spectrum with 
(n + m - 1) data points. 

In 1975 Tai et al (1975a) finished the construction of a 
doubly coded Hadamard transform spectrometer. 1 ,e spectro- 
meter -has 15 entrance slits and 255 exit slits, which can 
simultaneously obtain 15 spatial spectra, each having 255 
spectral elements. Tai et_ al (1975b) went on to give an ana- 
lysis of the errors in Hadamard spectrometry caused by imper- 
fect masks. 

Besides coding the radiation at both the entrance and the 
exit aperture, one can go one step further and use a two di- 
mensional mask at the entrance aperture (Harwit, 1973-)-> 

This yields a two dimensional picture at the entrance aperture, 
where each spatial point at the entrance has its own spectrum. 
To put it a different way, one obtains a two dimensional pic- 
ture of the source at the entrance aperture for each color of 
the spectral elements. 
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Harwit (1973) experimentally verified the ' operation of 
imaging spectrometry, and Swift et al (1976) constructed the 
first Hadaraard imaging spectrometer. 

There are other discussions of Hadamard transform spec- 
trometry in the literature, mostly of theoretical aspects. 
Nelson and Fredman (1970) give a more complete theoretical 
treatment of Hadamard matrix encoding. They also rediscovered 
a theorem due Initially to Hotelling (194*0 showing that the 
Hadamard matrix is the best design for a singly coding mask. 
Sloane and Harwit (1976) show the connection between Hadamard 
spectrometry and the mathematics of weighing designs in sta- 
tistics. 

There have been various comparisons of Hadamard transform 

spectrometry with other spectrometry. Larson et al (197*1) 

* 

makes a theoretical comparison of singly multiplexed Hadamard 
transform spectrometers and scanning spectrometers. Tney pre- 
sent a general mathematical framework for the comparison of 
relative performance and also verify their prediction by com- 
puter simulation of various characteristic spectra. Their re- 
sults show that where the noise level is constant and indepen- 
dent of the incident photon flux, the determined multiplex 
advantage is /N/2 , as predicted by Fellgett (1953.) . This is 
usually the case in a low energy region, such as the infrared. 
For a noise level that is signal-dependent, such as in the T3Y 
energy region, the detector is characterized by an output with 
statistics approaching a Poisson distribution and variance 
therefore proportional to the input signal. In that case the 
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HTS technique will be advantageous only for spectra that are 
characterized by a few- well-defined and intense peaks on a 
very low intensity background. For spectra with high back- 
ground, for dense spectra, or for spectra having very weak 
spectral features, the HTS will have no advantage over the 
conventional single slit (SS) technique. 

Hirschfeld and Wyntjes (1973) compare Fourier transform 
and Hadamard transform spectrometry,' They also describe 
various limitations of Hadamard -transform spectrometry. This, 
paper was followed by an exchange of notes between Decker 
(1973) and Hirschfeld and Wyntjes (1973) in the journal Applied 
Optics in which some of these limitations are disputed. 

These papers concern themselves with a number of practical 
matters on which opinions can vary. Here we mention these con- 
troversial papers mainly for completeness. Their contents will 
be discussed further below. 

Wyatt and Esplin (197*0 analyzed the effect of band width 
on noise equivalent power (NEP) for multiplex spectrometry 
with cryogenically cooled, cooled-background extrinsic long 
wavelength infrared detectors. They find that the NEP is 
directly proportional to band width, so multiplex schemes that 
require increased band width are not of real advantage. They 
further conclude that doubly encoded systems that are based on 
m + n - 1 measurements would have a real throughput advantage 
' Various other aspects of Hadamard matrices and Hadamard 
transform spectrometry which have not been mentioned above are 
covered in .articles by: Baumert, Pratt et al (1969)* Hirschy 


et al (1971), Allen et al (1972,1973), Kowalski et al (1973), 
Planky et al (197*0, Oliver et al (197*0. 

In this thesis Chapter II will describe the mathematical 
properties of Hadamard matrices and their application to' spec- 
troscopy. Chapter III describes the Hadamard transform spec- 
trometer, and gives results on laboratory performance. 

Chapter IY gives a comparison of Hadamard transform and Fourier 
transform encoding In spectrometry. The output of an HTS is 
fed Into a mini computer. The computer performs a real time 
inverse Hadamard transform to recover the spectrum. Chapter V 
describes the algorithm and programming of Inverse Hadamard 
transform. Chapter VI discusses observational results and 
their interpretation. 


CHAPTER II 
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HADAMARD MATRICES 


(A) Weighing Designs 

In order to understand the' mathematical advantage of 

Hadamard transform encoding, let us look at the following 

examples (Sloane et al, 1976 ). 

Suppose four objects are to be weighed, using a spring 

balance which makes an error e each time it is used. Assume 

2 

that e is a random variable with mean zero and variance a . 

First suppose the objects are weighed separately. If the 
unknown weights are ^i, hs ^33 #4> the measurements are tjj, 
n 2 , ti 3 s T14, and the errors made by the balance are ej , e 2s e3, 
et* , then the four weighings give four equations: 

m = $1 + ei r\2 = ^2 + e 2 

ti 3 = i|>3 + e3 ri4 “ ^4 + e 4 

The best estimate of the unknown weights are the measurements 
themselves: 

a 

♦ 1 ** ni 0 ti + ei 

A 

Tfi 2 = ri2 = ^2 + e 2 
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These are unbiased estimates: 

Ei|> i = ip i 

a 

Ei |»2 = 1 P 2 (E denotes expected value) 

with variance or mean square error 

* 2 2 2 
E(ipj — ipj ) ~ Eo i — a 

On the other hand, suppose the balance is a chemical 
balance with two pans, and the four weighings are made as 
follows : 


ni 


'h 


*2 

+ *3 


ip4 

+ 

©l 

t2 

= 


“ 

^2 

- ^3 

- 

t|> 4 

+ 

e2 

T13 

' = 

V i 


^2 

- ^3 

- 

K 

+ 

e 3 


= 


- 

^2 

- if»3 

+ 

^4 

+ 

e 4 (2-1) 

means 

that 

in 

the 

first weighing 

all 

four 

objects are 


placed in the. left hand pan, and in the other weighings two 
objects are in the left pan and two in the right. (Note that 
the e are independent of the weights on the balance. This 
point is crucial). Xt is easy to solve for ih s $3* Kj 
as long as the . coefficient matrix for tp is not singular. 

Thus the best estimate for ipi is 





I 
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♦ l 


^(ni + ti2 + n 3 

if»l + 


+ m*) 

+ e 3 + e 4 ) 


. o 

The variance of Ce, here C is a constant, is C times the 
variance of e, and the variance of a sum of independent ran- 
dom variables is the sum of the individuals variances. 

A A A A 

Therefore the variance of ifii(and also of ijia, 1(13, 1^4) is 
,22 
4g _ g 

16 4 


Weighing the objects together has reduced the mean square 
error by a factor of 4. In effect the signal to noise ratio 
(SNR) 9 which is given by the root mean square (rms) error is 
reduced by a factor of 2. 

Finally, suppose the balance is a spring balance with 
only one pan, .so only coefficients 0 and 1 can be used. A 
good method of weighing the four objects is: 


ni - 


^2 

+ $ 3 

+ *4 

+ 

ei 

T12 = 

¥1 

+ ^2 



+ 

e 2 

ns = 

* 1 


+ $ 3 


+ 

e 3 

T14 = 

4*1 



+ ^4 

+ 

*4 


In this case the variances of i|>i, $39 are — g— , -q—j 

22 •> y 

—5—, -respectively, a smaller improvement than in the pre- 

y y 

vious case. 

The theory of weighing designs is of immediate interest to 
multiplex optics, .since the simultaneous measurement of the 
Intensities of different bundles of rays is completely ana- 
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logous to the simultaneous weighing of different groups of 
weights. In measuring the intensity of radiation passed 
through slits in a mask, we are effectively ’•weighing’ that 
radiation. 


(B) General Mathematical Formulation 

One can put the problem into a more general form. 

Let be the 1* unknown, be the 3 measurement, be 


the error associated with the 3 


th 


measurement. Let w*^ be 


ii. it. 

the weighing coefficient of the f 1 measurement with the i n 


unknown. Then 


n 


3 


*31*1 



i=l n 

3=1* * • • -n 


(2-3) 


In matrix notation: 

n = W f + e ■ (2-4) 

With the notation < > for ensemble averages, the error e^^ has 
the following properties: 

(1) <e ± > = 0 

(2) e^ is independent of ^ 

(3) <e i e j > * 0 if errors are assumed to be 

uneorrelated. 

= a 2 if 1=3 

The problem now is the following: (i) For a particular coding 
matrix W, what should be the decoding matrix A? i. e. What Is 
A such that j£=A where is an unbiased estimate of 


V 
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(ii) What is the best choice of W (or A) that will minimize 
the error, of measurement i.e. What is W such that 

e = * f-i ' V* ” 


is a minimum. 

In the absence of noise, i.e. n=0, it is clear from (2—4) 

that 

& * - w" 1 ^ 

and therefore 


A = W" x 

and 

i = £ il “ t 

In the presence of noise. 


* = A/tl 

= A W ^ + A _n 


and with the assumed properties of the noise <£>= one ob- 
tains 

<tjj> = A , W <_£> + A<e> 

« A|i 


Assuming no prior knowledge of the unknowns, one may use the 

a 

unbiased condition <£>= This again implies 
A - W“ l 

So with the assumed properties of noise and unbiased con- 
dition, the decoding matrix is just the inverse of the coding 
matrix. 

One still has to find a coding matrix which will minimize 


the uncertainty, e . 


The second question can be solved in the following ways 

f; Vi 

Let be the i measurement in the absence of noise, then 
for each measurement 


*1 


*11+1 + *12+2 

+ ••• + W. ib + e* 

in n i 


■= 

"i + e i 


*3 

ss 

A ji n l_ + A 32 i, 2 

+ "\ + V>n 



A jl <+ 1 + e l 5 + 

A 32 ^2 + e 2> + ••• ' 



+ A jn (n n + e n> 



= 

(A^i^i + ••• + 

A JnV + < A j2 e 2 + *'* + A ^n e n^ 


=2 

+ noise 



The mean square of the noise term corresponding to the j 
unknown is therefore 

*3 = <*3*1 + •" + A !n> ° 2 

" i 2 a 2 (2-7) 

where 

A.j = CA Jl + * A jn^ C 2 " 8 ) 

and Aj represents the improvement in the SNR for the weighing 
design, compared to the SNR for individual weighings. 

Hence, the problem of maximizing the signal to noise 
ratio becomes the problem of minimizing e^ ; or A^ (Nelson and 
Fredman (1970)). 

Sloane e£ al (1969) 'independently developed an expression 
for b/c 2 where 
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e/ff 2 = Trace W“ 1 {W“ 1 ) T 

- Trace A- A T ' . (2-9) 

Equation (2-8) and (2-9) amounts to the same thing be- 
cause it can be seen very easily that 

* 

m n ? 

Trace A A - 

The question of minimizing e had been answered by Hotelling 

(Hotelling, 194*0 and rediscovered by Nelson and Fredman. 

Hotelling has shown that for any choice of mask W with |W. }<1, 

2 J 
the are bounded by e i > , and that it is possible to have 

e i =^r for i«l,... N if and only if a Hadamard Matrix of the 

order N exists (by taking W=H^) . This leads to the discussion 

of the Hadamard Matrix. 

(C) Hadamard Matrix 

A Hadamard matrix of order N is an NxN matrix H^ of +l’s 
and -I’s which satisfies: 

= HI n (2-10). 

where I N is an NxN unit matrix. 

A Hadamard matrix has following properties: (Golomb(1964) ) 

(1) Its row vectors (or equivalently, its column vectors) 
are mutually orthogonal. 

(2) The Hadamard’ properties will not be disturbed by:. 

a. Interchanging rows, 

b. Interchanging columns, 

c. Changing the sign of every element in a row, or 


IT 


d. Changing the sign of every element in column. 
These properties enable the first row and column of 
every Hadamard matrix to be normalized to contain only +1*8. 
If G represents the remaining M x M matrix (M= N-l), then H 
can be partitioned into 


1111 

1 


(M= N-l) 


L 1 ' J 

It is conjectured that Hadamard matrices exist for all 
multiples of four. Further, if one of the following conditions 
is also satisfied, . 


(1) N « p + 1 

(2) N = P(P + 2) + 1 

(3) n = 2 m 


F prime 

P and P + 2 prime 
m an integer 


,th 


then G can be made cyclic. That is, the (j + 1) row can be 
generated by shifting the y row one position to the left. 
For example, when N-8, we have matrices of the form shown in 
Fig. 1-1. Note that H and G are symmetrix. 

Another choice for W is the matrix S obtained from G by 
replacing +l T s by 0’s and -1’s by -rl’s. 

The properties of 'the H, G and 2 are discussed by Sloane 
et al. If rows I and j are any two rows of H, G or S, It can 
be shown that their dot product is: 


In H n : row ' i • ' row j 


0 i/j 
N i=j 


In : row i • row j ® 


-1 i^j 
M 1=3 


In S M : row i * row j 


WH 1£3 

M/2 1=3 


The inverse of each matrix is: 


% 1 N — N 1 % X 


ji { r; T ^ . q“^- 

m+i'% * %i 




where J is a M x M matrix consisting entirely of -1 ' s and 
N» M+l. 

Table 2-1 gives the value of A for different matrices. 
The matrix I_ represents the weighing scheme weighing each 
object separately. This corresponds to a conventional single 
slit spectrometer or to a wedge filter monochromator. • 


Table 2-1 


MATRIX A 

ELEMENTS OF A 

A" 1 

(FOR LARGE 

I 

1,-0 

1 

1 

H ' 

1,-1 

/N 

/N 

G 

1,-1 


if 

S 

1, 0 

^(2-f)” 1 

2 


cosine sq.uared 
functions 

N 2 
¥ N+l 

ffi 

Jts 


* F is the Fourier Transform case which will be disscussed in 
Chapter iv. 
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If the number of measurements N is a multiple of 4, and 
the matrix coefficients are +1, the best weighing scheme is 
the Hadamard matrix H. e will be reduced by a factor =— 

compared to weighing the unknown separately. This is the 
maximum advantage a weighing scheme can obtain with weighing 
coefficient l^jl-l* If N is. not- a. multiple of. 4, or .if the 
weighing coefficients are 0’s and l*s, it is not possible to 
simultaneously minimize and some other ■ criterion 

must be used (Sloane et al, 1976). Also the errors are uni- 
formly larger than for the H-matrix, as shown for the G and 
S matrices in Table 2-1 above. 

It is' interesting to see that a spectrometer using the 
Fourier Transform, such as a Michelson interferometer, has a 
multiplex advantage a factor of /S’ lower than the H-matrix 
and a factor of /2 lower than the S~matrix encoding instru- 
ment. 

Following are some computer simulations for S-matrix 
transformations with various inputs (See Fig. 2-l(a) to (c)). 


. INPUT 

1. Constant 

2. Hadamard code: representing 

single line emission 

3. Single line: 1 at the 1 st 

element and 0 for the ’rest . 
This 2’epresents an unknown 
impulse coming in during the 


OUTPUT 
Constant 
Single line 

Hadamard code. 

Note, unlike the mono- 
chromater, the error 
propagates to other 


INPUT 


OUTPUT. 






20 


observation. 

4 . Sine wave . 

5. Square wave. 

The Input is not a perfect wave 
because we have an odd number of 
elements. The input values have 
amplitudes 0 or 1 for each element. 

6. Straight line at a slope 1/255. 
This may represent a shift in 
baseline. 


elements. 

sine wave with dif- 
ferent phase. 

Hot a perfect square 
wave. 


Refer to figure (2-1C). 


(D) Optical Realization of Hadamard Encoding 

We have made use of (modified Hadamard) S-raatrices in two 
optical instruments: One is a Hadamard transform spectrometer 

having an encoding mask at the exit aperture. The other is a 
doubly encoded HTS which has encoding masks at both the en- 
trance and exit apertures. 

S codes can be used for both the entrance and exit masks 
for the HTS, with +1 standing for an open slot through which 
radiation is transmitted, and with 0 standing for a closed 
slot where radiation is blocked. 

The cyclic property- of the S-raatrix is very desirable, 
for then only a single mask 2N-1 slots wide need be constructed 
Successive encoding positions are generated by stepping the 
mask one slot width along its length. This avoids the con- 
sturction of N masks with N 2 slots. 
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Although H and G matrices have better coding efficiency 
than the S matrix, they introduce technical difficulties when 
used for spectrometer coding schemes. To utilize H and G ma- 
trices, one must measure the reflected as well as the tran- 
smitted radiation. In this case +l's represent reflecting 
slots and -l's represent transmitting slots. Therefore a mi- 
nimum of two detectors must then be used, one in a subtracting 
mode, the other in the normal mode. The use of two detectors, 
however, increases the noise. Furthermore, the H matrix, with 
all elements +1 in the first row, makes the dynamic range of 
the detector system change by about a factor of two. Also its 
lack of the cyclic property does not allow one to generate the 
rest of the masks by the simple stepping technique mentioned 
above. 

1. Hadamard Transform Spectrometer (HTS) 

For the singly encoded HTS, we use the following optical 
arrangement (figure 2-2). Radiation passing through the single 
entrance slot is rendered parallel and directed to the dis- 
persing element. The dispersed radiation is then de-collimated 
and focused upon the multi-slot mask at - the- exit aperture. 

The spectral elements transmitted by the mask pass through 
suitable post-optics and are collected onto a detector. One 

then makes N (in our case N=255) measurements by sequentially 

* • 

stepping the mask H times. The Inversion procedure ip=S £ 
recovers the spectrum. Figure (2-3). gives a 255 cyclic S- 
matrix code for the exit mask. 

Theoretically, the mean square error in the spectrum 


Single entrance tilt 


Exit mask 
In position j 



Schematic representation of HTS using the 
S code. 








ii ii j i ii imi in i iiiiiiiiiiiiiiMifiiiiiiiiifiaiiiiiiiiiHiii mu 


00000 00101 
01001 11110 
01011 11110 
01001 01101 
00101 11010 
01101 01000 
01001 00011 


10001 11101 
10101 oino 
11110 01101 
00011 00111 
11110 11011 
00100 11101 
10001 


ooooi mil 
00001 10001 
11011 10010 
00111 10001 
11100 00110 
10010 01001. 


11001 00001 
01011 00110 
10100 10100 
101 10 00010 
10011 01011 
loooo ooiii 


Figure 2-3 . The 509 exit slit mask and the 255-element 
S code. 


2 ? 


given by a monochromator is or 2 . For an HTS using the S code, 

>69). 

_ VN 


1 ? 

this error is j|(2-^) 2 <j 2 (Sloane et_ al, 1969). Hence the rms 
gain in S/N for the HTS is G= — — 5 - . For N=255, 


2-1 
* N 


G^8.0 (figure 2-4 gives Decker’s results. The experimental 
gain was measured as 8.0±0.3). Note that the scale in figure 
(2-4) are in arbitrary units, and the zero point appears to be 
shifted between parts (a) and (b). 

2. Doubly Encoded Hadamard Transform Spectrometer (DBTS) 

Figure (2-5) is a schematic representation of an optical 
system which has a number of entrance as well as exit slits. 
Instead of passing radiation through only one entrance slit, 
a mask M slits wide is placed at the entrance aperture. 

Radiation passed into the spectrometer through different com- 
binations of open and closed slits. The dispersed radiation 
at the exit plane is analyzed in the same fashion as in the 
HTS. Encoding is accomplished by sequentially stepping one 
of the masks through its N different positions for each posi- 
tion of the other mask. 

In a DHTS the entrance aperture is modulated by a P x P 
S matrix-. 

Let e=e . be the P x P matrix whose rows represent P 
ir 

different entrance masks. e ir =l for open slots and 0 for closed 
slots (liiiP, lir£P). Similarly let represent the exit 

mask. When the entrance mask is in position i and the exit 
mask is in position 3, the detector measurement u.y is 


Intensity in Arbitrary Units 
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Mercury Emission Spectrum - Monochromator 



Mercury Emission Spectrum - HTS Spectrometer 



Wavelength, (>t) 


Figure 2-4 . Comparison spectrum of the mercury emi- 
ssion lines in the 1.4-1.8ym region: (a, bottom) 
as obtained in the Hadamard-transform mode, (b, top) 
as obtained under identical conditions using the 
same optical system as a conventional monochroma- 
tor. This figure is taken from Decker, 1971b. 


Entrance mask 
in position I 


I 


Exit mask 
in position ] 


Figure 2-5. Schematic representation of DHTS 
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P N . 

r=l 1=1 ir rs * 6j 


+ V 


ij 


( 2 - 11 ) 


where is the spectral element produced by radiation pas- 
rs 

Vi j- 

sing through the r vn entrance slot and the S a exit slot, v 


ith 


ij 


is the noise in the (i,j) measurement; it has the following 
properties : 


<v ij > = 0 


“ 0 s ikV 


If the instrument has no optical magnification, the spec- 
trum of radiation that passes solely through the r th entrance 
slot to first order is shifted by r spectral channels from the 
spectrum passing solely through the first entrance slot. 

Hence, only P+N-l distinct spectral elements exist: 


♦ 


-(P-l)»* 


• • • j ^ q $ ^1 * ° " V* 


N-l 


where 


r rs 


♦ 


r-s 




(t=r-s ) 


In matrix notation one may write 


jl 


T 

e i X + v. 


( 2 - 12 ) 


Employing the same analysis as one does for the HTS, i.e. 

A 

using the unbiased condition and the properties of v. 
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one obtains 

A 

where 





-1 


1 


*0 



*1 


^P-l ^F-2 ■ ’ * 




■N 


^-N-P+2 


a 

Each row i of jp represents a spectrum at the exit mask 
for radiation that enters the instrument through the i en- 
trance position. Hence the j diagonal gives a one-dimensional 
spatial picture across the entrance aperture for the spectral 
element j . 

A 

One may obtain an average spectrum ^ by averaging all the 
elements in each diagonal. 


* 


i 



P 

l 

r=l 




r s r-t 


t 1 0 


1 

5Rt"[ 


P 

s 

r=l 


*r 


s r-t 


t < 0. 


Harwit at al (1970) shoived that if both the entrance and 
exit masks are S matrices and we define 


a 


2 

t 
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Then 



1 / 


16 N 2 -l 2 
(N+l) N-]~t | ° 


l6o 2 

(N- | t | )N 2 for N large 


where 

t = -(N-l), (N-l) 

If the total mean square error for the unknown is 


e 



a 


2 

t 


( 2 - 13 ) 


where one sums only the central element-, then for the S-code 


N large 

-i N 

where o z = constant ^ • 

There are two points that should be made about the DHTS. 

(1) It has not been shown that Hadamard codes are the best 
codes for such an instrument. In fact some evidence suggests 
that Hadamard codes are not precisely optimum for this "two- 
ended” operation (Harwit et al a !974b) . (2) For PxN data points 

the spectrum yields only P+N-l spectral elements 3 plus a one- 
dimensional image. It is also possible to reconstruct the 
(P+N-l) elements with (P+N-l) data points only (Harwit et al, 
1974b). Fig. (2-6) gives a 15 element cyclic S-matrix code for 


°(i I ) 

* N 2 


the entrance mask. 
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Table 2-2 compares three different grating spectrometers. 
The first column represents the conventional single entrance 
and exit slot instrument. N measurements are made in time T, 
with a mean square error a 2 in each. The second column is 
for a singly multiplexed instrument with an exit mask S, and 
is taken from Sloane et al (1969). The last column is for the 
doubly multiplexed system, using Equation (2-13) for a 2 and 

U * 

has been multiplied by a factor of N to allow for having to • 
make N 2 measurements in time T. 

Table 2-2 

NO MASK • SHT DHT 

,-1 ^ N 

a 1 5* 

* /227T 

(E) Errors in Hadamard Spectroscopy 

During the manufacture of the masks , whether by deposi- 
tion of metal or by removal of metal through an etching process 
it is possible to obtain a systematic error that leaves each 
of the opaque portions of the mask either too wide or too • • 
narrow by a -fixed amount. This will cause a systematic vari- 
ation in signal passing through the slit. For example if the 
open slit is too narrow by a fixed amount e, the light passing 
through an open slit position will be 

I when the open slit is bounded by two open 

slits. 


I 0 (l-e) when the open slit Is bounded by one open 
slit and one closed slit- 

I Q (l-2e) when the open slit Is bounded by two closed 

slits. 

A similar analysis holds for open slits that are too wide, 
except that the minus sign in these expressions is replaced 
by a plus sign. 

The spectrum of a single (spectral line resulting from’ 
such imperfect masks) Is remarkably simple (Tai et aJL, 1975 b). 
Independent of the particular S-matrlx mask to be used, there 
are always precisely four false blips present in the final 
spectrum. The amplitude of these blips Is always the same for 
a fixed narrowing or widening of the transmitting slits. Two 
of the blips always surround the main spectral line, and a pair 
of adjacent blips always are some distance removed from that 
line. For the 255 element S-matrix, these two ar^ located 24 
and 25 elements away to the left of the parent line. The am- 
plitude of the displaced blips is positive when the transparent 
slits are too wide and is negative when the slits are too nar- 
row. In contrast, the two 'blips surrounding the parent line 
always are positive. -Figure (2-7) shows the negative features 
accompanying the 1.7 pm mercury vapor doublet and the computer 
simulation of a pure spectral line input and its distorted spec- 
trum. For the general case the reader may refer to Tai et al 
(1975 b). 
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Figure 2-7. (a) Spectrum of the 1.7vm mercury vapor doublet 

showing negative peaks to the left at the emission peaks; 

(b) Shows the response we would obtain to a single spectral 
line with a perfect mask; (c) Shows the response for a single 
line with the radiation simulated as passing through a mask 
with slits too narrow because each opaque mask element pro- 
trudes into the adjacent transparent slot by a tenth of a 
slot width; (d) Shows the effect of simulating slits that are 
systematically too wide. Note that the main spectral line 
has been placed in different positions for the synthetic . uns 
(b),(c), and (d). 










CHAPTER III 
INSTRUMENTATION 

(A) Multislit Spectrometer 

A conventional spectrometer has four essential elements, 
an entrance slit, a dispersive device such as a grating or 
prism, a set of imaging optics, and an exit slit. 

Such a spectrometer has two important parameters. The 
first is the "resolution” R, which is a measure of how well 
the spectrometer can separate two neighboring lines. The 
second parameter is the "throughput" E. This is a measure, of 
the light gathering capability of the system. The two para- 
meters, R and E, are lumped together into what is called "lu- 
minosity" L, defined as (Vanasse, 197 * 0 - 

L * E • R - 


For a conventional grating spectrometer having a grating 

of area A given by V7H, where V7 and H are the width and height 
& 

respectively of the grating, the throughput is determined by 
the product of A with the solid angle r subtended by the slit 
at the collimating mirror (or lens). The solid angle is given 
by* 

n B Si 
F 2 

where w and 1 are the width and height respectively of the slit 
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and P is the focal length of the collimating mirror 

E w • H • (3-1) 

F 2 


The resolution of this instrument is 


R = 



« nN 


n 


¥ 

• — 

d 


(3-2) 


where X is the wavelength, AX the closest wavelength that can 
he separated, n the order, N the total number of lines on the 
grating, W the width of the grating and d the spacing between 
rulings . 

Since d sin a - nX (3-3) 
substitute (3-2) into (3-3) 


| sin 


nx 


(3-*> 


if the slit width is limited by diffraction, which is the mi- 
nimum slit width, then 


w ^ 


FX 

¥ cos a 


substituting (3-5) into (3-*0 one gets 


(3-5) 
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* Tjl 

R 'v- — v an a (3-6) 

w 

Comparing equation (3-1) and (3-6) one sees immediately 
that, for a fixed grating area W*H, and fixed optical system, 

E is proportional to the slit width w and R is inversely pro- 
portional to the slit width w. This means that an increase in 
luminosity of the system by increasing the slit width is made 
at the sacrifice of resolution, and vice versa. 

Prom (3-1) and (3-6) one obtains 

L ^ E • R 

* W • H * — 

P 


Another feature of a conventional spectrometer is that it 
transmits only one narrow spectral* range of light to the de- 
tector, and all other spectral elements are wasted. As a re- 
sult, the instrument is inefficient. 

Within the past two decades there has been much research 
done in an effort to design new spectrometric systems with a 
view to maximizing the luminosity L, and to observe a number 
of spectral elements simultaneously. This can provide a mul- 
tiplex advantage, or a wide aperture advantage. Two quite dis 
tinct modulation techniques have been employed in the past. 

The -first depends on the wave nature of radiation, and makes 
use of interferometry. The Fabry-Perot interferometer, 
Michelson interferometer, and Mach-Zehnder interferometer 
(Jacquinot, 1954, I960; Vanasse and Sakai, 1967) are instru- 


ments of this type. The other technique employs dispersing 
spectrometers in which entrance and exit slits are replaced by 
opaque . or transmitting masks. Golay's .multislit spectrometer, 
Girard's Grill spectrometer and Hadamard spectrometers (Harwit 
et al, 1974a) are representative of these instruments. 

A spectrometer, whether interferometric os mask-multiplexed, 
yields a multiplex advantage mainly for detector noise or ampli- 
fier noise limited applications. In these cases it can be shown 
that for N spectral elements, one can achieve of the order 
improvement in the overall spectral signal-to-noise ratio, 

S/N, over a conventional spectrometer (Chapter II). 

For photon noise limited applications, the multiplexing 
advantage is cancelled by the N-fold increase in the photon 
noise attributed to the N-fold increase in the energy falling 
onto the detector. Nevertheless, for photon noise limitations, 
the throughput advantage can still be realized. The large 
throughput will become a disadvantage when the noise is back- 
ground noise whieh increases faster than the noise due to the 
source (Harwit ejfc al, 19711a) • 

In this chapter we will describe the experimental study 
of a Hadamard transform spectrometer (HTS) and calibration in 
the laboratory. Figure (3-1) is the flow chart of the whole pro* 
cess, starting with radiation from the telescope and ending 
with the output of the computer. Each component will be des- 
cribed. 







Radiation from telescope 


Chopper 

Square wave modulates 
radiation 


Chopper radiation 


HTS 

encodes modulated 
radiation by 
Hadomard code 


Encoded chopper radiation 


Alignment unit 


Cryogenic System 
Converts chopped 
radiation into AC 
voltage signal 


A.C. voltage .at chopped frequency 


Princeton Applied Research 
Synchronus demodulator 
changes AC. voltage to 
DC. voltage 


D-C. voltage proportional to incoming 
A.C. signal 


Voltage tc frequency converter 
generates pulses at a frequency 
proportionol to voltage 


Frequency proportional to DC input 


Monsanto Counter 
_ Counts .the number .of pulses 
Generated in a tixed time interval 


Outputs the number of pulses 


Digital Computer 

Performs inverse Hadomard transform 


Final spectrum 


The flow chart of the data taking process, 
starting with radiation from the telescope 
and ending with the output of the compu- 
ter. 











(B) The Opties 


1. Spectrometer 

Pig. (3-2) shows the basic spatial design of the spectro- 
meter. It works in' the Ebert-Fastie mode. Radiation passing 
through the entrance aperture S falls upon the spherical mirror 
M which, In turn, renders it parallel and directs it to the 
grating G. The dispersed radiation is collimated by the other 
half of the spheroid and focused upon the exit aperture S’.. 

A 255-Blot encoding mask is located at this position. The exit 
focal plane Is positioned in such a way that it bisects a 90° 
corner reflection. The corner reflector returns the radiation 
through the spectrometer again and displaces the beam from 
the center of the principal plane to one side. This reverse 
process dedisperses the beam and allows It to be brought to a 
focus at the entrance plane. The dimensions of the focused 
image are roughly the same as the dimensions of the entrance 
aperture (Decker 1971). These procedures allow one to use a 
smaller detector. 

A diagonal mirror at the entrance directs the dedispersed 
radiation to the liquid helium cooled post optics. 

The entrance mask can be a single entrance slit with any 
width between zero to 1.5 mm for the 1x255 program, or it can 
be a fifteen S-matrix code with each slit having width 0.1 mm 
for the 15x255 element program. In normal use the height of 
of entrance slit is 3.5 mm. It can be Increased up to 10 mm. 

is a spherical mirror with a 49 • 5 cm focal length. 

On its back it is held in place with three teflon-tippled 



Post Optics 


Figure 3-2 . Optical path through the spectrometer. 

The dedisperser and exit mask are shown 
rotated by 90°. 
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screws. Three teflon-tippled springs bear on the front edges 
of the mirror. This design allows one to make slight adjust- 
ments in the position of M when aligning the instrument. The 
central part of the mirror is blocked off to reduce stray ra- 
diation. 

G is a 75 m x 75 J® grating with 20 lines/mm and blaze 
angle 5 0 11*, The corresponding blaze wavelength at first order 
is 9*03v. Figure (3-3) is the calculated grating efficiency as 
a function of slit position s at two wavelengths . At 8y the 
energy imaged within one slit width is B0% of the total. For 
1% the energy within one slit is 52%. The grating is mounted 
on a yoke which allows it to be adjusted in three mutually per- 
pendicular directions. It is located at 0.82 focal length from 
the primary mirror M. 

All mirrors inside the spectrometer are silver coated with 
a protective coating of S1C> 2 . The reflectivity of silver coat- 
ing at 10y is better than 97%. 

For a multiplexing spectrometer which has N entrance and 
N exit slits j one wishes to image entrance slits S n 3 S 2a 

t 1 1 

S onto exit slits S- , S 0 , S such that S-, is imaged onto 

n l J d* J n l 

t T 

...... a S n onto S n at a particular wavelength X. Let 6 be 

the angle subtanded by S. Where 6 and 6' a.^e measured from 
the center of M. The grating equation for imaging S onto S' 
is 


Sin a + Sin 8 


mX 

a 


(3-7) 


in 

-sr 



t 



Linear Distance from Center of Central Slit (mm) 


Figure 3-3 . Grating diffraction efficiency as a function of slit position 
~ at 8 and l4ym. 
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Differentiating (3-7) with respect to a gives 

dB _ _ Cos a (3-8) 

da Cos B 

The minus sign indicates that a and B change in opposite direc- 
tions. da is the width of the entrance slit and dB is the width 
of exit slit, so unless a=B. 

In our spectrometer, the exit slit width is 0.1024 mm, 2.4# 
larger than the entrance slit width. This is the effect of 
anamorphic dispersion- a magnification produced by the grating 
(3-8), when the lower limit on the slit width is set by dif- 
fraction. The total number of spectral elements that can be 
observed simultaneously is limited by the optical aberrations 
of any particular optical system, which set a limit on the 
total useful width over which the spectrum can be displayed. 

2. Post Optics 

The post optics consist of a liquid helium cooled Arsenic- 
doped silicon (As: Si) detector, with appropriate optics for 
focusing the radiation onto the detector (Pig. 3-4). Radiation 
enters the evacuated dewar through a barium fluoride window,' 
passing through a filter with pre-selected band-width. The 
filtered radiation then passes through a cooled barium-fluoride 
filter and is focused onto the detector inside the housing by 
a gold coated mirror. A light baffle Is partitioned in front 
of the detector housing. 

Belov; is a brief discussion of each of the cryogenic 


components 



Figure 3-** . Liquid 
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(a) Barium Fluoride Window 

The barium fluoride window is 2 mm thick and 0.75” in dia- 
meter, It is u&ed to cut off wavelengths longer than l4y. It 
has a transmission efficiency of around 90£ out to lly before 
it starts to cut off. At l4.29y (700 cm ), its' transmission 
efficiency is 50$. 

(b) Filters 

A liquid helium cooled filter wheel with 8 filter posi- 
tions is housed inside the dewar. Table 3-1 gives a brief 
description of each filter. 




Table 3-1 


Filter position 

Range 

Band-width 

Peak transmission 
efficiency 

1 

no filter 


2 

closed 
(blocked by 
aluminum foil) 


3 

8y - 

9.9y 1.9v 

.7155 

4 

8.7v " 

11. ly 2.4y 

83S 

5 

lly 

12. 4y 1 . 4y 

82/5 

6 

11. ly - 

13. 8y 2.7u 

8655 

7 

8 . 4y - 

15y 6.0y 

855? 

8 

glass 

shorter 

than 2y 


Filter positions 1 and 2 are for testing purposes. 

The filter wheel is held in place by a spring-loaded screw 
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The wheel was originally connected to the outside world through 
a stainless steel rod and could be changed to different posi- 
tions by turning the rod. It was found that the- stainless 
steel rod conducts too much heat from the outside into the 
helium containing can. When the stainless steel is replaced 
by a G-10 Glass Epoxy Lamitex rod, the holding ’time for the 
liquid helium of the dewar Increases from 9 hours to 15 hours. 

(c) Field Mirror 

A gold coated mirror with focal length 7 nun and f/o«4l is 
used as a field mirror to focus the radiation onto the detector. 
The reflection efficiency for gold mirrors at lOy is over 99 % • 
The mirror has 3 degrees of freedom of adjustment, one trans- 
lational and two rotational adjustments. 

(d) Liquid Helium Cooled Barium Fluoride Filter 

This barium fluoride window is also used to cut off the 
radiation longer than l4u . Although the filter barium fluo- 
ride window cuts off radiation longer than l4y from the out- 
side world. It will emit radiation of Its own because it Is 
at room temperature. Since all the interference filters have 
a long wavelength leak between 20y to 26 y, and since the de- 
tector will cut off radiation longer than 24u only, there is 
still radiation from 20y to 24y that gets into the detector 
as' background radiation. The insertion of the liquid helium 
cooled barium fluoride filter eliminates this peak. It was 
found to cut down the background radiation by a factor of four. 

(e) Detector 

An arsenic doped silicon detector with dimension 1.2 mm x 


4 


.1 
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3.2 mm is positioned inside a housing, which has a baffle at 
its entrance. 

Arsenic doped silicon is a N-type extrinsic semiconductor. 
When the detector absorbs radiation, free carriers are provided 
for the conduction band, thus changing the reslstence of the 
detector. The following discussion follows the work of Putley. 
For further details, one can refer to Putley (1964) and Kittel 
( 1966 ). 

Let c be the conductivity of the detector 
e be the electric charge of the carrier 
t be the life time of free carriers 
N be the density of free carriers, and 
y be the mobility of free carriers 
Then, 

c D = NePy (3-9) 

and 

Acjp = AJexyn • 

= eryrj (3-10) 

where 

ACp is the change in conductivity of the detector 

AJ is the number of photons incident in unit time 

AP is the radiation power incident 

v is the frequency of the incoming photons, and 

n is the quantum efficiency, i.e. nmuber of electrons 

freed per photon. • 
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Since 



where R^ is the resistance of the detector 
£ is the length of the detector, and 
A is the area of the detector 

Then 


( 3 - 11 ) 


AR 


D 


_£ 

A 

£ 

A 


AP n 
hv N 2 eyx 


( 3 - 12 ) 


Let 


V. 


B 


Rr 


( 3 - 13 ) 


where I D is the detector current, and 

V B is the bias voltage across the detector 


Then 


AX 


D 


B 


D 


iR D 


Let 


V, 


V. 


B £ AP 


n 


( 3 - 14 ) 


R^j- A hv N 2 eyT 


Vl 


R, 

— — v. 


r t 


B 


( 3 - 15 ) 


52 



i 

! 

j 

[ 

| where V is the voltage across the load resistor 
I o 

) 

1 H t is the load resistance 

! L 

Then 

! iV o - “d*L 

i 

! 

( 

R L V B hv* T T,epT (3-16) 

From equation (3-16) one can calculate AP from AV Q . 

(f) Procedures for Alignment of the Optics 

i) Place all components in their respective positions 
and line them up visually. Be sure there is no mechanical 
binding in the mask and in the driving mechanism. 

ii) Using a laser s put the spot from the entrance slot 
on the middle of the grating. It is suggested that only one 
entrance slot be used. 

iii) Adjust the grating tilt until the line of dispersed 
dots exits at the proper position at exit. As the grating is 
rotated 9 this line of spots should remain level, not displaced 
normal to itself. 

iv) Put one half of the dedlspersing mirror combination 
in place. Use a T square to line it up roughly. At this point, 
use a mercury emission lamp with proper f-number to simulate 
the beam coming from the telescope. A number of colored image 
of the single entrance slot will be seen at the exit. 

v) Adjust the spherical mirror for coarse adjustment 
and the position of the dedispersing mirror as a fine adjust- 

f ment to bring the image to a focus on the exit plane. 

i'- 

ll 

vi) Adjust the grating position in its yoke until the 

Is ■ 

1 : ■ , 
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color image from the mercury lamp is parallel to the exit slit 
length. 

vii) Adjust the angle of the dedisperser so that the 
color image reflected by it is perpendicular to the exit mask. 

vlii) Put in the other half of the dedisperser and line 
it up at an angle so the radiation falls back upon the grating. 
Adjust with fine adjustment screws so the grating is fully 
illuminated. The image of the grating will appear on itself 
when viewed from the diagonal, 45° * mirror which diverts the 
radiation to the dewar. 

ix) Place the dewar on the spectrometer using V* spacers 
to represent the thickness of the dewar bottom cover. Adjust 
the 45° mirror so that the dedispersed image (with color) falls 
on the center of the filter on the filter wheel. 

x) Turn the filter wheel to position one (no filter), 
so that radiation can fall on the field mirror. Adjust the 
field mirror until the dedispersed radiation impinges upon the 
detector. Be sure that all the light falls onto the detector. 

Be sure that the mirrors accept all the radiation. 

xi) Insert the housing. Be sure that the incoming ra- 
diation is clear of the housing. 

xli) Put in the liquid helium shield, and the radiation 
shield. Put on the nose. Use GE varnish and aluminum foil 
to reduce openings in the baffles so that they will transmit 
only the bright white fringes. 

xiii) When the spectrometer is on the telescope, maximize 
the signal by tilting and rotating the dewar. 


(C) The Electronics 


The data-taking process is controlled electronically to 
ensure a smooth process. The operator only needs to turn the 
switch on. The spectrometer ..will then automatically take date 
in, process it, and stop at the end of the transform indicated 
by the operator. All the operator has to do in the whole obser- 
vation is to keep the astronomical object in the beam. Figure 
3-5 shows the block diagram of electronic and computer set 
up. The following are brief descriptions of the electronic 
parts Incorporated in the system. 

1 . Alignment Sensor 

The alignment sensor (figure. 3—6) is used to synchronise a 
Monsanto electronic counter, and the computer with the spectro- 
meter. The circuit is shown in figure- (3-6). The exit mask is 
continuously moving. When the exit mask is at its starting 
position, i.e. the first 255 slots are at the exit aperture, 
a light pulse goes through an alignment slot on the exit mask 
and is detected by a photo-cell on the other side of the exit 
mask. Two transistors amplify the output light curve and two 
IC chips change the light curve into an alignment pulse. This 
alignment pulse, through the drive unit, readies the counter 
for counting, readies the computer to accept data, and to turn 
on an indicator light showing that the system is taking data. 
One can adjust the starting position of the exit mask by ad- 
justing the intensity of the light. After 255 data points have 
been obtained, the exit mask is at its other end, and another 
alignment pulse turns off the counter and the indicator light. 

















The computer after taking 255 readings checks that the indi- 
cator light is off. This step ensures the synchronization 
of computer and spectrometer. 

2 . Drive Unit 

A drive unit (Fig. 3-7) is used to drive the entrance 
and exit masks of the spectrometer. It can be set to be 
adjusted by pulse streams at 100 Hz, 200 Hz, or 400 Hz. The 
unit performs the following functions: 

(a) It receives an alignment pulse from the sensor 
circuit ( 18 ) and generates a pulse to reset the counter for 
counting (17 ) . The pulse will also ready the computer for 
taking the data. 

(b) It drives the exit mask in a continuous mode at a 
displacement rate of one slot for every ill pulses(l 8 ). 

(c) After each set of ill pulses the unit instructs the 
the mini-computer to read the integrated signal off the counter 
and then reset the counter for the next data integration. 

(d) After 255 reset pulse the unit advances the en- 
trances mask by one slot by sending the entrance mask advance 
motor 40 pulses ( 11 ). 

3 . Preamplifier 

The circuit (Fig. 3-8) shows a transimpedance amplifier 
implemented with a Burr Brown operational amplifier. The cir- 
cuit has the advantage of high speed, low susceptibility to 
microphonics, and detector operation at constant voltage with 
a high resistance load resistor. 

Neglecting the voltage noise and current noise In the 























figure 3 — 8 . Basic dstsctor bias and prsamplifisr 
' * ” circuit. 
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first order approximation 3 the current through the detector. 
Ip, also goes through the load resistor R^ because the input 
Impedance of the operational amplifier can be taken to be very 
large, so 


V o " X D R F 


V. 


B 


Rp 


where V Q is the output voltage, 

V„ is the constant bias voltage, and 
R Is the load resistance 


(D) Laboratory Calibration 

■* 

1. Calibration of Post Optics 

Liquid nitrogen is used to calibrate the efficiency and 
sensitivity of the post optics. Liquid nitrogen is in a dewar 
with black paper along the wall to simulate a blackbody. 
Through a chopping device, the post optics will alternately 
see radiation from the liquid nitrogen and from the room. The 
difference between these two radiations gives the A.C. signal. 
D.C. measurements are obtained by putting liquid nitrogen 
directly in front of the dewar window. The following are the 
results of the calibration: 

Dewar profile: f H = 7*6 

f v = 10.3 
8.7t - ll.lv 


Wavelength region: 


Band-width: 
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2.4v 


Bias voltage: 

Load resistor at LHg. temperature: 
A.C. power: 

A.C. signal: 

A.C. noise: 


15u 

1.2 Mn 
-10 

9x10 watt. 
272.7 mv 
4 .4yv 


(NEP) A.C. detector: 
(NEP) A.C. system: 
A.C. responsivity : 
D.C, power: 

D.C. signal: 

D.C. noise: 


5.2xl0” 13 

z 

1.25xl0~ 12 H ^ 

I. 2 amp/watt. 
3.05x10“^ watt. 

II. 32V 
4yv 


(NEP) D.C. system: 
D.C. responsivity: 


1.47xl0" 12 
2.58 amp/watt. 


Background noise: 


NEP Blip /2 P BG hv 


=3.5xl0 -13 


where P BQ is the D.C. power 

v is the frequency that is assumed to be 10H . 

z 

The system is a factor of 3*57 away from background limited. 

2. Test of the Spectrometer 

A mercury vapor lamp with emission at 1.7y is used with 
the spectrometer to test the computer program. A PbS detector 
and an appropriate blocking filter isolates the 1.7y doublet 
of mercury. The slit width and length for each mask is 0.15 mm 
and 3.5 mm. Fig. (3-9) is the mercury vapor spectrum at 1.7w 
for the 1 x 255 mode. Fig. (3-10) is the mercury vapor spectrum 


i 


Hg Vapor emission line 
single entrance, 255 exit slit 
Integration time 4.0 min 


yv^v^vv/\/V A^Aj/yv, 




mercury vapor 1.7pm line using the 
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for 15 x 255 mode. Figure (3-10a) shows the eighth of a series 
of fifteen individual spectra. Figure (3-10b) shows an ave- 
rage of all fifteen spectra, and (3-10c) shows all 15 spectra, 
each spectrum being displaced vertically from the next. The 
diagonal pattern near the right-hand edge represents a dis- 
placement of the peak between successive spectra. This re- 
presents the actual shift in spectral range between adjacent 
spatial elements. 

Once it was clear that the instrument with the computer 
program worked properly in the 1.7y region, the spectrometer 
was tested with the cryogenically cooled, arsenic-doped sili- 
con detector. The transmission spectrum of polystyrene with 
a soldering iron as the source was obtained. The shap e of the 
filter profile and the transmission spectrum of polystyrene 
showed that the instrument worked in the lOy region. 

The wavelength calibration of the spectrometer is obtained 
by comparing the polystyrene transmission spectrum obtained by 
the spectrometer and the spectrum obtained by a Perkin-Elmen 
monochromator. The wavelength calibration is then checked 
against the moon spectra at 9.5v where the atmosphere has 
strong absorption features. 
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Figure 3-10 » 

Calibration spectra obtained 
for the mercury vapor emi- 
ssion lines at 1.69vm and 
l.Tlym: (a) The eighth of a 
series of fifteen individual 
spectra obtained; 


of all fifteen 


(c) A display of the fifteen 
spectra, each spectrum being 
displaced vertically from 
the* next. The diagonal pat- 
tern near the right-hand 
edge represents a displace- 
ment of the peak between 
successive spectra. This re 
presents the actual shift in 
spectral range between ad- 
jacent spatial elements. 








CHAPTER IY 


COMPARISONS BETWEEN FOURIER TRANSFORM AND 
HADAMARD TRANSFORM SPECTROSCOPY 


Since the Michelson interferometer spectrometer (MIS) and 
Hadamard transform spectrometer (HTS) both have the multiplex 
advantage and the advantage of large through-put , there are a 
number of comparisons between them In the literature. In this 
chapter the comparisons are carried out in four different as- 
pects: mathematically, computationally, optically and mecha- 

nically. 

(A) Mathematically . 

Fourier transforms and Hadamard transforms can be viewed 
as using two different weighing schemes. Appendix A gives the 
mathematical analysis of the coding error for Fourier transform. 

Let £ be the path difference in a two beam interferometer. 
P(v) is the power at wave number v (i.e. the spectral density 
function), v here is taken to be the inverse of wavelength, 
then S(£), the power received for path difference £, is: 

S(£) = I P(v) cos 2 (2ir£v)dv . A-l 

0 

= l/2P rt + l/2/“p(v) cos ( 4 tt£v ) dv A-2 

In the Fourier transform each spectral element can be 
viewed as having a phase modulation by a cosine squared term. 

Its argument depends on the stepped distance and the particular 
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* .* 

wavelength. In the Kadamard transform each spectral element 
is modulated by a step function of 1 and 0 for a S-matrix 
coding. 

One can ask oneself whether these modulating or encoding 
schemes are equally efficient? 

In table 2-1 we stated that for large N the multiplex ad- 
vantage, or the efficiency, of the H-matrix coding of elements 
1 and -1 Is vTT, the S matrix coding of element 1 and 0 is ^ 

C. > 

and for a single detector MIS the Fourier coding is /N/% . 

This last figure is based on calculations shown in Appendix A. 
The other two values were described by Sloane et al (1968). 

The H-matrix is an orthornormal matrix. The element -1 means 
"subtract" the radiation, while the +1 element means "add” the 
radiation. No radiation Is wasted and thus the coding scheme 
has the highest efficiency. For S-matrix coding half the 
slits are open, 1, letting light pass through; and half the 
slits are 0, blocking the light. Half the radiation is there- 
fore wasted each time and one can intuitively see why the 
efficiency of the S-matrix is only half that of the K-matrix. 

Although the S-matrix is less efficient, it has the 1m- 

til 

portant advantage that it is cyclic; that is, the (i+1) 

til 

column of the S-matrix is obtained by shifting the i column 
cyclically one place downwards. Instead of constructing a 
mask of N 2 siits for N spectral elements, one constructs only 
one mask with 2N-1 slits. Such a mask has two advantages. 
First, the cost of mask construction is reduced by and 

the design' of the advancd mechanism is considerably simplified, 
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since the total weight of the mask also decreases as *\,N/2. 
Secondly, it can be self-supporting and therefore permits the 
construction of a spectrometer which requires no transmission 
materials. In operation the mask is stepped one slit width - 
along the length of the mask- for each successive encoding po- 
sition. 

In the Fourier case, equation (A-2) of Appendix A 

SCO = 1/2P + l/2/"P(v)cos(iJir!;v)dv (A-2) 

o 0 

shows that half the power goes into 1/2P , the first term on 
the right hand side, which is not modulated at all. This re- 
duces the efficiency by a factor of two as in the S-matrix case. 
The other half of the power is modulated by a cosine term. 

Cosine modulation gives a factor 1//2 because cosines functions 
do not form an orthronormal set themselves. The total Fourier 
modulation efficiency is therefore 1/2^/W/^. Mathematically, 
the S-matrix is a factor of /2 better in SNR than the Fourier 
transform. The true Hadamard code, which cannot be realized 
experimentally, as yet, is a factor of /8 better that the 
Fourier code. 

* \ 

(B) Computer Requirements 

Both MIS and HTS require a digital computer to decode the 
data. However, in the HTS case this reduces to nothing more 
than a series of additions and subtractions. Hence, as much 
as an order of magnitude in computer time can be gained over 
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the Fourier decoding procedure required by MIS (Decker, 1971). 
In addition, HTS do not have the large zero path-length spike 
which is characteristic of MIS, and hence can be operated 
with a substantially loader dynamic range. 


(C) Optics 

HTS attempts to "liberate” the grating Instrument from 
its inferior position and offers a possibility to convert a 
conventional scanning spectrometer into a multiplex instrument 
at a moderate cost. However, it is also the grating and optics 
that limit the capabilities of HTS. 

1. Resolution 

The resolution of a grating instrument is 


R 


X 

AX 

mN 


(4-1) 

(4-2) 


where m is the order of diffraction. N is the total number of 
rulings in the grating. W is the width of the grating, d is 
the separation of lines. 

The MIS introduces variable path differences between two 
interfering beams. The resolution is determined by the maxi- 
mum permissible path difference in the interferogram. 


Av “ 


1 

2x 


(4-3) 


1 

n 


(4-H) 


where x is the maximum path difference between interfering 
beams . 

E is the displacement of one of the two mirrors from 
the white light fringe position 9 and 
Av is the increment in wavenumber. 


Since 

R 


A 

AX 


1 

AAv 


(Wertz a, P.5) 


2x 

X 


(4-5) 


A MIS can have a resolutuion as high as *vlO 


2. Slit Width 

In the HTS, the minimum usable slit width is determined 
by the diffraction pattern 


w * ^ F (4-6) 

where P is the focal length of the imaging mirror. 

It has been argued that a boxcar profile is a poor match 
to a sine diffraction pattern(Hirschfeld s eb al, 1973 9 Mertz, 
1976 b). In the. presence of diffraction, the transmission at 
each point of the mask will be a complex function of the spec- 
tral distribution, the mask position, and the relative width 
of the nearby transparent and opaque slits. 

One way to correct this is to make the slits, wider than 
the diffraction limit, allowing the mask's transmission to 


approach the geometric optical limit. This way, however, not 
only the resolution of the instrument is reduced, but the to- 
tal number N of spectral elements that can be observed simul- 
taneously, also decreases. 

There is another way to look at the same problem. The 
grating is an operator that changes the frequency domain into 
the spatial domain, so that intensity as a function of fre- 
quency, after passage by the grating, becomes a. function of po- 
sition. With diffraction effects included, the intensity has 
a new functional dependence on distance. The Hadamarcl mask 
code and the subsequent decoding process just translate this 
spatial distribution function back into its spectral domain. 
Therefore, an intensity pattern which is complicated by the 
diffraction pattern, after the Hadamard coding and decoding, 
should still show up the same intensity pattern. Since one 
knows the diffraction pattern for a given optical system, the 

diffraction effect on the spectral intensity distribution can 
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be computed and corrected, so that the sine diffraction effect 
would not make the slit width any larger than the diffraction 
limit. One point, however, should still be noted. The cor- 
rection that would have to be applied is wavelength dependent 
because diffraction is wavelength dependent . 

A similar problem exists in MIS (Stewart P.296) because 
the moveable mirror must be stopped when its maximum displace- 
ment is reached. 

The sidelobes in the interferometric case is of the form 


S(wO = — — — «««» (i|-7) 

to -10 1 

(Stewart P.295) 

which has considerably stronger side lobes than the diffraction 
pattern. Here, T is the time for the mirror to travel from one 
end to the other, and to is the central frequency. Various 
schemes of apodization have been introduced to compensate for 
the side lobes in the interferometric case. 

3 . Multiplex Number 

The total number of elements N that can be observed simul- 
taneously yields the multiplex advantages. In HTS, the total 
aperture size is limited by aberrations largely due to off axis 
radiations. The aperture width is (Hirschfeld, 1973) 

w = P • S (4-8) 

w 

where F is the instrument focal length and S /2 is a factor of 

w 

the order of O.OS' 1 ' 0.], that describes how far off axis one 
can go before the aberration pushes the individual slit width 
up to the point where no further gain in W is possible. 

Since the minimum slit width, determined by diffraction ef- 
fects, is ^l,22Af, the total number N is 

FS 

N * — £ (4-9) 

1. 22Af 

and is of the order of 10 . 


For our HTS at Cornell, we have 
F = 49.5 
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S^O.08 (assume ^ medium value) 

X= lO.Oym 

f= 7.5 

1H250 

Actually, however, Mertz and Flamand (1976) suggested that 
values >>0.1 may be realized in practice. 

The MIS has a very large effective value cf N. This is 
its main gain. The total wave number range observed in MIS is 

v max " v min = h 

1 

TZs 

where a is the step size of the mirror, o is the resolution in 

g 

wave number. N for the MIS can be 10°. 

Spectral Range 

The MIS also has a broad free spectral range. Its range 
is limited by the beam solitter efficiency which usually varies 
approximately as the cosine of the wavelength. HTS free spec- 
tral range is usually about one grating order. Its free spec- 
tral range can be increased by using order sorting, but this 
increases the technical difficulty.. Although the HTS has a 
smaller spectral free range, it can be set to recover only 
those spectral bands of particular interest throughout any • 
spectral regions (Decker, 1971). This is impractical with a 
MIS. 
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5. Throughput 

The throughput is defined in section 1-A as the product 
of aperture A and angular acceptance ft 
E = Aft 

For a MIS it can be shown that ftR=2v and 


e mis 



2tt 

R 


where R is the resolution of the Instrument. A Is the area 

m 

of the interferometer mirror. Typically, A^l cm 2 , and for 
R^IO 3 , E mis ^6 x 10~ 3 . 

For the HTS, from equation (3-1) one obtains 


e hts 



w^h\ 

2 ' 
F* 


where h and w are slit height and width and F is the focal 
length. For a double multiplexed spectrometer, the throughput 
of MIS and DHTS are about the same, of the order 10~ 2 cm 2 . 

For the same throughput, the HTS may have a worse system trans- 
mission because the HTS requires a dedispersing process. 

DHTS has the additional advantage that one can construct 
a one dimensional picture of the source. 


(D) Mechanical Requirements 

The HTS mask can be made self-supporting hence no beam- 

; 

splitters or transmission optics, are required. Furthermore, 
in the MIS, construction tolerances usually involve dimensions 


and motions that have to be maintained to within fractions of 
wavelengths. For a HTS, the corresponding tolerances are frac- 
tions of a slit width, and these tolerances are normally two 
orders of magnitude more relaxed, so that this Instrument will 
be more suitable for rugged applications and less costly. 

It is clear from the above comparison that the MIS has 

the advantages of highest resolution, very large multiplex 

number and free spectral range. The HTS has the mechanical 

advantages and computational advantages for large N. The HTS 

•3 

can have on the order of 10 J spectral elements and a resolu- 
tion sufficient to resolve the rotational lines of many mole- 
cules. For most IR astronomical observations this will be 
sufficient. Furthermore, its potential for modifying the 
existing sacnning spectrometer at a moderate cost make this a 
very worthwhile field for further study. 


CHAPTER V 


PROGRAMMING FOR HADAMARD TRANSFORM 
SPECTRAL DATA REDUCTION 


An 8K Computer Automation minicomputer model L.S.I. or 
model Alpha-l6 can be used to interface with the output of a 
Monsanto scalar counter which digitized the output of the de- 
tector used with the Hadamard transform spectrometer- The com- 
puter processes each data point as soon as it receives it and 
when the data gathering run is completed, the computed spec- 
trum is also ready within a fraction of a second. The final 
spectrum can be displayed on a cathode ray tube for quick vi- 
sualization, or printed on paper by a teletype machine for 
more detailed analysis. Also, it can be stored on paper tape 
for future use. 

There are two inverse transformation programs: 1 x 255 

for the single entrance slit and 255 exit slit Hadamard trans- 
form spectrometer, HTS, and 15 x 255 for the 15 entrance slit 
and 255 exit slit instrument, DHTS. 

(A) HTS Program 

The inverse HTS program processes the raw data obtained 
by the combination of a single entrance slit and 255 exit slit. 
This program is in double precision format. Two areas in the 
memory are reserved by the program to store the final spectrum. 
The firal spectrum can be stored in either the plus beam or the 


75 


76 



minus beam area. The plus and minus beams are arbitrarily 
named. 

Data can come in at a rate of 10 data points per sec,, 

5 data points per sec., or 2.5 data points per sec., depending 
on how the clock driving the spectrometer is set. Since a com- 
plete pass has 255 points, each pass takes 25.5 sec., 51 sec., 
or 102 sec., depending on the data rate.- Each pass yields one 
spectrum. One can take as many passes as one wants until one 
is satisfied with the SNR of the spectrum. 

The whole program (Appendix B) is linked by the following 
subprograms: COMMAND, TRANSFORM, INPUT/OUTPUT, READ, CLEAR, 

PUNCH, GRAPH, DISPLAY, MATHMATICAL PACKAGE and MASK. The func- 
tion of each subprogram is described briefly in the following 
sections. 

1. COMMAND: This subprogram performs two functions. 

(a) It commands the computer to do one of the following func- 
tions: TRANSFORMATION, CLEAR, READ, PUNCH, GRAPH or DISPLAY. 

(b) If two distinct spectra are stored in two different beam 
areas, the COMMAND program can take the difference and ratio 
of the two spectra. 

2, INVERSE TRANSFORM: This is the most important subprogram 
in the whole program. It will take data points from the counter, 
transform them and enter them into either the plus or minus beam 
areas, or it will read the data points from the paper tape into 
one beam area, transform them and enter them into the other 
beam area. This program is called the inverse transform, since 
it inverts the trar.s format ion performed by the coding mask, and 
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yields a spectrum. The program also performs the Hadamard 
transform, i.e. transforms the spectrum back to raw data, from 
either input . 


The algorithm is based on the following idea: Let ^ be 

the j tJl spectral element and w^ be the xfeight of the j th ele- 


mask. w.. eauals 1 for transmitted radiation 


ment of the i^* 1 
and 0 for blocked radiation. Each measurement then has a 
value 


255 

z s 

i=l 


ij*j + 


where v i is the random detector noise satisfying the properties 
mentioned in Section II-B. S_ is the 255 x 255 matrix, is 
the i^ data point entered into the computer. The computer's 
job is to decode to reconstruct the original spectral values 
. Therefore 


255 

I 

i=l 


S ji "i 


* _1 
where is the unbiased estimate of ^ and S is the inverse 

of the S matrix. 

According to the relation 


S" 1 = |(2S - J) 

one obtains S” 1 by keeping all +l's in the S-matrix and re- 
placing all 0’s by -1. The matrix obtained in this way is the 


2 

Inverse matrix of S except for a constant factor 255 ** w ^ ic ^ 

only gives a different normalization. 

To reconstruct the spectral values ^ one needs to add or 

subtract each measured value to different bins, according 

to whether S ” 1 is plus or minus. Fig. (5-1) is a flow chart 

for the 1 x 255 transform program. 

By a Hadamard transform we mean a program that transforms 

* 

the spectrum back to its raw data „ This procedure is useful 
because by inspecting the raw data display which usually appears 
quite smooth, any bad data point can be easily identified, 
and for example, replaced by the average of its two adjacent 
data points. This procedure will improve the final spectrum. 

The Hadamard transform turns out to be extremely easy.. 

All one has to do is to change one statement in the inverse 

-1 

transform program. When is -1, instead of negating the 
data, one just sets it to zero. 

Data points are taken both with the exit mask moving in 
a forward and in a reverse direction. The inverse transform 
program takes care that when the exit mask moves in the for- 
ward direction, the spectrum is transformed into the plus beam 
area. When the exit mask moves in the reverse direction, the 
■final spectrum is stored in the minus area. Not adding the 

The notation here may be a bit confusing. We use the S ^ to 
transform the raw data into an intensity spectrum. The in- 
verse transform uses S to transform the intensity spectrum 
back into raw data. 


























forward and backward spectrum eliminates a degradation of the 
final spectrum due to any asymmetry between the data taking for 
different directions of motion of the mask. 

3. IN?UT/OUTPUT : This program links up the computer, the 

teletype and high speed reader. It consists of the following 
functions: Keyboard Input, Paper Tape Input, Output to Tele- 

type, Output Text from Buffer, Output Floating Point Number, 
Wait. for Execute Signal, Command Error Exit, Carriage Return- 
Line Feed. 

REAP, CLEAR, PUNCH : This program reads the. spectrum 

from the paper tape Into either plus or minus beam areas for 
further manipution, or punches the spectrum out from the beam 
area; it also can clear the beam area. The speed of teletype 
for reading Is 100 words per sec. Hence it takes about 12 min. 
to read the spectrum. Punching has the same rate. 

5. GRAPH: This program plots the graph on teletype paper, 

with its numerical value in floating point format. This pro- 
cedure offers one the chance to Inspect the spectrum in detail 
if it is needed. It takes about fifteen to twenty minutes to 
finish a spectrum, depending on the complexity of the spectrum. 

6. CRT : A cathode ray display is interfaced with the out- 

put of the computer. It takes about 2 sec. to display the 
spectrum, with a factor of 5 higher resolution than the graph 
printed by the teletype on paper. One can display the spec- 
trum at the end of any pass to see how good it is. 

7. MATHMATICAL PACKAGE : This package is supplied by the 

Computer Automation library tape, with a little modification 
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from our own on its double precision part. 

8. MASK : This part contains the matrix, the 255 ele- 

ments exhibited in the first mask position plus an additional 
25*1 elements representing- the further • cycling of this mask. 

(B) DHTS 

The doubly encoded Hadamard transform program processes 
the data obtained by the various combinations of fifteen en- 
trance slots and 255 exit slots. It is a single precision 
program. It can accept data at a rate of 5 data points per 
sec. and 2.5 data points per sec. on'i* because it takes a 
longer time to process each data p ’ ..ne whole transform, 

takes about 1^ minutes. The final am consists of 15 se- 

parate spectra, representing a one-dimensional color picture 
across the spectrometer entrance aperture. Each separate pec- 
trum contains 255 spectral elements. The program can co-add 
all fifteen separate spectra yielding a sum spectrum with im- 
proved SNR. 

The program consists of the following subprograms: 

COMMAND, TRANSFORM, DATA, INPUT/OUTPUT, READ, CLEAR, PUNCH, 
GRAPH, DISPLAY, ENTRANCE MASK, EXIT , ppendix C). C 

Since most of the subprograms ar** TVv-o: ie same function 

as their counterpart in the 1 x 255 c pt written in 

single precision format, their descr J — l1 not be repeated 

here. Only the TRANSFORM, DATA, and programs will be 

discussed because they are different from those in the 1 x 255 


scheme. 
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1. INVERSE TRANSFORM : The program can accept data eigher 

from the counter or from paper tape. In order to eliminate 
any asymmetry due tc the different directions of motion of the 
exit mask, the computer will accept data only when the mask is 
moving in a given direction, either forward or backward, de- 
pending on the operator. If one wants to save time, one can 
still choose the mode in which the computer will accept data 
in both directions of mask motion. Hence, the program provides 
six modes for operation: accepting data from the counter, in 

the forward, backward, or both directions, and accepting data 
stored on the paper tape, in the forward, backward and both 
directions. Figure (5-2) shows the flow chart for 15 x 255 
program. 

Let both the entrance and the exit masks be linear arrays 
encoded by Reed-Muller codes. Then the matrix of spatial- 
spectral elements ifi is related to the matrix of measurements 
n by 

s_ jp S - jn 

To obtain the spatial-spectral Information about the viewed 
scene we solve this equation by premultiplying the data matrix 
by s" 1 and postmultiplying by S -1 

. _ -1 0-1 
jjj - s_ ji S 

Now consider the element It Is multiplied only by 

-1 

elements of the first column of s_ ; and in turn it multiplies 
only the elements of the first row of S . To a given spectral- 






















spatial element tp^ , It therefore contributes an amount . 

—l -1 

s 4 -Ti, n S ,4 . But the elements s*. and s . . all have values, 
il 11 lj il ij 

either +1 or -1, and the result is that each element ip^ of 
the matrix receives a contribution or "H-n ^ rom the 

reading, ri^. 

This procedure is generally valid. Any reading tk ^ will 
.make additive contributions that can only have values or 

-n^ to each element ip^ of the ip matrix. 

For real time decoding we therefore need the following: 

(a) A memory that consists of bins containing the con- 
tributions to the elements ip . accumulated up to any given 

1 J 

time t in the cycle of measurements. For a device that can 
resolve m spatial and n spectral elements, this memory re- 
ruires mn bins and of the order of mn memory words. 

(b) For each acquired reading we perform a series of 
additions of values either +n k£ or -n ki , one to each of the 
ip.y memory bins. But before that can be done, we need to de- 
cide on the assignment of + or - needed for a given bln. 

This is done in the following way. 

We store the sequence .of + and - signs in one column of 

s” 1 and in one row of S -1 and in one cycled permutation of 

each of these vectors. Let us designate these signs by their 

“1 —1 

positions in these two vectors, as and , respectively, 

i=l,... m; j=l, ... n. (Since each of these sequences is cyclic 
it can, respectively, be brought into its kth and ith cycling 
position after a measurement ri^). The elements of the two 
vectors then are multiplied in all possible combinations to 
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give a matrix having mn elements. 




s'V 1 

s i 


i=l* • * * m; j— !• • • • n 


Each element z^ is either + or - depending only on whether the 
signs s" 1 and Sj 1 are similar or dissimilar for a particular 
combination of i and J values. 

The additions +n k£ or -n^ to the bins are made as 
successive elements, Z^ are computed, so. that the elements 
Z^ need never be stored. Figure (5-3) shows the relation of 
Z^ to a superarray containign the set of all elements that 
are constructed at various stages of the computation. 

When only a restricted number of spectral elements are 
of interest, we need to compute elements representing only 
selected j values. This might be useful, for example, if only 
certain atmospheric C0 2 absorption lines needed to be studied, 
and the spectral elements between were of lesser interest. In 
that case only 2^+1- 1 £,+^-1 e ^- emen ts corresponding to given j 
values need to be used, and the computing time decreases as 
p/n, where n is the total number of available spectral elements, 
and p is the actual number of interest. 

—l —1 

One starts with the inverse of the codes, s and S , for 
the entrance and exit masks stored in the computer. One stores 
509 elements of the exit mask, i.e., the 255 elements exhi- 
bited in the first mask position plus an additional 25^ ele- 
ments representing the further cycling of this mask. Simi- 
larly one stores twenty-nine elements for the entrance mask. 
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Figure 5-3 . Matrices generated by the computer during the 
reduction of the spectral data. 
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representing the first fifteen elements used, plus the further 
cycling of fourteen elements. 

For each reading one essentially makes use of the ma- 
trix (figure 5-3) making use of elements k to k+14 of the stored 

entrance code and elements i to £+25*1 of the stored exit code. 

—1 -1 

This matrix consists of + and - signs. When s^ and have 
the same sign, both being + or both being -, the matrix posi- 
tion is assigned a + sign, and the reading is added to 

-1 

the accumulatively stored value of . If the elements 

— 1 

and have dissimilar signs, a - sign is assigned to ij and 
the reading is subtracted from the stored- values. This 
whole process takes *100 msec, and is carried out yrhlle the 
succeeding intensity measurement is being made. 

In prat ice, we start with the first spatial element, i=l 
and add or subtract the contributions to all the values, 
successively going from j=l to 5=255* We then repeat this pro- 
cedure for i values going from 2 to 15* This whole procedure 
is carried out while the exit musk is moving from posiiton Z 
to Jl+1 depending on whether the exit mask is moving forward 
or back. The entire process is then repeated for the next 
reading • When the exit mask reaches its 255th posi- 

tion, J l remains unchanged, but the entrance mask moves from 
the position k to k+1. ■ 

For odd values of k, the exit mask moves in the direction 
of increasing l values, and for even values of k, it moves to- 
ward decreasing values. In short, the exit mask moves back 
and forth as readings are taken. After the entrance mask has 
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moved, through all. Its fifteen posit ions * .and the total- 3825. 
readings have been taken and' added onto or subtracted from the 
elements, the run is completed. 

2. DATA ; Instead of processing a data point immediately 
as it comes in, this program stores the data in the memory, so 
one can display the raw data points first, correct them if there 
are any obvious bad points, and then transform them. This pro- 
' gram serves the same function as the inverse transform in the 
1 x 255 system. 


3. DISPLAY : The display program allows one to display the 

information in a number of different ways: 

(a) One can call for the spectrum corresponding to any 
one of the entrance slit positions and display it individually. 

(b) One can display the sura of the different spectra. 

In order to do this, one has to take into account that the 
spectrum for a given entrance slit is .displaced by one spec- 
tral position from adjacent entrance slit position. In other 

-1 

words, the wavelength for element corresponds to the wave- 
length for element j+1* because 0 f the slightly displaced 

light paths through the spectrometer. 

(c) Finally one can display all fifteen of these spectra 

simultaneously, with the zero baseline of each spectrum- vert I- 

* 

eally displaced from the next one. While this format is some- 
what crowded, it does allow a quick; comparison of the indivi- 
dual spectra. 
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(0) Correction Program ; 

This program is shown in Appendix D written in BASIC lan- 
guage. It corrects the error introduced by the imperfect mask. 
Instead of correcting the mask which in practice is not pos- 
sible, the program corrects the final spectrum. That is much 
easier. 

As seen in section II-D, for any spectral line I Q , the 

i 

distorted spectrum shows a line I 0 =I 0 (l-e), two positive blips 
with amplitude (1/2) (el^/l-e) adjacent to the line on both 
sides, and two negative blips with same amplitude, i.e. (e/2) 
(I^/l-e) at 24, 25 elements to the left. The correction pro- 

t 

gram takes the intensity of every element, I Q , multiplies it 
by e/2, adds it to the elements 24 and 25 positions to the 
left, and subtracts it from the two adjacent elements, one on 
each side of the line. The final spectrum is complete except 
for a different normalization factor. This is a linearized 
correction procedure valid only for small values of e, e<<l. 


CHAPTER VI 

ASTRONOMICAL OBSERVATION 
(A) Correction Procedure 

The correction of the spectra for telluric absorption and 
for instrumental response is a critical procedure. The correc- 
tion is carried out by comparing the source spectrum (either 
stellar or planetary) with a lunar or solar spectrum which is 
taken on the same day at an airmass as close as possible, to 
the star. The following procedures are used in the data re- 
duction: 

(1) Correct the raw star spectrum and lunar spectrum (or 
sun) for negative dip due to the imperfect mask as described 
in section II-D. 

(2) Correct for different entrance slit width if neces- 
sary because different slit widths will give different reso- 
lution. 

(3) Most of the astronomical infrared sources to be ob- 
served are weak sources, hence a positive offset is always 
added to the signal to prevent the signal becoming negative. 
(The ..electronics are confused by negative signals.) The cor- 
rection procedure shown previously also take out this offset. 
One can take out the offset from the raw spectrum if one knows 
how large the offset Is. Another way to correct this is by 
subs tract ing a constant intensity from the stellar spectrum 
until the ratio of the stellar spectrum over the Moon spectrum. 


around" any large telluric absorption reigori, such, as the ozone 

■v ■ 

band, is optimally corrected. By optimal correction we mean 
that the atmospheric feature appears as neither a positive, nor 
a negative hand. The Moon is a strong Infrared signal and does 
not require an offset, so the result can be used as a calibra- 
tion for base line. 

CM) Align the stellar spectrum and the lunar spectrum by the 
telluric absorption feature. The final stellar spectrum is ob- 
tained by taking the ratio of the stellar spectrum and the lu- 
nar spectrum and multiplying it by the black body temperature ■ 
of the Moon. The lunar temperature is obtained by noting the 
phase angle of the lunar east limb where it has usually been 
observed, and extrapolating the temperature from the value gi- 
ven by Linsky (1973). This procedure assumes the lunar in- 
frared emissivity at 8-1 4 ym as unity, which is not true. Murcray 
et al* s (1970) results for the lunar emissivity of 8-l4ym re- 
gion are shown in figure (6-1). The observation was done from 
a balloon. The strong feature centered at 9.6|im is a result 
of telluric ozone absorption. Murcray* s result has not been ' 
used in our analysis because in the region to be discussed, 8.5vm- 
l4ym, the Moon’s emissivity is about constant except for the 
ozone absorption. ■ 

(B) Observation of ot-Orionis - 

• The observations of ct-orionis were carried out with the 50” 
infrared telescope at Kitt Peak National Observatory, Arizona, 
in May 1974. The beam size is 18 sec. x 47 sec. The Kitt Peak 
50” telescope bolometer was used with the Hadamard spectro- 
meter. The dewer has a band pass of 8 - l4ym. The spec- 







trometer operated in the 8-llym region with a resolution X/AX 
around 500. The chopping frequency was 10 cycles per second. 

Three runs were taken. Each run consisted of 10 scans of 
the sources. Two lunar spectra were taken on the same day 
for correction purposes. For the lunar temperature we used 
383°K. Figure (6-2) shows the raw spectra of a-orionis and 
the Moon. The a-orionis spectrum is the sum of two indepen- 
dent runs. . 

Figure (6-3) shows the ratio spectrum of a-orionis corrected 
for lunar temperature. Except for the region immediately around 
the ozone band all the telluric absorption features are gone. 

The region between 9.35um to 9-7vm is unreliable because the 
ozone band has a very low transmission. 

a-orionis is a late type super-giant with temperature 
around 3000°K. A stellar continuum corresponding to a 3000°K 
blackbody is also shown in figure(6-3) s normalized to arbitrary 
units. The broad .emission feature around 10ym, which is due 
to silicate emission, is clearly shown in the spectrum. This 
lOym ’silicate' emission feature of a-orionis has been previous- 
ly discussed by others. 

Woolf and Ney (1969) renormalized Gillett et al’s spectra 
(1968) of a-orionis and interpreted the 10ym emission as coming 
from circumstellar dust which absorbs starlight and reradiates 
at infrared wavelengths. Since the emission is far more sharply 
peaked than a black body, the wavelength dependence of the 
emission probably closely mimics the wavelength dependence of 
the opacity of the material. In the same paper, Woolf and Ney 


and the 






also proposed that it :5 s silicate grains which one is observing 
in the circumstellar dust cloud. This suggestion has been 

strengthened by high resolution spectra of Gamown et al (1972) 

— 1 —1 

with resolution X/AX-250 from 850 cm to 1100 cm . The emi- 
ssivity of silicate grains should have a second peak near 20ijm. 
This second emission feature was observed by Low and Swamy (1970) 
in narrow-band photometry of a-orionis. Another supporting 
piece of evidence for the silicate model of the a-orionis dust 
cloud comes from observations of silicon monoxide (SiO) . 

Silicon monoxide is expected to be among the most abundant 
molecules present in the atmospheres of cool stars of normal 
composition. It is also a reagent in the condensation mecha- 
nism thought to produce circumstellar silicate grains (Mass 
efc al, 1970). The presence of silicon monoxide In a-orionis 
was confirmed by Cuduback et_ al (1971) • They observed SiO ab- 
sorption features around 4um. 

The silicates are expected to form from the material 
ejected by cool stars. Gilman (1969) calculated what solids 
would condense from gas of stellar composition as it moved 
away from a star and cooled. This is critically dependent on . 
the ratio of oxygen to carbon in the gas. These elements 
first combine to form carbon monoxide ; the subsequent develop- 
ment depends on which of the two is left over when all the other 
has been used up in this way. If carbon predominates s graphite 
will be the principal condensate , or under certain circumstances 
silicon carbide. If oxygen wins , the grains which should form 
are the silicates of calcium, magnesium, aluminium and iron; 


combinations of these are responsible for the 10 and 20pm spec- 
tral features. Aluminium and calcium silicate may be rare be- 
cause of the low cosmic abundance of aluminium and calcium. 
Woolf and Ney (1969) expected magnesium silicate, MgSiO^, with 
some iron silicate, FeSiOg, to be most abundant. Por recent 
work on dust grains one can refer to Salpeter's paper (197*1 a) 
on the theory of nucleation and dust grains in carbon-rich 
stellar atmosphere and his paper (1974 b) on formation and flow 
.of dust grains in cool stellar atmospheres. 

Laboratory spectra exist for silicate absorption features 
(Day, 1974). Any fine feature of the astronomically observed 
silicate emission may be washed out by different particle sizes 
and shapes, uncertainty in temperature and mixture of composi- 
tion. Gammon et al (1972) examined the excess in XY Seg and 
0 Get and concluded that the type of silicates involved are 
basic rather than acidic. Day (1974) synthesized the amorphous 
magnesium silicate and obtained an absorption band quite si- 
milar to the material causing the interstellar 10pm absorption 
feature. He suggested that the exi stance of disordered struc- 
tures seems a more reasonable expectation than crystalline 
terrestrial-type minerals. Por the momnet , the nature of the 
silicates is certainly at an unsettled stage. 

Penman (1976) had measured the middle infrared reflectivi- 
ties of five silicate minerals, and used Kramers-Konig analysis 
to obtain the optical constants of the samples. He then used 
the optical constants in Mie computations of the absorption 
properties of very small mineral grains. The final absorption 
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cross-section spectra compared to the observed lOym silicate 
of- the source W3/IRS5. All the calculated spectra have sharper 
features than the astronomical features due to the application 
of Mie theory. However, the hydratedsllicate, Chloritite 
(hydrated Mg/Pe/Al silicate) and serpenlenite (hydrated Mg/Pe 
silicate) fits the observed astronomical positions correctly. 
They fall almost exactly at the center of the astronomical 
features. ' : 

To the author's knowledge only two spectra of a-orionis 
in 8- 14 yin exist in the literature. Gillett et al (1968) 
(figure 6-4a) obtained results in the wavelength region from 
2.8 to l4ym, with a resolution X/AX=50. Treffers and Cohen 
(1973) (figure 6-4 b) obtained a spectrum from 8-l4ym, and lh 
the 20ym region with resolution 1000. Gillett f s and Treffers 
and Cohen's spectra are shown in figure (6-4). Compared with 
our result, all show the lOyra emission feature if Gillett* s 
black body curve is lowered instead of being drawn tangent to 
the observed data at 10ym. Our spectrum shows a rather rapid 
dip at wavelengths beyond 10ym while Treffers and Cohen show 
a slower nearly constant decline. Further high resolution 
observations should clear this matter up. Our own instrument 
now operates at a resolution similar to that of Treffers and • 
Cohen, and if used on a telescope as large as the 120 inch 
Lick Observatory reflector that they used, sufficiently high 
signal to noise ratios should be obtained. 
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(a) Low resolution spectrum taken by Gillett 
et al(1969). Different symbols represent 
spectra taken oh different nights, (b) High 
resolution spectrum taken by Treffers and 

Cohen (1973). 







(C) Observations of Jupiter 


The observations of Jupiter were also carried out with 
the 50" infrared telescope at Kitt Peak. The beam size is 
18" The Kitt Peak bolometer assigned to the 50” tele- 

scope with band pass 8-l4ym was used with the Hadamard spec- 
trometer. For our Jupiter observations, the spectrometer 
operated in the 10.8-13.4ym region with an effective resolu- 
tion X/AX around 250. The chopping frequency was 10 cycles, 
per second. 

Two runs were taken. Each run consisted of twelve scans 
of the whole Jovian disk. Two lunar spectra were taken on the 
same day for correction purposes. For the lunar temperature 
we used a value of 383°K. Figure (6-5) shows the raw spectrum 
of Jupiter and the Moon. Both the Jovian and lunar spectra 
are the sums of two independent runs. Figure (6-6) shows the 
ratio spectrum of Jupiter corrected for lunar temperature. 

The arrows labeled by HjO show the position of telluric water 
vapor features. Most of the telluric features have been can- 
celled properly. 

Jupiter is covered by clouds which in the visible part of 
the spectrum are seen from earth. Current models show three 
distinct cloud layers (figure 6-7, Ingersoll). The lowest 
layer is water ice, with maximum density at about 270°K. The 
middle cloud is solid ammonium hydrosulfide '(M^SH) at about 
200°K. Lewis and pr inn (1970) suggested that ultraviolet ra- 
diation from 2200 to 2700 1 is not absorbed by He j CH^ and 
absorbed a little by NH^. Therefore, the radiation in this 













region may reach the ammonium hydrogen cloud and photolyze hy- 
drogen sulfide there into hydrogen polysulfides (H 5 S ) , ele- 

^ X 
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mental sulfur, and ammonium polysulfides 
these species are yellow,* orange, or brown and may explain the 
color of zones. However, Sagan and Salpeter (1976) suggested 
that even under the most optimistic assumption that every HgS 
photo dissociation event leads to polymeries, the implied op- 
tical depth falls short by two orders of magnitude from mat- 
ching the observed values. Moreover, pure polymeric sulfur 
fits the observed optical properties of the Jovian red chro- 
mophores only poorly (Rages and Sagan, 1977 ). The upper cloud 
is solid ammonia at around 150°K. Solid ammonia is whitish and 
probably forms the white zones on the Jovian disk. The color 
of the Great Red Spot may be due to high altitude ultraviolet 
photolysis of phosphine (PH^) into Pglfy and amorphous red phos- 
phorus. The total depth of the cloud system is about 70 Km, and 
the pressure range probably runs from about 0.5 bar at the top 
to 4.5 bar at the cloud base. It is the cloud tops and above 
where the 10um infrared radiation originates. Our spectrum 
measures , a ..color temperature of .l35°K which is consistent with 
Ingersoll’s picture. The radiation should come from the cloud 
tops because the Jovian atmosphere at 10.5 to 13vm has appre- 
ciable opacity, as discussed in the next paragraph. 

♦ 

The main constituents of the Jovian atmosphere are hydro- 
gen molecules, helium molecules, methane and ammonia, with mi- 
nor constituents .hydrogen sulfide, water, ethane and acetylene. 
Table 6-1 shows their observed abundance ratio by number. 




TABLE 6-1 


Solar composition atmosphere (fraction by number) 


Species 

(1) 

(2) 

H 2 

0.886 

0.870 

He 

0.112 

0.128 

HgO 

1.05 x 10“ 3 

8.80 x 10~ 4 

CH 4 

6.30 x 10“ 4 

6.17 x 10" 4 

NHq 

1.52 x 10“ 4 

1.49 x 10" 4 

mJ 

h 2 s 

,2.90 x 10" 5 

2.56 x 10-5 


(1) Weidenschilling and Lewis (1973)* 

(2) Podolak and Cameron (197*0 ; Cameron (1973). 
The table is taken from Ingrosell. 
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The abundances of hydrogen, methane, ammonia, and helium seem 
consistent judged from solar atomic abundances. For detailed 
information one can refer to McElroy (1973). 

The opacity due to hydrogen molecules is caused by pre- 
ssure induced dipole absorption. The hydrogen molecule has no 
permanent dipole moment, and consequently, no permanent dipole 
spectrum. Gaseous Hg, however, has a weak pressure - induced 
dipole spectrum which absorbs significantly over the long path 
lengths and low pressure of the Jovian atmosphere. The induced 
dipole moment results from two distinct physical process 
(Kranendonk and Kiss, 1959). The first takes place when the 
permanent quadrap ole moment of one molecule induces a dipole 
moment in another molecule by virtue of the neighbor's pola- 
rizability. This is a long range interaction. The second 
physical process takes place when the overlap forces of the 
two adjacent molecules cause an asymmetrical distortion of 
their electronic charge clouds. The net induced dipole moment 
is modulated by the relative translational and rotational mo- 
tion of the colliding pair and this modulation produces the 
absorption of infrared radiation. The translational spectrum 
is predominant at long wavelengths with its peak at 100pm at 
100°K (Trafton and Munch, 1969). In our wavelength region 
(10.5 - 13um) its contribution to the opacity is negligible. 

The rotational hydrogen collisional spectrum, however, has its 
peak at 17pm and contributes a continuous opacity in our wave- 
length region (Th. Encrenaz, 1972). 

The helium molecule also has no permanent dipole moment. 


Its opacity comes from the collision with the hydrogen mole- 
cules and resembles the H 2 -H 2 collision process. The colli- 
sion Is less important due to the smaller abundances of helium. 

Ammonia is an Important source of opacity at lOym under 
Jovian atmospheric conditions. The lOym band of ammonia arises 
from transitions through the v 2 mode. In the v 2 mode of vibra- 
tion, the nitrogen atom oscillates vertically relative to the 
plane of the hydrogen atoms (figure 6- 8a). The nitrogen atom 
is able to penetrate through the potential barrier to the other 
side of the hydrogen plane. This Inverted position leads to 
inversion splitting of the levels of the ammonia molecules. 

The splitting generates both symmetric and antisymmetric energy 
levels with a given vibrational quantum number.. The lOym band 
of ammonia arises from transitions from the ground vibrational 
state to the first excited state In the v Q mode (figure 6-8b). 
Another transition, from the first excited symmetric vibra- 
tional state to the second excited asymmetric state is also In 
lOym range, but the contribution due to this- "hot band" is small 
for the low temperature in the Jovian atmosphere. 

The 'ammonia is clearly seen in abosrption in the spectrum. 
The centers of the bands are shown by the arrows labeled NH^. 

The ammonia absorption has been observed by different groups. 

Gillett et al (1969) (figure 6-9) observed Jupiter from 
2.8-l4ym with low resolution A/AX c 50. Briefly, their results 
show the following: The spectrum has a depression at 3.3ym 

caused by CH^. Solar heating of the upper atmosphere via this 
band results in warming of the upper atmospheric layers. This 
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Figure 6-8 . (a) Schematic representation of the NH^ 

molecule. The components of angular momentum and 
the motion in the vibrational mode are also shown 
(b) Energy levels of the v 2 vibrational mode of am- 
monia. Superscripts a and s refer to the antisymme- 
tric and symmetric levels which arise due to inver- 
sion splitting. ■ 




proceeds until energy is radiated at the same rate via the 7*7wti 
hand of CH^. Ammonia absorption around 10pm was also detected. 
Judging from the CH^ emission, these authors were the first to 


suggest a temperature inversion on Jupiter caused by solar 
heating of the 3.3m band of CHn. They also showed that the 
NHg band at 10pm is saturated, and calculated that the Hg abun- 
dance at 12.5pm, assuming a temperature of 125°K, is 12 km-atm. 
with a pressure P„ <v<l/4atm./ 

• - 11 2. v " 

Ait ken and Jones (1972) obtained a Jovian spectrum from 

8 - 13pm at a resolution X/AX^l43 (figure 6-10). The ammonia 
absorption band is again seen. They estimated that the ammo- 


nia abundance in the band is about 2.7 cra-atm. and a lapse rate 


r- at 13pm given by 
^20km. 


H=-30K, where H is the scale height 


The most recent published infrared spectrum is by Lacy 
et al (1975) who used the Lick Observatory 12 0 M telescope. 


■ 

High resolution spectra were obtained at .8 90 cm (11.24pm) 
with x/ax=1780. Medium resolution data were observed from 
1000 cm -1 to 350 cm" 1 (10 * 12.75pm) with X/AX from 250 to 333 
(figure 6-11 a,bY. The authors also calculated synthetic 
spectra, -assuming that NH^ and H 2 are the only sources of opa- 
city.: Their conclusions from comparison between observed and 

computed spectra follow: All of the prominent lines in their 

observed . spectrum are saturation bands- broadened to a width 
many-" times the pressure — broadened line width. The observed 
135°K continuum is primarily formed by the wings of the NH^ 
line . The H 2 opacity may be important if is unsaturated 


as 



Figure 6-10 . (a) Room temperature absorption spectrum of 

ammonia, p=0 » 06 atmos , w=0 .6 cm atmos . 

(b) Brightness temperature; (c) Surface 
brightness of the central region of Jupiter 
from 8 to 13. 5vm. (Taken from Aitken and 
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near 135°K. A pressure of 0.125 atm. at 135°K is required to 
form the continuum. The minimum temperature in their synthe- 
tic model is 118+5°K while the observed minimum temperature is 
123°K a about* 5 °K larger than the derived temperature due to 
incomplete resolution of the features. The lapse rate at 135°K 
is 7.5±2.5°K/SH. Gillett et al (1969) estimated the lapse rate 
at the NH 0 saturation level is 4°K/SH. A discrepancy occurs ■ 

J ; v 

in the eomparision of the medium resolution data between 870 
and 890 cm*" 1 . In this region the Jovian spectrum seems to be 
depressed by about 2°K relative to the calculated curve. The 
authors suggested that it may be due to an as yet unidentified 
minor constituent of the Jovian atmosphere. 

Our spectrum has about the same resolution as the medium 
spectra of Lacy et, al' s and so the two spectra can be compared. 
The line positions match well. The vertical matching shows a 
drift towards longer wavelength. Figure (6-lla) is obtained 
by matching points A and B. The short wavelength side matches 
but the end of the long wavelength side is about 1.5 times 
higher. Also our spectrum seems to match the' theoretical curve 
better at 870 to 890 cm - ' 5 ’. Figure (6-llb) is obtained by mat- 
ching A 1 and B* . Then the long wavelength side matches better 
than before but the short wavelength side is different. Also 
now our spectrum matches better with Lacy et al’s observed 

— I 

result at 870 - 890 cm A conclusion about the discrepancy 

at 870 cm"” 1 between Lacy's observed and calculated spectra can 
not be reached at present until there is a better way for mat- 
ching our and Lacy's et, al' s spectrum. Also it is not clear 
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Figure 6-lla. Spectrum of the N and S polar regions of Jupiter at 3-1} cm 

— resolution divided by the spectrum of the Moon* Data points 

are shown as solid circles, and the solid line represents 
the best fitting synthetic spectrum calculated from the 
model calculated by Lacy et al (1975) (The graph is taken 
from Lacy et. al ) . The dotted curve is our observed spec- 
trum by matching Lacy's spectrum at points A arid B. 
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whether the difference in matching of our spectrum and Lacy 
et al* s between long wavelength and short wavelength is real 
or not. There could be several reasons for the difference. 

It could be due to the inaccuracy of the end of the spectra, 
because the short wavelength side is the end of our spectrum 
and the long wavelength side is the end of Lacy's spectrum; 
or it may just be due to the matching technique. More effort 
is needed to clearify this point. 

The -bsorption due- to is much less important beyond 
12y, so one may be able to use the Hg opacity to estimate the 
lapse rate in that region. The lapse rate can be estimated 
by the equation 


W - T B (X 2> 



a(X 2 ) 


where T^CX^) 'the brightness temperature at X^ 

T b (X 2 ) is the brightness temperature at X 2 

r is the lapse rate 

a(X^) is the absorption coefficient at X^ 

a(X 2 ) is the absorption coefficient at X 2 

H is the scale height 

If one chooses. X^=11.95, X 2 =12.3^ with measured bright- 
ness temperature 1.^=133.94, T 2 =133. 64, a(X^), a(X 2 ) are taken 
from Calpa and Ketebaar (1957), one obtains a result rH=-43.2°K 
for an H 2 opacity dominated atmosphere. Gillett et al (1969) 
calculated the adiabatic lapse rate (rH)ad--42°K for an Hg at- 
mosphere, while Aitken and Jones (1972) measured a value of 
FH=-30°K from their spectrum. Our value seems closer to the 
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value calculated "by Gillett rather than to Aitken’s. It Is 
emphasized here that the estimate is based on the assumption 
that the H- opacity is the dominant opacity at wavelengths 
11.95v and 12.34y, which may not he true. 

Methane has a strong emission hand at 7.7u hut does not 
have any band structure in our region. Methane has two im- 
portant contribution to the overall thermal structure of the 
Jovian atmosphere. The first one is that methane has an ab- 
sorption band at 3.3tm which absorbs solar energy and rera-; 
diates at 7.7tm producing an inversion temperature layer with 
a maximum temperature of 150°K at an altitude 160 - 200 km. 
Secondly, methane is photo dissociated by ultraviolet light 
in the upper atmosphere. This results in products such as 
ethane (CgHg), acetylene (CgHg), and ethylene (CgH^). Strobel 
(1973) estimated that column densities of CgHg, ^2^2* ^2^i 

0*1 *i £ 

above the cloud top are approximately 10 , 3 x 10 A , 

IK _2 

3 x 10 "cm respectively. CgHg and C 2^2 were observed . 

by Ridgway (1973) using the 60" Kitt Peak solar telescope in 
the 750 - 875 cm” 1 (11.42 - 13. 33um)range with resolution 
A/AX= 770, (figure 6-12a). The lines are shown in strong emi- 
ssion at the 140°K temperature. The apparent lines are super- 
positions of many lines, each group corresponding to a sub- 
band. Ridgway calculated that the mixing ratios are N(CgHg)/ 
N(H 2 )=4 x 10" 3 and N(C 2 H 2 )/N(H 2 )=8 x 10" 5 . The ratio N(C 2 Hg)/ 
N(C 2 H 2 )«|0 where Strobel predicts about 200. Combes et al 
(1974) ’ s (figure 6-12b) observation confirms the presence of 
very strong emission lines of C 2 H 2 and CgHg. The abundance of 
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Figure 6-12 . (a) Thermal emission spectrum of Jupiter cor- 

rected for absorption in the earth’ s . atmosphere observed by 
Ridgway (1973). The dashed line is the predicted form of 
the H 2 continuum, (b) The ratio of the Jovian spectrum to 

the atmospheric absorption spectrum observed by Combes et 
at (197*0 . The solid and dashed lines are the blackbody 
curves at 135°K and 120°K respectively. 









ethane estimated from Ridgway’s spectra depends strongly on the 
distribution of gas temperature, which is not well-determinded. 

If the temperature in the mesospheric inversion layer turns out 
to have a maximum value of 150°K, the observations Indicate 
(Ridgway 197*0 2 x 10” 2 gm cm' -1 of ethane in this high tem- 

perature region. Sagan and Salpeter ( 1976 ) estimate the column 

_o _2 

density of ethane molecules to be 3 x 10 gm cm by assuming 
that ethane is produced by photolysis' of methane by solar ul- 
traviolet photon and destroyed mainly by eddy diffusion Into 
the troposphere, followed by pyrolysis in deeper, hotter layers. 
This would be in very serious conflict with the observations, es- 
pecially since only a small fraction of the theoretical column 
density refer to the hotter inversion region. 

The ethane emission band Is also shown in our spectrum 
(figure 6-6). In the figure the indicated emission line posi- 
tion was extrapolated from Ridgway’s spectrum, and the laboratory 
observed position by, Smith (19*19) are also shown for comparison. 
Our positions agree with Ridgway’s reasonable well, while Smith’s 
seem displaced from ours by O.Olym, possibly due to a uncertain- 
ty in position calibration. Only a portion of the CgHg spectrum 
can be seen in our spectral coverage. The abundance of C^Hg Is 
not estimated here because the absolute amplitude of our spec- 
trum is not well calibrated. 

Our spectrum contains both ammonia and ethane features 
while other observers have not shown both. This will be use- 
ful because one can compute the synthetic spectra including 
both. Ammonia will provide us with information about the top 
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of the cloud layer while ethane provides us with the informa- 
tion about the inversion layer. A synthetic spectrum including 
both ethane and acetylene would be interesting because the in- 
clusion of these new gases would affect the models, especially 
around the inversion. Acetylene would appeal? around 13ym. 

The absorption of solar radiation by CH^ at 3-3ym used to be 
thought to be radiated solely by CHj^ at 7. 7 yin. The 7.7ym emi- 
ssion intensity is a critical test of a temperature inversion 
model and the emission intensity calculated by Wallace et al 
(1974) is within 255S of the value observed by Gillett et al 
(1969) • If ethane and acetylene do contribute to emission in 
the thermal infrared, there must be some additional source of 
solar absorption in order to produce the observed inversion tem- 
perature. Additional absorption at this altitude, perhaps due 
to particles, is suggested by the -low ultraviolet albedo of 
Jupiter in the wavelength region 2100 to 3600 A (Wallace, 

Caldwell and Savage, 1972). 

Terrile and Westphal (1976) had imaged Jupiter at high 
spatial resolution at 8 - 14ym. All images reveal a belt and 
zone structure similar to visible photographs. In the 8 - l4ym 
broad-band data, belts appear to be about 2°K hotter than the 
zone. The lowest belt- zone contrast is found in the hydrogen 
opacity dominated region at 12.5 ym s while images at 9.5 ym have 
the greatest contrast. This is consistent with the dynamic 
picture that zones are rising columns of air and belts are 
sinking columns of air. Ammonia gas being carried upward in 
Zones will freeze’ out and form a thick cloud on top of the zones. 
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giving a low infrared temperature to the zones, and the crys- 
talized NH^ particle will be carried down to the deep atmos- 
phere in the belt where they will get sublimated. One can look 
deeper into clouds in the belt because of the lack of the am- 
monia cloud On top of it and therefore see a higher infrared 
temperature'. Although there are a number of high resolution 
spectra of good spectra for methane. The observation of methane 
will be interesting not only because it provides us with in- 
formation about the inversion layer, it is also useful to find 
out the temperature profile. If the temperature profile of 
the Jovian atmosphere is known, one can use it to find the am- 
monia abundance profile.. Ammonia itself is not a very good 
tool for probing the temperature profile because it has a low 
vapor pressure and its variation is very sensitive to tempera- 
ture changes. Observations of methane at 3.3ym and 7.7ym, 
should be able to accomplish this. 

The Hadamard transform spectrometer described in this 
thesis would be able to make these observations, with small 
modifications that would permit observations to be made at these 
wavelength. In addition, observations should be undertaken at 
8.0 to 9.5 microns where neither ammonia nor methane have strong 
absorption .features. At these wavelengths one would be ob- 
serving the clouds. In this 8.0 to 9.5 micron region our in- 
strument should be able to image the bands and zones of Jupiter, 
to probe for spectral differences and cloud features. Such ob- 
servations should increase our understanding of Jupiter T s cloud 


structure 


(D) Observations of Mercury 


The observations of Mercury were made with the newly built 
Cornell 25" telescope at Mount Pleasant, Ithaca, New York. 

The telescope has a focal ratio f/13.5 * The spectrometer’s 
acceptance . beam size is 7Tb” x 78” . The dewar described in 
section III-B was used with the spectrometer. The spectrometer 
operated in the lD.5^13um region with a resolution x/ax around 
300. The chopping frequency was 10 cycles per second. 

The observational procedure was carried out a little dif- 
ferently from the observations of a-orionis and Jupiter. 

First, Sun spectra were used for correction spectra rather than 
lunar spectra. There are no known' molecular lines in this re- 
gion. Its temperature at ll.lOym is 5030°K(Saildy & Goody). 

The 'observations were made on August 3, 1976. At that time the 
Sun was about an hour away from Mercury, and was observed 
through roughly the same air mass as Mercury. Secondly, Mer- 
cury is so faint in broad day light that we were unable to see 
it in visible light. The- way to find Mercury was the following: 
We pointed the telescope in the correct region and scanned for 
the infrared signal. The signal is so strong that one can see 
it go off scale on the synchronous demodulator. The pointing 
accuracy of the telescope is 6 sec. of time In right ascension 
and 25 sec. of arc in declination. Thirdly, since we could not 
see Mercury visually for tracking, we adapted a* different me- 
thod for tracking Mercury. Since our computer programming is 

r 

set up in such a way that the data taken when the mask is 
moving forward and moving backward are stored in different 


areas, we only took data when the mask was moving forward. 

When the mask was moving backward we maximized the signal to 
assure correct pointing and waited for the next forward pass. 
An> noise introduced when moving the telescope for maximizing 
the signal would have gone into the backward-pass data bins ■ 
and those data points were thrown away anyway. Since each pass 
takes 51 sec. only. Mercury remained at essentially the same 
position during the forward data taking pass. 

Two runs of Mercury and two runs of the Sun were taken. 

4 

Each run consisted of ten scans of the sources. Figure (6-13) 
shows the raw spectra of Mercury and the Sun. Since Ithaca 
has a lot of moisture in the air during the summer time, the 
correction for atmospheric features is more difficult than at 
Kitt Peak and is done in a different way. A constant was added 
to the Mercury spectra such that atmospheric features in the 
Mercury/Sun ratio spectrum were minimized. This step ensures 
that the atmospheric features are largely corrected. The Mer- 
cury spectrum which has a constant added to it was then multi- 
plied by another constant to make its amplitude as close to 
that of the solar spectrum as possible. The solar spectrum 
was then subtracted from the modified Mercury spectrum. The 
multiplication of the Mercury spectrum by a constant assured 
that atmospheric features in the two spectra had similar am- 
plitudes before the subtraction step. This difference spec- 
trum was then added to a '‘perfect 11 solar spectrum which is 
calculated according to the > blackbody function appropriate to 
the solar temperature. What one gets from these procedures is: 
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Fleure 6-13. The raw spectra of Mercury and the Sun. 





Pinal Mercury Spectrum 

«•. observed corrected Mercury spectrum - observed solar 
spectrum + perfect solar spectrum 
= ("perfect 11 Mercury spectrum 4- noise in tbe Mercury spec- 
trum) - ("perfect” solar spectrum + noise in the solar 
spectrum) + "perfect” solar spectrum 
- "perfect” Mercury spectrum + (noise in the Mercury spec- 


trum - noise in the solar spectrum) 

Any systematic noise such as emission and absorption due to the 
sky or to the telescope will be subtracted away. The advantage 


of this method is .that the final spectrum is obtained through 


subtraction rather than by division. Division, in the low sig- 
nal portion of the spectrum, produces deceptive high noise • 
•spikes in the ratio spectrum. The method we have used tends 
to eliminate these. 

Since this will be the first high resolution Mercury spec- 
trum obtained, we will calculate Mercury’s disk integrated in- 
frared temperature and compare this temperature with the observed 
one. Morrison and Sagan (1967) had calculated the infrared 
brightness temperature of the center.-of-disk as a function of 
phase angle and heliocentric longitude, but there is no cisk 
integrated infrared temperature available in the literature. 

I 1 In the following we will discuss the factors that may affect 

the brightness temperature, and then present a method of cal- 
culating the disk integrated infrared brightness and compare it 
with our observation. A good review of thermophysics of ’Mer- 
j cury is given by Morrison (1970). 
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In 1965 a Fettingill and Dyce (1965) used radar to discover 
that Mercury has a rotation period of 59 days, two thirds of 
the orbital period* instead of an 88 day synchronous rotation 
around the Sun. This implies that Mercury has a solar day* 
on the planet, about 176 terrestrial days long, equal to two 
orbital revolutions in three rotations. Mercury’s non-synchro- 
nous rotation leads to time-dependent thermal emission of the 
planet due to the diurnal variation of the insolation. This 
diurnal variation would not happen for a synchronously rotating 
Mercury. The diurnal variation changes the brightness tempera- 
ture as a function of both phase angle and heliocentric longi- 
tude. It also allows a measurement of the thermal properties 
of Mercury’s surface. . 

Because of the high eccentricity of Mercury’s orbit, 

(e- 0.2), the diurnal cycle of insolation is markedly different 
from longitude to longitude, and can differ by a factor of 2.5. 
The eccentricity enters in two ways. First, the variation in 
distance from the Sun produces a solar constant that varies by 
more than a factor of 2 from perihelion to aphelion. Second, 
the changing orbital angular velocity causes the apparent speed 
of the Sun across the sky to vary; near perihelion the angular 
velocity of revolution actually slightly exceeds the angular 
velocity of rotation, and the apparent planetocentric solar mo- 
tion is retrograde (figure 6-14, Soter and Ulrichs, 1967). 

The two effects of the eccentricity reinforce one another, with 
the larger flux coming at a time when the angular rate of the 
Sun across the sky is largest. The two longitudes (180° apart) 



The relative positions of the Sun are marked at 11 day 
intervals with the planet held as a fixed reference. 
Planeto-graphic longitude are indicated for Mercury. 
(Taken from Sot er and Ulrichs, 1967). 
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that see the Sun overhead at perihelion receive more than two 
and a half as much energy per period as the longitude 90° away, 
where the Sun is always small and rapidly moving while near the 
zenith. 

Besides the insolation geometry, a possible atmosphere on 
Mercury will also affect the thermal emission. C0 2 is a major 
product for a possible secondary atmosphere, furthermore, since 
C0 2 could have been photodissociated and reduced to CO by pre- 
ferential loss of oxygen. Pink et_ al (1973) had set up a search 
for a possible Mercury atmosphere of C0 2 and CO. They set up 

_ii 

an upper limit 1.0 x 10 mb surface pressure for C0 2 and 2.0 x 
10~^mb for CO. Mariner 10’ s results also suggest that Mercury 
has no atmosphere although it may have a thin layer of He and 
other inert gas trapped by Mercury’s magnetic field. 

The optical observations of Mercury show that the integral 
spectral reflectivity of Mercury is quite similar to that for 
the integral moon (McCord and Adams, 1972). The Bond albedo 
of Mercury (de Vancouleurs , 1964) is 0.058. 

In our model calculation we assume the following things. 
This model has been described by Murdock (1974): 

1. The emission from the dark .side at the phase angle we 
observed is negligible. On the day we made our observations, 
August 3, 1976, the illuminated portion- amounted to 0.973 of 
the total disk and the dark portion was 0,017. The dark side 
temperature is about 110°K, which at 10nm has a flux 6.7181 x 
10" 8 watt-cm" 2 -u" 1 -sr" 1 . The flux for 500°K at 10ym is 7.1007 x 
10 watt-cm -v -sr , so the contribution from the dark side 


is negligible. 

2. At infrared wavelengths, the radiation that reaches the 
observer originates very near the surface, so the infrared tem- 
perature is assumed equal to the insolation temperature. Soter 
and Ultichs’ (1967) results show that the day time temperature 
is independent of the thermal properties of the surface material 
and determined largely by the insolation temperature. 

3. The infrared emissivity we assumed was 0.9, which is the 
lunar value. We choose this value since Mercury’s surface may 
be similar to the lunar surface. 

In figure (6-15) we choose two coordinate systems on Mer- 
cury surface. The unprimed system Is the ’’solar system” with 
the z-axis pointing towards the Sun. A Is the subsolar point. 
The hemisphere above the plane BCD facing the Sun is the illu- 
minated part. The primed system is the ’’earth system’ with the 
z’-axis pointing towards the earth. A’ is the subearth point. 
The hemisphere above the plane B’C’D’ facing the earth is the 
portion that is being seen from earth. The two systems are 
different by an angle a with the x-axis as the common axis. 

If the surface is in equilibrium with sunlight and cannot 
conduct heat away, the subsolar point temperature Is: 


S(l-A) 1/4 
ceR 2 


(6-1) 


6 -2 

where S Is solar constant at earth, equal to 1.360x10 erg cm s 
A is the Bond albedo for Mercury, assumed to be 0.058. 

o Is the stefan-Boltzmann constant equal to 5.67x10”^ 

- 2-1 
erg cm s 
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Figure 6-15 • Two coordinate systems on the surface of 

Mercury, The unprimed system is the "solar 
system" with the Z-axis pointing towards the 
Sun, A is the subsolar point. The primed 
system is the "earth system" with the Z’-axis 
pointing towards the earth. A’ is the sub- 
earth point. 
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e is the infrared emissivity at lOym from the surface, 
assumed to be 0.9 

R is the distance between Mercury and the Sun in astro- 
nomical units 

Mercury's subsolar point temperature varies with distance from 
— 1/2 

the Sun as R and therefore is a function of heliocentric 

longitude due to the eccentricity we discuss above. 

In the unprimed system the temperature distribution on 
the surface will be concentric isothermal bands around the sub- 
solar point. The temperature of a band at colatitude e is given 
by 

T(8) * T Q cos 1/4 0 (6-2) 

The total intensity at any wavelength region will be composed 
of the contributions at that wavelength from many isothermal 
regions each with its own apparent area. 

Since one is interested in the flux coming to the earth 
one will have, 

dF£ = I x (e»,<j>» ) dA* 

= r 2 sine 1 d0 , d4> ! 

= 5 X (0’ ) cos 6 * r 2 sin0 ’ d8 1 d<j> 5 


(6-3) 

(6-4) 

(6-5) 


where dF* is the flux at wavelength. X coming from the area dA* 

A 

I x is the Intensity at ( e ' a 4> T ) 

dA' is the differential area on Mercury 

r is the radius of Mercury 
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E^(e ’ ,$’ ) 1^(8 1 a 4> t )/cosfl 1 Is the component of flux 

that radiates toward the 

i 

earth. 

3^(6’ ) Is a complicated function of ( 0 T a ^ T ) bacause the 

temperature distribution is a complicated function of (8 , ,<f’ , )» 
However, one can convert the system to the "sun coordinate" 
where s^(e,$) is a simple function. 

Since * dA* = dA (6-6) 

cosQ 1 - cosecosa - sines in <j>sina (6-7) 

equation (6-5) becomes 


where 


where 

and 
so • 

The integral (6-9) can be separated into two parts, for 0<e_ 
where 8. » t? - et (6-10) 

is the limit of the cap shared by both the Sun and the earth, 
l.e. the limit of integration over de can 4> go from 0 to 2w. 
For e>0., the limit of the integral over d is constrained by 


dF^ - 3^(e) C cos0cosot “ sines in <j>s inti] r 2 sin0d8d«j> 

V 6) = ^ e 0 2 /kT(6) ~ (6_8) 

Cjs* 1.1909 x lO^watt cra” 1 1 i" 1 sr“ 1 

C 2 = 1.4388 x lo\iK 

T(e) is given by (6-2) 

F^ = r 2 /de/d(j> H^(e)(sinecos6cosa - sin 2 esin<j>sina) 

(6-9) 
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the physical .condition that some part on <j> that is illuminated 
by the Sun can not be seen from the earth. 

So (6-9) becomes 


3? ^ 

* X 


jB 2v 

r 2 / de/ dtf) S , ( 6 ) (sin ecos ecos a - sin 2 a sin$sina) 
o ' o . A 


TT <j> ( 0 ) 

+ r 2 / T de/ 2 dtji =,( e) (sinecosecosa - sin 2 esin<i>siha) 


<h( e) 


( 6 - 11 ) 


To find $i(e) and $ 2 (e) one notices that $ is given when e Ts =90°. 
From (6-7) 

cose* = cosecosa - sinesinijisina 


e* 

and 

$ ■ 


j =$ sin<i> * cotacote 

— 1 

sin (cotacote) 


(6-12) 


(6-ljs) 


Since $ will also be symmetric about, the y-axis 9 equation (6-11) 
can be rewritten as: 

TT ' 

rj T “Ct PIff * * 
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A 0 0 A 
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2 “ a sin 2 esin<j.sina) 

(6-14) 


After evaluating the integral one obtains the following: 
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2irr 2 cosa 

* ( 6 ) sin 2 e 
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is the mean of ( E^(e) ) in the interval of 9 from 0 to 

-a <, The same meaning applies to the third term of (6-15) • 
As a check of equation (6-15), if a=0, that means when 
subsolar point and subearth point coinside s ,let S( 0)= constant 
evaluating (6-15) gives: 

F * trr 2 S 

If a- j j that means the s.ubsolar point and subearth point are 
90° apart. With 5(0) assumed to be constant, equation (6-15) 
gives 



which is as one expects since one is seeing half of Mercury. 

Equation (6-15) can be readily integrated on' a computer. 

It is applied to our case with the following physical parameters: 
Date: August 3 S 1976 . 

Phase Angle: 53° 





I 13*J 

» i . ■ 

' t 

Radius vector: 0.414 A.U; 

Orbital longitude: 19 8. 09° 

Mercury perihelion point: 0.3075 A.U* 

Subsolar temperature at perihelion point: 700°K 

Subsolar temperature at ct=53 : 603 K 
Equation (6-15) was computed on a LSI mini computer at each 
wavelength., from 10.6 to 13.2ym. The program and the result are 
shown in the Appendix D. The integral was divided into twenty 
steps. H ( 6 ) was evaluated from equation (6-8) 


A -a 5 e c 2 /kT(e )^ 1 

where T(e) = T cos 1 ^^ 


The calculated spectrum is a measure of color temperature 
and is used to compare with the observed spectrum. Figure (6-16) 
shows the final Mercury spectrum conected for solar temperature, 
with a number of blackbody slopes shown to match. The calculated 
spectrum (cross) matches the blackbody temperature 525 °K, which 
also matches the observed ‘spectrum. We concluded that the best 
fit lies in the 500°K region. Murdock (197*0 measured a effec- 
tive brightness temperature at 10.8um at the same phase angle 
to be around 650°K. Our results disagree with his results. 


o 


Relative Flux 



Figure 6-1 6 . The final Mercury spectrum, corrected for solar tem- 
perature, with a number of blackbody slopes shown to 
match. 


APPENDIX A 


ESTIMATE OP CODING ERROR FOR FOURIER 
* TRANSFORM SPECTROMETRY 

Let 5 be the path difference in a two beam interferometer. 
?(v) is the power at wavelength v s (i.e. the spectral density 
function) v here Is taken to.be 1/x. S(?) is the power received 

for path difference Then (p .96 Stewart) 

* * , 

S(£) «* /“P(v) cos 2 (2tt5v) dv A-l 

0 

« 1/2P + l/2/*P(v) cos (-4*^) dv A-2 

0 0 

The reciprocal Fourier property is (Morse and Feshback P.45*0 
that If 

FU) = M- /”cosUv)f (v)dv 
* 0 

Then f(v) = /*°eos(5v)F(5)dS 

0 

which implies, neglecting the constant term, that 

P(v) « l6/"cos(%5v)s(5)de A- 3 

0 

Suppose we take measurement at (N+l) equally separated 
step* s in the variable £. Let the step length be x» then 

5 ~ ht 

In general x is chosen such that 
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i ; A-4 
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which, If v . <<v . effectively implies sampling twice per 

mxn ihex 

cycle (Stewart, p.303) . 

One can now write the integral for P(v) as 
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Now consider frequency v= + m$ 


where 
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and m is an integer. m=0,l* • ■ • N 

The 6P(v) is the power in the resolved spectral band width 
6 about frequency v 
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Writing this out in matrix form* with 9 1 =4irrv in ^ 8 £ =4Trc(v mi r> +S) 


and with ir ( ) =6 P ( v jnin ) we have 
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The SNR improvement is the reciprocal of this quantity. 
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0066 

3408 


JSS 

ST0P 

SS D0WN? 

01 19 

0067 

43 07 


SSN 

; 07 

N0j. STOP DOWN? 

0120 

0068 

F206 


JMP 

ST0P 

YES 

0121 

0069 

FA56 


J ST 

AOC 

N0j. SAVE DATA 

0122 

006A 

E24D 


LDX 

N0SM 


0123 

0063 

■ OHO 


ZAR 



0124 

006C 

FIOO 


JMP 

DISP 


0125 

006D 

B241 

RTTR 

LDA 

DIR 


0126 

0063 

F642 


JMP 

NEXT 


0127 

006F 

£243 

STOP 

LDX 

MI 

. • 

0128 

0070 

C60f 


LAP 

r 


0129 

0071 

F100 


JMP 

DI SP 


0130 

0072 

C63F 

RT0P 

LAP 

‘? * 


0131 

0073 

F9 00 


J ST 

BTC 


0132 

0074 

01 IF 


DATA MSG 


0133 

0075 

F900 


J ST 

IKB 


0134 

0076 

COCE 


CAI 

*N * 

A30RT LAST PASS? 

0135 

0077 ’ 

F203 


JMP 

N0“ 

N0 

0136 

0078 

C0D9 


CAI 

»Y» 

YES 

0137 

0079 

F203 


JMP 

YES 

YES 

0133 

007A 

F608 


JMP 

RT0P 

WRONG ENTRY 

0139 

0073 

FA44 

N0 

J ST 

ACC 

SAVE DATA 

0140 

007 C 

F203 


JMP 

Y2 


0,141 

007D 

C6AC 

YES 

LAP 

'* * 

PRINT ABORT 

0,142 

007E 

F9 00 


J ST 

BTC 


0143 

007F 

0128 


DATA 

ABT 


0144 

OOSO 

3234 

Y2 

LDA 

CT 

SET COUNTER 

0145 

003 1 

9A30 


STA 

SCT 


0146 

0082 

0350 


ARP 


RESET SUBSCRIPTS 

0147 

003 3 

99 00 


STA 

XB 


0148 

0084 

9900 


STA 

xc 


0149 

003 5 

E230 


LEX 

NPBrf 


0150 



*ADD 

T0 BEAM ARRAY 


0151 

0036 

EE33 

.ALP 

STX 

IBP 


0152 

0087 

B400 

' 

LDA 

§0 

GET DATA 

0153 

0038 

S401 


LDX 

ex 


0154 

0039 

F9 00 


J ST 

DPFLT 

FL0AT IT 

0155 

008 A 

9A21 


STA 

MSIC1 

SAVE 1 T 

0156 

0033 

EA21 


STX 

M5K2 


01*57 

003 C 

F9 00 


J ST 

FAD 

ADD T0 BEAM 

0158 

008 D 

OOAC 


DATA 

MSK1 

* 

0159 

003 S 

0000 

BEAM 

RES 

2> O' 


0160 

009 0 

D900 


IMS 

XB 

BUMP SUBSCRIPTS 

0161 

0091 

D9 00 


IMS 

XC 


0162 

009 2 

E63F 


LDX 

IBP 



0163 

009 3 

C202 


AX I 

2 

0164 

0094 

DA ID 


IMS 

SCT 

0165 

009 5 

F60F 


JMP 

ALP 

0166 



* PRINT 

PASS 

COUNT 

0167 

009 6 

F9 00 


J ST 

CRLF 

0168 

009 7 

C702 


LAM • 

2 

0169 

0098 

9A23 

• 

STA 

C0N1 

0170 

0099 

Z21G 


LDX 

NPEM 

0171 

009 A 

0110 

MK 

ZAR 


0172 

009 B 

E401 


LDX 

@1 

0173 

009 C 

F9 00 


J ST 

DPFLT 

0174 

009 D 

9A0E 


STA 

MSK1 

0175 

009 E 

EAOE 


SIX 

MSK2 

0176 

009 F 

F9 00 


J ST 

0FPA 

0177- 

00A0 

OOAC 


DATA 

MSK 1 

0178 

00A1 

0110 . 


ZAR 


0179 

00A2 

F900 


J ST 

0TL 

018 0 

00A3 

0130 


DATA 

RUNS 

018 1 

00A4 

E212 


LDX 

NMEM 

0182 

00A5 

DA16' 


IMS 

CON 1 

018 3 

00A6 

F60C 


JMP 

MK ' 

0184 

00A7 

F7A7 

T2 

RTN ■ 

TRANS 

0185 



* DATA 

STORAGE 

018 6 

00A8 ' 

FAOO 

Z CT 

DATA 

-1536 

0187 

00A9 

0000 

CNT 

DATA O; 0 

* 

OOAA 

0000 




0188 

OOAB 

1400. 

NT 

DATA 

: 1400 

0189 

OOAC 

0000 

M 5K 1 

DATA 

O' 

019 0 

OOAD 

0000 

MSK2 

DATA 

0 

019 1 

OOAE 

OOFF 

H255 

DATA 

255 

0192 

OOAF 

0000 

DIR 

DATA 

0 

0193 

OOBO 

0000 

DOT 

DATA 

0 

0194 

00B1 

0000 

DCT2 

DATA 

0 

019 5 

0032 

0000 

SCT 

DATA 

0 

0196 

0033 

18 00 

N I 

DATA 

: 18 00 

0197 

0034 

18 02 

IB 

DATA 

: 18 02 

0198 

0035 

FEOO 

CT 

DATA 

-512 

0199 

0056 

1400 

NPEM 

DATA 

Y 1400 

0200 

0037 

1600 

NMB-i 

DATA 

* 1600 

0201 

0038 

0000 

N0EM 

DATA 

0‘ 

0202 

0039 

0000 

THEM ‘ 

DATA 

0 

0203 

003A 

9002 

PL S 

DATA 

: 9002 

0204 

D03B 

9202 

MINS 

DATA 

: 9 202 

0205 

003C 

0000 

CON 1 

DATA 

0 

0206 

003D 

opoo 

TETR 

DATA 

0 

0207 

OOBE 

1000 

BPNT 

data 

: 1000 

0208 

OOSF 


MASK 

REF 


0209 



* 



02*10 



* ADD INPUT , 

ARRAY TO 

•021 1 



* 



0212 

OOCO 

08 00 

ACC 

ENT 

, 

0213 

00C1 

B604 


LDA 

TETR 

0214 

00C2 

208 3 


JAM 

£+4 

0215 

0003 

360A 


LDA 

THEM 


DONE? 
NO ' 


TEMP ARRAY 


1 



0216 
0217 
02 IB 

0219 

0220 
0221 
0222 

0223 

0224 

0225 
022 6 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 
024 0 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 
0261 
0262 
02.63 

0264 

0265 

0266 

0267 

0268 
0269 
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00C4 

0210 


CAR 



00C5 

F20 1 


JMP 

S+2 


00C6 

B617 


L DA 

DIR 


00C7 

COOO 


CAI 

0 


OOC8 

F202 


JMP 

S+3 


OOC9 

E613 


LDX 

NPEM 


OOCA 

F201 


JMP 

5+2 


OOCB 

EG 1 4 


LEK 

NMEM 


OOCC 

EEl4 


STX 

N0EM 


00 CD 

DC01 


IMS 

@1 

BUMP COUNTER 

OOCE 

C202 


AXI 

2* 


OOCF 

EA1 2 


SIX 

TBP 

SAVE TEMP POINTER 

OODO 

C7FF 


LAM 

255 

DATA C0UMT ' 

OOD1 

9E28 


STA 

CNT 


00D2 

E61F 


LEK 

NI 

. INPUT POINTER 

00D3 

C202 

Al 

AXI 

2 

• BUMP C0UNTER 

00D4 

EE8 1 


SIX 

IBP 

+ SAVS 

00D5 

B618 


LDA 

TETR 

• 

00D6 

2188 


JAL 

T5 

, 

00D7 

3400 ' 


LDA 

@0 


OODS 

E401 


LDK 

@1 


00D9 

1328 


LLX 

1 ' 


OODA 

IBS 6 


LLR 

7 


00D3 

1 3A3 


LEX 

1 


OODC 

1356 


LLA 

7 


OODD 

10D6 


ARA 

7 


OODE 

F202 


JMP 

T6 


OODF 

3400 

T5 

LDA * 

eo ' 


OOEO 

E401 


LDX 

61 


OOE1 

F9 00 

T6 

J ST 

DPACC 

* 

00E2 

0000 

1 at* 

DATA 

0 


OOE3 

C202 


AXI 

2 

BUMP TEMP 

00E4 

EE02 


5TX 

TBP 


00E5 

E69 2 ‘ 


LDX 

IBP 


OOE6 

0110 


EAR 



00E7 

9 COO 


STA 

@0 


00E8 

9 C01 


STA 

91 


OOE9 

DE40 


IMS 

CNT 

D0NE? 

OOEA 

F617 


JMP 

Al 

N0 ' 

OOEB 

F72B 


RTN 

ACC 



* 


00EC 08 00 
00ED 4006 
00EZ 45 Cl 
OOEF F60 1 
OOFO 0E00 
OOF 1 49 Cl 

OOF2 F60 1 
00F3 0F00 

00F4 40C4 

OOFS 4005 
00F6 F70A 


*WAIT.F0R ALIGNMENT PILSE 


TURN 0N BIT 
Cl D 

SSN Cl 
JfrlP s- 1 
SS*1 

SEN : Cl 
JMP 5- 1 

sum 

SEL : C4 
CIS 
RTN 


DI -SABLE AUT0 
LIGHT OFF? 

N0 ' * 

YES BYTE 0N 
LIGHT ON? 

NO 

YES, BYTE 0FF 
CLEAR FLAG ' 
ENABLE AUT0 


TURN0N 


0270 



* 




0271 



* INPUT 

FROM 

MON SANTO 


0272 



* CONVERT BCD T0 BINARY 

027 3 



* 




0274 

00F7 

08 00 

M0N S 

BUT 



0275 

OOFS 

B633 


LDA 

TETR 


0276 

00F9 

208F 


JAM 

T3 


0277 

OOFA 

C6FF 


LAP 

255 


027S 

OOF3 

8E4B 


ADD 

DCT • 


0279 

OOFC 

1050 


ALA 

1 


028 0 

OOFD 

£644 


LEX 

TREM 


026 1 

OOFE 

38 02 


JXN 

S+3 


0282 

OOFF 

8 AOS 


ADD 

EMTR 


028 3 

0100 

F201 


JMP 

£+2 


0284 

0101 

8A05 


ADD 

DPTR 


028 5 

0102 

0048 


TAX 



028 6 

0103 

3400 


LDA 

@0 


028 7 

0104 

£401 


LEX 

§1 


0288 

0105 

F9 00 


J ST 

DPFIX 


0289 

0106 

F70F 

, 

RTN 

MSN S 


029 0 

0107 

1002 

DPTR 

DATA 

: 1002 


029 1 

0108 

1202 

H-1TR 

DATA 

: 1202 


029 2 

0109 

49 C6 

T3 

SEN 

: C6 

FLAG 

029 3 

010A 

F60 1 


JMP 

£- 1 

NO 

029 4 

0103 

5 AC 6 


I NX 

:'C6 

YES* INPUT TO X 

029 S 

01 oc 

2843 


JXZ 

£-3 

IGNORE ZEROES 

029 6 

OlOD 

C704 


LAM 

4 * * 

SET DIGIT COUNT 

029 7 

OIOS 

9 £65 


sta 

-CNT 

* 

0298 

01 OF 

0110 


ZAR 


CLEAR TALLY 

0299 

Olio 

9 £63 


STA 

MSK2 

‘ 

0300 

Oil 1 

F206 


JMP 

M2 


0301 

0112 

1350 

M 1 

LLA 

1 

X10 

0302 

0113 

9E66 


STA 

MSK2 


0303 

0114 

1351 


LLA 

2 


0304 

0115 

8 £68 


ADD 

.MSK2 


0305 

0116 

9 £69 


STA 

MSK2 


0306 

0117 

0110 


ZAR 



0307 

0118 

1303 

M2 

LLL 

4 

GET NEXT DIGIT 

0308 

0119 

8E6C 


ADD 

MSK2 

ADD TO TALLY 

0309 

01 1A 

DE71 


IMS 

CNT 

LAST DIGIT? 

03.10 

0113 

F609 


JMP 

Ml 

NO 

0311 

one 

0048 


TAX 



0312 

01 ID 

0110 


ZAR 


* 

0313 

oils 

F727 


RTN 

MON S 

YES 

0314 

* 


*TEXT 

STORAGE 


0315 

01 IF 

SD3A 

MSG 

DATA 

: SD3A 


0316 

0120 

C1C2 


TEXT 

•ABORT LAST RUN? • 

« 

0121 

CFD2 




f 


0122 

D4A0 


• 




0123 

CCC1 






0124 

D3D4 



* 



0125 

A0D2 




• 


0126 

D5CE 






0127 

BFAO 







0317 0128 8 06 A ABT DATA sSDSAi 

0318 . 0129 B1A0 TEXT • 1 RUN ABORTED/ * 

012A D2D5 
0123 CEAO 

0 12C Cl C2 ... 

01 2D CFD2 
012E D4C5 
012F C4AC 

0319 0130 A0D2 RUNS TEXT * RUNS KEPT* 

0131 D5CE 

. 0132 D3A0 

■ 0133 CBC5 

0134 D0D4 

0320 0135 SD8A DATA : SDSAj 0 

0136 0000 

0321 0137 8D5A M0DE DATA :SDSA 

0322 0138 CDCF TEXT *M0N 5ANT0/ TAPE 0R INVERSE? 

0139 CED3 ... 

013A C1CE 

013B D4CF ' 

0 1 3C ACD4 
01 3D Cl DO 
0 1 3E C5A0 
0 1 3F CFD2 
0140. A0C9 
' 0 14 1 CED6 

0142 C5D2 

0143 D3C5 

0144 3FA0 

0323 END 
0000 ERRORS 




150 


0001 




NAM 

IKB, IER. BIPT,XEQ 


0002 




NAM 

OTT, ERR. CRLF 


0003 




NAM 

0TL..0FPA 



0004 

0000 



REL 

0 



0005 



* 





0006 



* PANIC SWITCH 



0007 



* 





0008 

0000 

FBOO 


J ST 

*CMND 



0009 

0001 


cm d 

REF 




0010 



* 





001 1 



*KEY BOARD INPUT 



0012 



* 





0013 

0002 

03 00 

1KB 

ENT 




0014 

0003 

4033 


SEL 

7, 0 

AUTO- ECHO 


0015 

0004 

4039 


SEL 

7, 1 

KBD M0DE 


0016 

0005 

59 39 


RDA 

7, 1 

READ 0N FLAG 


0017 

0006 

403C 


SEL 

7, 4 

RESET 


0018 

0007 

F705 


RTM 

IKB 



0019 

0003 

05 00 

IER 

ENT 




0020 

0009 

FE07 


J ST 

IKB 



0021 

OOOA 

CODF 


CAI 

! DF 

BACK ARROW? 


0022 

0003 

FFOA 


J ST 

* CMN D 

YES 


0023 

OOOC 

COS A 


CAI 

: 8 A 

LINE FEED? 


0024 

OOOD 

FFOC 


J ST 

Hs CMN D 

YES 


0025 

OOOE 

F706 


RTN 

IER 

NEITHER. RETURN 


0026 



* 





0027 



* PAPER 

TAPE 

; INPUT 



0023 



^ DSO 

= 0 

FGR TTY 



0029 



* 

1 

F3R HSR 



0030 






- 


0031 

OOOF 

03 00 

BI FT 

ENT 




0032 

0010 

5S01 

BIP2 

I SA 


READ SWITCHES 


0033 

001 1 

13D0 


LRA 

1 

DSO UP FOR TTY 


0034 

0012 

220A 


J0S 

HSR 

DOWN F0R HSR 


0035 

0013 

49 3B 


SEN 

7, 3 

TTY BUSY? 


0036 

0014 

F604 


JMP 

BIP2 

yes 


0037 

0015 

403A 


SEL 

7, 2 

NO, STEP READER 


0033 

0016 

48 39 

VT 

SSN 

7* 1 

FLAG? 


0039 

0017 

F203 


JMP 

IT' 

•YES 


0040 

0018 

0150 


IAR 


NO, BUMP C0UNT 


0041 

0019 

2149 


J AZ 

BIP2 

RESTART IF TIME 

UP 

0042 

001A 

F604 


JMP 

WT 

ELSE CHECK FLAG 

AGAIN 

0043 

00 IB 

53 33 

IT 

INA 

7, 0 

INPUT FR0M TTY 

* 

0044 

001C 

F70D 


RTN 

BI PT 

& RETURN 


0045 

001D 

49 33 

HSR 

SEN 

6, 3 

HSR BUSY? 


0046 

001E 

F60E 


JMP 

BI P2 

yes 


0047 

00 IF 

4032 


SEL 

6; 2 

NO, STEP READER 


0043 

0020 

4S 35 

UK 

SSN 

6, 5 

FLAG? 


0049 

0021 

F203 


JMP 

IH 

YES 


0050 

0022 

0150 


IAR 


NO.- BUMP COUNT 


0051 

0023 

2153 


JAZ 

. BI P2 

RESTART IF TIME 

UP 

0052 

0024 

F604 


JMP 

UH 

ELSE CHECK FLAG 

AGAIN 

0053 

0025 

53 35 

IH 

INA 

6, 5 

INPUT FROM HSR 




■ 

\ 







J 

— t- 


151 . 


0054 

0026 

F717 


RTN 

BIPT 

cl RETURN 

0055 



* 


* 


0056 



-IS WAIT 

FOR EXECUTE SIGNAL 

0057 



* 




0058 

0027 

0800 

XEQ 

ENT 



0059 

0028 

FE20 


J ST 

IER 

INPUT 

0060 

0029 

COB D 


CAI 

i 8 D 

CARRIAGE RETURI'J ? 

0061 

002A 

F703 


RTN 

XEQ 

YES, RETURN 

0062 

002B 

F603 


JMP 

XEQ+ 1 

NO, GET MORE 

0063 



* 




0064 



* OUTPUT T0 TTY 


0065 



* 




0066 

002C 

0800 

0TT 

ENT 



0067 

002D 

4030 


SEL 

7*4 

RESET INTERFACE 

0063 

002E 

6D3B 


WRA 

7* 3 

WRITE 0N NOT BUSY 

0069 

002F 

49 3B 


SEN 

7* 3 

D3NE? 

0070 

0030 

F601 


JMP 

S- 1 

NO 

0071 

0031 

F705 


RTN 

0TT 

YES 

0072 



* 




0073 



* COMMAND ERROR EXIT 


0074 



* 




0075 

0032 

03 00 

ERR 

ENT 



007 6 

0033 

C6DF 


LAP 

: DF 

PRINT ARROW 

0077 

0034 

FE08 


oST 

0TT 


0078 

0035 

FF34 


J ST 

*CMND 

RESTART COMMAND 

0079 



* 




003 0 



* CARRI AGE- RETURN, LINE 

FEED 

008 1 



* 




008 2 

0036 

08 00 

CRLF 

r'jT 



008 3 

0037 

068 D 


LAP 

: 8 D 

CR 

0084 

0038 

FEOC 


J ST 

0TT 


008 5 

0039 

068 A 


LAP 

:8A 

LF 

008 6 

003A 

FEOE 


J ST 

0TT 


0087 

003B 

F705 


RTN 

CRLF 


0088 



* 




0089 



* OUTPUT TEXT 

FROM BUFFER 

009 0 







009 1 

0030 

08 00 

0TL 

ENT 



009 2 

003D 

SAOF 


ADD 

CAI 

MAKE COMPARE INSTRUCTION 

009 3 

003E 

9 A06 


STA 

0T2 

&SAVE IT 

0094 

00 3F 

9A09 


STA 

0T3 


009 5 

0040 

E704 


LDX 

*0TL 

GET TEXT POINTER 

009 6 

0041 

PE05 


IMS 

0TL 

SET RETURN ADDRESS 

009 7 

0042 

B400 

0T1 

LDA 

@0 

GET WORD 

0098 

0043 

11D7 


RRA 

8 

PRINT FIRST BYTE 

0099 

0044 

FE1S 


J ST 

0TT 


0100 

0045 

coco 

0T2 

CAI 

0 

LAST ONE • 

0101 

0046 

F70A 


RTN 

0TL 

Y ES, RETURN 

0f02 

0047 

1157 


RL A 

8 

PRINT SECOND BYTE 

0103 

0048 

FE1C 


. J ST 

0TT 


0104 

0049 

0000 

0T3. 

CAI 

0 

LAST ONE? 

0105 

004 A 

F70E 


RTN 

’ 0TL 

YES RETURN 

0106 

004B 

0128 


IX R 


BUMP POINTER 


0107 

004C 

F60A 


JMP 

0T1 i 

L00P 

OIOS 

004D 

GOOO 

CAI 

CAI 

o : 


0109 







0110 



* OUTPUT FLOATING POINT 

NUMBER 

0111 



* 




Ollfi 

004E 

03 00 

0FPA 

ENT 



0113 

004F 

E70 1 


UK 

#0FPA 

GET POINTER 

0114 

0050 

DE02 


IMS 

0FPA 

SET RETURN ADDRESS 

0115 

0051 

EA01 


STX 

0PT 

SAVE POINTER 

01 16 

0052 

F30E 


J ST 

* FAS 

CONVERT' T0 ASCII 

0117 

0053 

0000 

0PT 

DATA 

0 


OHS 

0054 

0059 


DATA 

BUF 


0119 

0055 

0110 


EAR 


SET HID FLAG 

01 SO 

0056 

FEi A 


J ST 

0TL 

PRINT NUMBER 

0121 

0057 

0059 


DATA 

BUF 


0122 

0058 

F70A 


RTM 

0FPA 


0123 

0059 

0000 

BUF 

RES 

S^ 0 


0124 

0061 


FAS 

REF 



0125 

0000 

ERRORS 



HID 




153 


0001 




NAM 

READ., CL 

EAR 

r 

0002 




NAM 

FUN CH,» GRAPH 

\ 

s ! 

0003 




EXTR 

IKB* 0TT., El PT 


0004 




EXTR 

0TL, CRLF^ FPL0T, 0FPA 

11 j 

i 

0005 




EXTR 

XAjXBjX 

C> FAD 

Y* 

0006 

0000 



REL 

0 


‘ j 

0007 



* 





0003 



* READ PAPER 

TAPE & ADD T0 BUFFER 

I 

0009 







| . 

0010 

0000 

08 00 

READ 

ENT 



I 

001 1 



* INITIALIZE VARIABLES 


L 

0012 

0001 

1 128 


RLX ' 

r 

SET INDEX BIT 

i 

0013 

0002 

1400 


S0V 



j 

0014 

0003 

1 IAS 


REX 

1 


is 

i 

0015 

0004 

EA33 


STX 

R5 

& SAVE PB INTER 

I 

0016 

0005 

EA34 


STX 

R5+2 


1 

% 

0017 

0006 

0523 


XRP 


RESET SU3SCRI FTS 

a 

t 

a 

0018 

0007 

E9 00 


STX 

XA 


] 

’ 


0000 





i 

0019 

0008 

E9 00 


STX 

XB 


! 

i 

* 


0000 





i 

0020 

0009 

E9 00 


STX 

xc 


] ■ 



0000 





5 

0021 

000A 

B236 


LDA 

TP 

SET INPUT BUFFER POINTER 


0022 

OOOB 

9 A38 


STA 

MPT 



0023 

OOOC 

5235 


LDA 

CT 

SET INPUT C0UNT 

: : -- 

0024 

OOOD 

9A37 


STA 

GNT 


I 

0025 



*SKIP 

LEADER 

& LABEL 


] 

i 

0026 

OOOE 

C63 A 


LAP 

:8A 

LINE FEED 

i 

, t 

0027 

OOOF 

F9 00 


J ST 

BTT 


\ 

t ; 



0000 





« ! 
i i 

1 .. 

0028 

0010 

F9 00 


J ST 

BIPT 


t 

1- . 



0000 





j 

0029 

001 1 

2141 


JAZ 

5- 1 

SKIP LEADER 

i : - 

0030 

0012 

F9O0 

R1 

J ST 

GTT 

ECH0 LABEL 




0000 






0031 

0013 

F9 00 


J ST 

BIPT 

READ TAPE 




0000 






0032 

0014 

C09 2 


CAI 

: 92 

END 0 F LABEL? 


0033 

0015 

F20 1 


JMP 

R2 

YES 


0034 

0016 

F604 


JMP 

R1 

N0 


0035 

0017 

F9 00 

R2 

J ST 

BIPT 

READ TAPE 




0000 






0036 

0018 

GOFF 


CAI 

: FF 

FILE MARK? 


0037 

0019 

F20 1 


JMP 

R3 

YES" 


0038 

00 1 A 

F603 


JMP 

R2 

N0 


0039 



READ 

L00P 


• 


0040 

001B 

F9 00 

R3 

J ST 

BI PT 

LOWER BYTE 


• 


0000 





j 

0041 

001C 

IBS 7 


LLR 

8 

SAVE 

s 

I 

0042 

00 ID 

F9 00 


J ST 

BIPT 

UPPER BYTE 




0000 





j 

0043 

001 E 

1B07 


LLL 

8 

REST0RE 

■ i 

si 

I 

1! 


i 


154 


0044 

00 IF 

£224 


LEX 

M PT ■ 

SAVE LOW BITS 

0045 

0020 

9 CO 1 


STA 

§1 : 


0046 

0021 

F9 00 
0000 


J ST 

BIPT 

LOWER BYTE 

0047 

0022 

IBS 7 


LLR 

8 


0048 

0023 

F9 00 
0000 


J ST 

BIPT 

UPPER BYTE 

0049 

0024 

1B07 


LLL 

8 


0050 



* CONVERT BASIC- FoP. TO 

CAI-F. P. 

0051 

0025 

2109 


J AZ 

zr ■' 

0=0=0' “ 

0052 

0026 

0048 


rr% 
1 . 


SAVE SIGN 

0053 

0027 

303 1 


JAP 

S+ 2 

ABS. VAL. 

0054 

0028 

0310 


NAR 


* * 

0055 

0029 

IBS 7 


LLR 

8 

REMOVE MSB 

0056 

002A 

1328 


LLX 

1 


0057 

0025 

SA17 


ADD 

D64 

FIX CHARACTERISTIC 

0058 

002C 

1B07 


LLL 

8 


0059 

002D 

13C0 


LAS 


RECOVER SIGN 

0060 

0G2E 

11 DO 


RRA 

1 


0061 

002F 

9514 

Z 1 

STA 

*MPT 

SAVE HI BITS 

0062 

0030 

DAl 3 


IMS 

MPT 

BUMP POINTER 

0063 

0031 

DAI 2 


IMS 

MPT 

TWICE! 

0064 

0032 

DAI 2 


IMS 

CMT 

MS RE? ' 

0065 

0033 

F61S 


JMP 

R3 

YES 

0066 



* ADD INPUT BUFFER T0 BEAM ARRAY 

0067 



* PAN I C 

SWITCH DISABLED 

FOR DURATION 

0068 

0034 

4006 


Cl D 

* 

DISABLE AUTO 

0069 

0035 

B20C 


L DA 

CT 

RESET COUNTER 

0070 

0036 

9 AOE 


STA 

CNT 


0071 

0037 

F9 00 

R4 

J ST 

FAD 

F. P* ADD 



0000 




' ■ 1 * 

0072 

0038 

0000 

R5 

DATA 

Oj : 9400.* 0 


0039 

9400 






003A 

0000 





0073 

0035 

D9 00 
0000 


IMS 

XA 

BUMP SUBSCRIPTS 

0074 

003C 

D900 

0000 


IMS 

XB 


0075 

003D 

D9 00 
0000 


IMS 

XC 


0076 

003E 

DAO 6 


IMS 

CNT 

MORE? 

0077 

003F 

F608 


JMP 

R4 

YES 

0078 

0040 

F740 


RTM 

READ 

* 

0079 



* STORAGE 



008 0 

0041 

1400 

TP 

DATA 

: 1400 


008 1 

0042 

FFOO 

CT 

DATA 

-256 


008 2 

0043 

0040 

D64 

DATA 

’64 


008 3 

0044 

0000 

MPT 

DATA 

0 


0334 

0045 

0000 

CMT 

DATA 

0 


OCGS 



* 




008 6 



CL* EAR 

BEAM 

ARRAY 


0037 



* 


• 


0088 

0046 

08 0 0 

CLEAR 

ENT 




00B9 

0047 

B605 


L DA 

CT : 

SET C0UMT 

009 0 

0048 

9E03 


STA 

CNT ' 


009 1 

0049 

0110 


EAR 


CLEAR 

009 2 

004A 

9 COO 

Cl 

STA 

GO 

HI BITS 

009 3 

004B 

9 CO 1 


STA 

©1 

L0 BITS 

0094 

004C 

C202 


AX I 

2 ’ 

BUMP POINTER 

009 5 

004D 

DE08 


IMS 

CNT 

DSN E? 

009 6 

004E 

F 604 


JMP 

Cl 

N0 

009 7 

004F 

F709 


RTN 

CLEAR 


0093 



* 




0099 



* PUNCH 

BEAM 

ARRAY - 

BASIC FORMAT 

0100 



* 




OiOI 

0050 

08 00 

PUNCH 

ENT 



0102 

0051 

FA22 


J ST 

LEAD 

L EADER 

0103 

0052 

B610 


L DA 

CT 

START COUNTER ^ 

0104 

0053 

9 EOE 


STA 

CNT 


0105 

0054 

F9 00 

PI 

J ST 

IK3 

ECHO LABEL 



0000 





0106 

0055 

009 2 


CAI 

: 92 

CTRL /TAPE? 

0107 

0056 

F201 


JMP 

P2 

YES *' 

0.1 OS 

0057 

F603 


JMP 

PI 

N0 

0109 

0053 

C6FF 

P2 

LAP 

s FF 

PUNCH FILE MARK 

0110 

0059 

F9 00 


J ST 

0TT 




0000 





01 1 1 

005A 

B40 1 

P3 

LDA 

©1 

PUNCH L0 BITS 

0112 

005B 

F9 00 


J ST 

0TT 




0000 





0113 

0050 

13P7 


LRA 

8 


01 14 

005D 

F9 00 


J ST 

0TT 




0000 





01 15 

005E 

B400 


LDA 

©0 

GET HI BITS 

0116 

005F 

C202 


AKI 

2 

BUMP POINTER 

01 17 

0060 

EE1C 


STX 

MPT 

AN D SAVE 

01 IS 



* CONVERT CAI 

-F. P. T0 

BASIC- F. P. 

0119 

0061 

210A 


JAZ 

‘ Z2 ' 

C= 0= 0 ' ’ 

0120 

0062 

0048 


TAX 


SAVE SIGN 

0121 

0063 

1350 


LLA 

1 

CLEAR A15 

0122 

0064 

IBS 7 


LLR 

8 

SPLIT 

0123 

0065 

9 622 


SUB 

D64 

FIX CHARACTERISTIC 

0124 

0066 

1400 


SO V 


INSERT MSB 

0125 

0067 

1150 


RL A 

1 


0126 

0068 

1B06 


LLL 

7 

RE- FORMAT 

0127 

0069 

1329 


LLX 

2 

RECOVER SIGN 

0128 

006A 

3201 


J0R 

£+2 

CORRECT FOR IT 

0129 

006B 

0310 


N AR 



0130 

006C 

F9 00 

22 

J ST 

0TT 

PUN CH HI BITS 



0000 





0131 

006D 

1B87 


LLR 

3 


0l$2 

006E 

F9 00 


J ST 

0TT 




0000 





0133 

006F 

E62B 


LEX 

MPT 

RECOVER POINTER 

0134 

007 0 

DE2B 


IMS 

‘CNT 

DON E? ‘ 

0135 

0071 

F617 


JMP 

P3 

NO 


0136 

0072 

FA01 


J ST 

LEAD 

YESj punch leader 

0137 

0073 

F723 


RTN 

PUNCH 


0133 



* PINCH 

S'* 0F 

L EADER 


0139 

0074 

03 00 

LEAD 

EMT 



0U0 

0075 

C732 


LAM 

50 


OUll 

0076 

9E31 


STA 

CNT 


0142 

007 7 

0110 


EAR 



0143 

0073 

F9 00 
0000 

L2 

J ST 

0TT 


0144 

0079 

DE34 


IMS 

CNT 


0145 

007A 

F602 


JMP 

L 2 


0146 

0147 

007B 

F707 

* 

RTN 

LEAD 


0143 



*PL0T 

DATA ARRAY 


0149 

0150 

0070 

03 00 

* 

GRAPH 

EMT 



0151 

007D 

EA06 


SIX 

T 

SAVE C0UNT POINTER 

0152 

007E 

C202 


AKI 

2 


0153 

007F 

EAOS 


sta 

PT 

SAVE DATA POINTER 

0154 

008 0 

0110 


EAR 


PRINT COUNT 

0155 

003 1 

F9 00 
0000 


J ST 

0TL 


0156 

003 2 

OOSA 


DATA 

CTO 


0157 

003 3 

F9 00 
0000 


J ST 

0 FPA 


015S 

008 4 

0000 

T 

DATA 

0 


0159 

008 5 

F9.00 

UOOO 


J ST 

CRLF 


0160 

003 6 

C7FF 


LAM 

255 

PLOT DATA 

0161 

003 7 

F9 00 


J ST 

FPL0T 


‘ 


0000 





0162 

0038 

0000 

PT 

DATA 

0 


0163 

0089 

F70D 


RTN 

GRAPH 


0164 

OOSA 

8 A3 A 

CTO 

DATA 

: 8 A3 A 


0165 

008 B 

C3CF 


TEXT 

’ C0UNT 



008 C 

£)5CE 






008 D 

D4A0 






008 E 

BDAO 





0166 

008 F 

0000 


DATA 

0 


Ol 67 
0000 

ERRORS 



Ew D 




0001 

0002 

0003 

0004 
000 5 
0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
0027 
0023 

0029 

0030 

0031 
•m32 

0033 

0034 

0035 

0036 

0037 
0033 


0039 
004 0 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 


& 

*• 

* 

* 

* 

* 

* 


FH-0T - FLOATING POINT PLOTTER 


CALLING SEQUENCE: 


L DA MDIM 
J ST *FPL0T 
DATA ARRAY 
( RETURN ) 





NAM 

FPL0T 

0000 



REL 

0 

0000 

08 00 

FFL0T 

ENT 


0001 

9 A4C 


STA 

NR0U 

0002 

9 A4C 


STA 

CNT 

0003 

B703 


L DA 

*FPL0T 

0004 

9 Air 


STA 

RP3 

0005 

9 A24 


STA 

RF4 

0006 

A249 


I0R 

BIS 

0007 

9A05 


STA 

RP1 

0005 

9A03 


STA 

RP2 

0009 

0350 


ARP 


OOOA 

9B4F 


STA 

*XA 

0003 

F204 


JMP 

M0VE 

OOOC 

F34E 

LP1 

J ST 

* FCP 

OOOD 

0000 

RP1 

DATA 

0; MAX 

OOOE 

0051 




000F 

2183 


UAL 

INC 

0010 

FB43 

HOVE 

J ST 

*FMV 

0011 

0000 

RP2 

DATA 

0#MAK 

0012 

0051 




0013 

D346 

INC 

IMS 

*)CA 

0014 

DA3A 


IMS 

GIT 

0015 

F609 


JMP 

LP1 

0016 

FE41 


J ST 

*CRLF 

0017 

FB40 


J ST 

* CRLF 

0018 

C6BD 


LAP 

i— » 

0019 

FB3F 


J ST 

£0TL 

001 A 

0061 


DATA 

SF 

OOlB 

FB43 


J ST 

*FDV 

001C 

0055 


DATA 

FSO^MAX^ SCALE 

00 ID 

0051 




001E 

0053 




00 IF 

F340 


J ST 

*0FPA 

0020 

0053 


data 

SCALE 

0021 

FB36 


J ST 

* CRLF 

0022 

FB35 


JST 

* CRLF 

0023 

FB3C 

LP2 

J ST 

*0FPA 

0024 

0000 

RP3 

DATA 

0 

0025 

DEO 1 


IMS 

RP3 

0026 

DEO 2 


IMS 

RP3 

0027 

C6A0 


LAP 

. » « 

0028 

FB2E 


JST 

• £0TT 

0029 

FB34 


JST 

*FMP 


158 


0050 

002A 

0000 

RP4 

DATA 

0, SCALE* MAX 


002B 

0053 





002C 

0051 




0051 

002D 

DE03 


IMS 

RP4 

0052 

002E 

DE04 


IMS 

RP4 

0053 

002F 

FB2D 


j ST 

*FIX 

0054 

0030 

0051 


DATA 

MAX* CNT 


0031 

004F 




0055 

0032 

B21C 


L DA 

CNT 

0056 

0033 

0048 


TAX 


0057 

0034 

305 5 


JAP 

P0S 

0058 

0035 

C6AD 


LAP 

>. i 

0059 

0036 

F320 


J ST 

£'0TT 

0060 

0037 

C6B0 


LAP 

*0* 

0061 

0033 

FB1E 


J ST 

£0TT 

0062 

0039 

F20E 


JMP 

OUT 

0063 

003A 

C6A0 

P0S 

LAP 

I f 

0064 

003B 

F313 


J ST 

£0TT 

0065 

003C 

23 09 


JXZ 

MARK 

0066 

003D 

C6B0 


LAP 

• o« 

0067 

003E 

FB IS 


J ST 

*0TT 

0063 

003F 

0503 


MX R 


0069 

0040 

EAOE 


STX 

CNT 

0070 

0041 

C6AA 


LAP 

<* • 

0071 

0042 

F201 


JMP 

$+2 

0072 

0043 

FBI 3 


J ST 

*0TT 

0073 • 

0044 

DAOA 


IMS 

CNT 

0074 

0045 

F 602 


JMP 

S-2 

0075 

0046 

C6AA 

MARK 

LAP 

* 

0076 

0047 

F30F 


J ST 

£0TT 

0077 

0048 

FBOF 

0UT 

J ST 

# CEL F 

0078 

0049 

DAO 4 


IMS 

MR0W 

0079 

004A 

F 627 


JMP 

LP2 

008 0 

004B 

FBOC 


J ST 

CRLF 

008 1 

0040 

DE4C 


IMS 

FPL0T 

003 2 

004D 

F74D 


KTM 

FPL0T 

008 3 

004E 

08 00 

NR0W 

KLT 


008 4 

004F 

03 00 

CM T 

HLT 


003 5 

0050 

8000 

B15 

DATA 

s 8 000 

008 6 

0051 

0000 

MAX 

RES 

2* 0 

003 7 

0053 

0000 

SCALE 

RES 

2* 0 

0088 

0055 

4348 

F50 

DATA 

: 4348* 0 


0056 

0000 




0089 

0057 


0TT 

REF 


009 0 

0058 


CRLF 

REF 


009 1 

0059 


0TL 

REF 


009 2 

005A 


XA 

REF 


009 3 

0053 


FCP 

REF 


009 4 

005C 


FI-i V 

REF 


009 5 

005D 


FIX 

REF 


009 6 

005E 


FMP 

REF 


009 7 

005F 


FDV 

REF 

* 

0098 

0060 


0FPA 

REF 



0099 0061 D3C3 SF TEXT 

0062 C1CC 

0063 S5C6 

0064 C1C3 

0065 D4CF 

0066 D2A0 

0067 BDAO 


*S CM. FACTOR =* 


0100 


EN D 


0001 




0002 



* CRT 

0003 



* 

0004 



* 

0005 



* 

0006 



* 

0007 



* 

0008 




0009 

0000 



0010 

0000 

08 00 

CRT 

0011 

0001 

C7FF 


0012 

0002 

9A5F 


0013 

0003 

9A5F 


0014 

0004 

49 C4 


0015 

0005 

F601 


0016 

0006 

40C7 


0017 

0007 

C202 


0018 

00C8 

0030 


0019 

0009 

A261 


00.20 

000A 

9 AOB 


0021 

OOOB 

9 AOE 


0022 

OOOC 

9A10 


0023 

OOOD 

9A13 


0024 

OOOE 

9A36 


0025 

OOOF 

0110 


0026 

0010 

9A54 


0027 

001 1 

9A54 


0028 

0012 

0350 


0029 

0013 

9 55B 


0030 

0014 

F204 


0031 

0015 

FB5A 

LP1 

00 32 

0016 

0000 

RP1 


0017 

0067 


0033 

0018 

218 3 


0034 

0019 

FB57 

M0VE 

0035 

00 1A 

0000 

RP2 


00 IB 

0067 


0036 

00 1C 

FB53 

INC1 

0037 

00 1 D 

0000 

RP3 


00 1 E 

0065 


0038 

00 IF 

308 3 


0039 

0020 

FB50 


0040 

0021 

0000 

RP4 


0022 

0065 


0041 

0023 

DB4B 

INC2 

0042 

0024 

DA3E 


0043 

0025 

F610 


0044 



* 

00*45 

0026 

FB4F 


0046 

0027 

FB4E 


0047 

0028 

C6BD 


0048 

0029 

FB4D 


0049 

002A 

0079 



FL SATING POINT PL0TTER 

CALLING SEQUENCE: 

J ST sfc CRT 
< RETURN ) 


NAM 

CRT 

REL 

0 

ENT 


LAM 

255 

STA 

N R0 W 

STA 

CNT 

SEN 

24* 4 

JMP 

S- 1 

SEL 

24* 7 

AXI 

2 ' 

TXA 


I0R 

B 1 5 

STA 

RP1 

STA 

RP2 

STA 

RP3 

STA 

RP4 

STA 

RP5 

ZAR 


STA 

M IN 

STA 

HIN+ 1 

ARP 


STA 

*XA 

JMP 

MB VE 

J ST 

*FCP 

DATA 

0*MAX 

JAL 

IN Cl 

J ST 

*FMV 

DATA 

0*MAX 

J ST 

*FCP 

DATA 

0*MIN 

JAP 

INC 2 

J ST 

*FMV 

DATA 

0* M I N 

IMS 

*XA 

IMS 

CNT 

JMP 

LP1 

J ST 

*CRLF 

J ST 

*CRLF 

LAP 

i — i 

J ST 

**0TL 

DATA 

SF 


160 


0050 

0051 

0052 


0053 

0054 


0055 

0056 

0057 

0058 

0059 

0060 


0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 


0073 

0074 

0075 

0076 

0077 

0078 

0079 
008 0 
008 1 
0082 
008 3 
003 4 

008 5 
008 6 
0087 
0038 
0089 

009 0 
009 1 
009 2 


* 


002B 

FB46 


J ST 

*FSE 

002C 

0067 


DATA 

M AX*MIN* MAX 

002D 

0065 




002E 

0067 




002F 

FE43 


J ST 

*FDV 

0030 

006C 


DATA 

F250*MAX* SCALE 

0031 

0067 




0032 

0069 




0033 

FB44 


J ST 

*0FPA 

0034 

0069 


DATA 

SCALE 

0035 

FB40 


J ST 

sfcCHLF 

0036 

FB3D 

V 

J ST 

♦ IMF 

0037 

0069 


DATA 

SCALE* MIN; MAX 

0038 

0065 




0039 

0067 




003A 

F33A 


J ST 

* FIX 

00333 

0067 


DATA 

MAX* MIN 

003C 

0065 




003D 

B227 


L DA 

MIN 

003E 

0310 


NAR 


003F 

308 1 


JAP 

S+ 2 

0040 

0110 


Z AH 


0041 

9A22 


STA 

X 

0042 

0350 


ARP 


0043 

9B2B 

LP2 

STA 

*XA 

0044 

FB2F 

J ST 

* FMP 

0045 

0000 

RP5 

DATA 

0* S CAL E* M AX 

0046 

0069 




0047 

0067 




0048 

DB26 


IMS 

*XA 

0049 

FB2B 


J ST 

*FIX 

004A 

0067 


DATA 

MAX* CNT 

004B 

0063 




004C 

E2 1 7 


LHC 

X 

004D 

B215 


L DA 

CNT 

004E 

49 C4 

D0T 

SEN 

24* 4 

004F 

F60 1 


JMP 

S- 1 

0050 

49 C2 


SEN 

24* 2 

0051 

r 601 


JMP 

S- 1 

0052 

6EC2 


01X 

24* 2 

0053 

2107 


JAZ 

D0NE 

0054 

318 3 


JAG 

P0S 

0055 

00A3 


EKR 


0056 

0150 


I AH 


0057 

F609 


JMP 

D3T 

0058 

0128 

P0S 

IXR 


0059 

00D0 


DAR 


005A 

F60C 


JMP 

•D0T 

0053 

3208 

D0NE 

L DA 

• V 
Jk 

005C 

SAIT 


ADD 

H256 


009 3 

005D 

9A06 


009 4 

005E 

DAO 3 


009 5 

005F 

F61B 


009 6 

0060 

FBI 5 


009 7 

0061 

F761 


0098 




0099 

0062 

0000 

NR0W 

0100 

0063 

0000 

GMT 

0101 

0064 

0000 

X 

0102 

0065 

0000 

MIN 

0103 

0067 

0000 

MAX 

0104 

0069 

0000 

SCALE 

0105 

006B 

8 000 

B 1 5 

0106 

006C 

447A 

F250 


006D 

0000 


0107 

006E 

0100 

K256 

0103 


' 

Si 

0109 

006F 


XA 

0110 

0070 


FCP 

01 1 1 

007 1 


FHV 

0112 

0072 


FSB 

0113 

0073 


FDV 

01 U 

0074 


FMP 

0115 

0075 


FIX 

0116 

0076 


CRL F 

0117 

0077 


GTL 

01 13 

0073 


0FPA 

01 19 



* 

0120 

0079 

D3C3 

SF 


007 A 

C1CC 



0073 

C5A0 



007C 

C6C1 



007D 

C3D4 



007E 

CFD2 



007F 

AOBD 


0121 



* 

0122 




0000 

ERRORS 




STA X 

IMS NR0W 
JMP LP2 
J ST * CRL F 
RTM CRT 

DATA 0 

DATA 0 

DATA 0 

RES a. 0 

RES 2 .» 0 

RES a, 0 

DATA : 8 000 

DATA : 447A, 0 

DATA 256 

REF 

REF 

REF 

REF 

REF 

REF 

REF 

REF 

REF 

REF 

TEXT ’SCALE FACTS R 


END 


0001 




i J AM 

0002 




EXTR 

0003 




MAM 

0004 




EXTR 

0005 

0000 

03 00 

VI EIJ 

ENT 

0006 

0001 

9A7 1 

DI SP 

STA 

0007 

0002 

C202 


AXI 

0008 

0003 

EA66 


STA 

0009 

0004 

EA66 


STX 

0010 

0005 

C7FF 


LAM 

0011 

0006 

9A6A 


STA 

0012 

0007 

9A64 


STA 

0013 

0008 

0110 


ZAR 

0014 

Uuu9 

9A63 


-TA 

0015 

000A 

9 A63 


STA 

0016 

OOOB 

9 A63 


STA 

0017 

OOOC 

9 A63 


STA 

0018 

OOOD 

8400 

ST CM 

LDA 

0019 

0002 

D260 


CMS 

0020 

oo of 

F20A 


JMP 

0021 

0010 

F205 


JMP 

0022 

0011 

B401 


LDA 

0023 

0012 

D25D 


CMS 

0024 

0013 

F20S 


JMP 

0025 

0014 

9A53 


STA 

0026 

0015 

F21 1 


JMP 

0027 

0016 

9A5S 

MAX 1 

STA 

0028 

0017 

8401 


LDA 

0029 

0018 

9A57 


STA 

0030 

0019 

F20D 


JMP 

0031 

001A 

3400 

MINI 

LDA 

0032 

0018 

025 1 


CMS 

0033 

001C 

F206 


JMP 

0034 

001D 

F209 


JMP 

0035 

001E 

8401 


LDA 

0036 

001F 

D24E 


CMS 

0037 

0020 

F202 


JMP 

0038 

0021 

F205 


JMP 

0039 

0022 

F204 


JMP 

0040 

0023 

8400 

MIN 2 

LDA 

0041 

0024 , 

, 9A48 
' B401 

' 

STA 

0042 

0025 


LDA 

004 3 

0026 

9A47 


STA 

0044 

0027 

E242 

FI CM 

LDX 

0045 

0028 

0202 


AXI 

0046 

0029 

EA40 


STX 

0047 

002A 

DA46 


IM S 

0048 

002B 

F61E 


JMP 

0049 



* 


0050 

002C 

B242 


LDA 

0051 

002D 

E242 


LDX 

0052 

002E 

F9 00 
0000 


J ST 


VI EV 



DPW : > DP SUB: , DPDI V: .» DPSH: 
DI SP 

RTTRj RT0P 

VEVR 

2 


RP1D 

RP2D 

255 

CNT 

POINTER 

NR0W 




in I TI AL 

IZE MAX AMD MIN 

MIN 
MIN+ 1 

max 

MAX* 1 
@0 
MAX 



MIN 1 

AKMAX 


M AX 1 

Al>MAX 


el 

Al=MAX 


MAX+ 1 



MIN 1 ’ 

A1=MAX 

A2<MAX+ 1 

MAX+ 1 

Al=MAX 

A2>MAX+ 1 

FI CM ‘ 
MAX 

§1 

MAX* l 
FI CM ‘ 

eo 

MIN 


• 

MIN 2 

AKMIN 


FI CM 

A1>MIN 


@1 

Al=MIN 


MIN+ 1 



MIN 2 

A1=MIN 

A2<MIN* 1 

FI Cl 

Al=MIN 

A2>MIN+ 1 

FI CM 

Al-MIN 

A2=MIN+ 1 

GO 

MIN 

@1 

M1W+ 1 
RP1D 



2 * 
RP1D 

*• 


CNT 
ST CM 

MAX 
MAX* 1 
DPSUB: 





1 


i 

0053 

002F 

006D 


DATA 

MIN : 

j 

0054 

0030 

1050 


ALA 

1 

* 

0055 

0031 

1329 


LLX 

2 

1. 

0056 

0032 

3201 


J0R 

S+2 

r 

i 

i 

0057 

0033 

0150 


IAR 


i; 

00 58 

0034 

13A3 


LRX 

1 

!; 

0059 

0035 

9 A39 


STA 

Max 

v L 

0060 

0036 

EA39 


srx 

MAX* 1 

1 

i 

0061 






1 

} 

0062 

0037 

B235 


L DA 

MIN 

j* 

0063 

0038 

E235 


LEX 

MIW+ 1 

1 

0064 

0039 

F9 00 


J ST 

DPDI V: 

i 



0000 




| 

0065 

003A 

006F 


DATA 

MAX 

i 

( 

0066 

003B 

9A2D 


STA 

TEM 

i 

0067 

0030 

1 0D1 


AHA 

2 

i 

0063 

003D 

0310 


WAR 


j 

0069 

003E 

S A2A 


ADD 

TEH 

\ 

0070 

003F 

10D5 


AHA 

6 

1 

\ 

0071 



* 



\ 

0072 

0040 

0310 


WAR 


[ 

0073 

0041 

303 1 


JAP 

£+ 2 

(■ 

0074 

0042 

0110 


ZAR 


i 

0075 

0043 

9A2E 


STA 

X 

i 

0076 






i 

0077 

0044 

E226 


L EEC 

RP2D 

\ 

0073 ' 

0045 

B400 

LPD 

L DA 

GO 

1 

0079 

0046 

E40 1 


LEK 

§1 

i 

003 0 

0047 

F9 00 


J ST 

DPDI V: 

\ 

i 



0000 




i - 

003 1 

0048 

006F 


DATA 

MAX 

j 

003 2 

0049 

9 A1F 


STA 

TEM 

i 

i ' 

003 3 

004A 

10D1 


ARA 

2 


003 4 

004B 

0310 


NAR 



003 5 

004C 

SA1C 


ADD 

TEM 


003 6 

Q04D 

10D5 


ARA 

6 


003 7 







0038 

004E 

£223 


LEK 

X 


0039 

004F 

49 C4 

DOT 

SEN 

24j 4 


009 0 

0050 

F60 1 


JMP 

S- 1 


009 1 

0051 

49 02 


SEN 

2 4; 2 


009 2 

0052 

F601 


JMP 

S- 1 


009 3 

0053 

6EC2 


'BIX 

24, 2 


009 4 

0054 

2107 


JAZ 

DONE 


009 5 

0055 

318 3 


JAG 

P0S 

$ 

It 

009 6 

0056 

00 A3 


DXR 


ji ’ 

0097 

0057 

0150 


IAR 


1 

y 

0098 

0053 

F609 


JMP 

DST 

S ' 

0099 

0059 

0123 

P0S 

IX R ’ 


n 

0100 

005A 

OODO 


DAR 


I ' 

0101 

005B 

F60G 


JMP 

D0T 

1 - 

0102 



* 


* 

a 

0103 

005C 

321 5 

DOME 

LDA 

> 



164 



0104 

005D 

3A16 


ADD 

H256 

0105 

005E 

9A13 


STA 

ic : 

0106 

005F 

E20B 


L EX 

RP2D 

0107 

0060 

C202 


AXI 

2 

OIOS 

0061 

EA09 


STA 

RP2D 

0109 

0062 

DA09 


IMS 

NR0U 

0110 

0063 

F61E 


JMP 

LPD 

0111 

0064 

B20E 


LDA 

VEVR 

01 12 

0065 

2101 


J AS 

S+P 

0113 

0066 

FlOO 

0000 


JMP 

RT0P 

0114 

0067 

FlOO 

0000 


JMP 

RTTR 

0115 

0068 

F768 


RTN 

VIEW 

0116 

0069 

0000 

THM 

DATA 

0 

0117 

006A 

0000 

RP1D 

DATA 

0 

01 18 

006B 

0000 

RP2D 

DATA 

0 

0119 

006C 

0000 

NR0V 

DATA 

0 

0120 

006D 

0000 

MIN 

RES 

2j 0 

0121 

006F 

0000 

MAX 

RES 

2j 0 

0122 

0071 

0000 

CNT 

DATA 

0 

0123 

0072 

0000 

X 

DATA 

0 

0124 

0073 

0000 

VEVR 

DATA 

0 

0125 

0126 

0074 

0100 

H256 

DATA 

END 

256 


0000 ERRORS 


0001 




NAM 

DPACC 

0002 




N AM 

DPFLT 

0003 




NAM 

DPFXX 

0004 




EXTR 

DPi\l: 

0005 




EX TP. 

DPM 

0006 

0000 



BEL 

0 

0007 



* 



0008 



* DOUBLE PRECISION ACCUMUL 

0009 



❖ 



0010 

0000 

08 00 

DPACC 

ENT 


0011 

0001 

9 Ai 0 


STA 

TEMP 

0012 

0002 

0030 


TXA 


0013 

0003 

E703 


LDX 

❖ DP AC C 

0014 

0004 

1200 


ROV 


0015 

0005 

8 C01 


ADD 

61 

0016 

0006 

S20C 


AND 

MASK 

0017 

0007 

9 COl 


STA 

@1 

0018 

0008 

B209. 


LDA 

TEMP 

0019 

0009 

3204 


JOR 

DC 1 

0020 

000A 

1200 


ROV 

* 

0021 

OOOB 

0150 


IAR 


0022 

OOOC 

3201 


JOR 

DC 1 

0023 

OOOD 

0110 


EAR 


0024 

OOOB 

8 000 

DC 1 

ADD 

60 

00' 5 

OOOP 

9 COO 


STA 

60 

0026 

0010 

DE10 


I'M S 

DPACC 

0027 

0011 

F71 1 


RTLv 

DPACC 

0028 

0012 

OODO 

TEMP 

DATA 

0 - ' 

0029 

0013 

7FFF 

MASK 

DATA 

: 7F ft 

'i 

sf C 



* 



0031 



❖ DOUBLE PRECI 

SION TO CAI- 

0032 



❖ 


“ - " * 

0033 

0014 

08 00 

DPFLT 

ENT 


0034 

0015 

38 01 


JXN 

DPF1 

0035 

0016 

2113 


J AZ 

DPF4 

0036 

0017 

BA 13 

DPF1 

EM A 

BITS 

0037 

0018 

0110 


ZAP. 

• 

0038 

0019 

BA1 1 


EM A 

BITS 

0039 

001 A 

9 Al 1 


STA 

SI GN 

0040 

001B 

308 1 


JAP 

DPF2 

0041 

001 C 

F9 00 


J ST 

DPI'!: 



0000 




0042 

001D 

1328 

DPF2 

LLX 

1 

0043 

001E 

1B00- 

DPF3 

LLL 

1 

0044 

00 IF 

DAOB 


IMS 

BITS 

0045 

0020 

3242 


JfiR 

DPF3 

0046 

0021 

IBS 8 


LLR 

9 

0047 

0022 

BA08 


EM A 

BITS 

0048 

0023 

0310 


WAR 


0049 

0024 

8 AOS 


ADD 

D160 

0050 

0025 

1356 


LLA 

7 

0051 

0026 

A204 


I OR 

•BITS 

0052 

0027 

BA04 


EM A 

SIGN 


167 


0053 

0028 

3205 


AND 

DPi-1 ' 

0054 

0029 

A202 


I SR 

SIGN 

0055 

002A 

F716 

DPF4 

RTN 

DPFLT 

0056 

0023 

0000 

BITS 

DATA 

0 

0057 

002C 

0000 

SI GW 

DATA 

0 

0053 

002D 

OOAO 

D 1 60 

DATA 

160 

0059 

002E 

3000 

DPH 

DATA 

: 8 000 

0060 

002F 

OS 00 

DPFIX 

ENT 


0061 

0030 

9E04 


STA 

SIGN 

0062 

0031 

8 2 1 C 


AMD 

DPP 

0063 

0032 

1 3D6 


LRA 

7 

0064 

0033 

9606 


SUB 

D 1 60 

0065 

0034 

9E09 


STA 

BITS 

0066 

0035 

B609 


L DA 

SIGN 

0067 

0036 

8213 


AND 

M AN T 

0063 

0037 

A218 


I0R 

SBIT 

0069 

0033 

3E0D 


EM A 

BITS 

0070 

0039 

2109 


JAZ 

DPFX2 

0071 

00 3A 

319 0 


JAG 

DPFX4 

0072 

003B 

SE10 


EM A 

BITS 

007 3 

003C 

1B07 


LLL 

B 

0074 

003D 

1B30 

DPFX 1 

LLR 

1 

0075 

003E 

DEI 3 


IMS 

BITS 

0076 

003F 

F602 


JMP 

DFFX 1 

0077 

0040 

1B00 


LLL 

1 

0073 

0041 

13A3 


LEX 

1 

0079 

0042 

F201 


JMP 

S+2 

008 0 

0043 

BE18 

DPFX 2 

EM A 

BITS 

OOS 1 

0044 

BE IS 


EM A 

SIGN 

003 2 

0045 

303 3 


JAP 

DPFX3 

008 3 

0046 

BE1 A 


EM A 

SI GN 

003 4 

0047 

F9 00 


J ST 

DPN; 



0000 




008 5 

0048 

F201 


JMP 

$+2 

008 6 

0049 

3E1D 

DPFX 3 

EH A 

SI GN 

003 7 

004A 

F71B 


RTN 

DPFIX 

0033 

004B 

0113 

DPFX 4 

Z AX 


0039 

0040 

1400 


S3 V 


009 0 

004D 

F71E 


RTN 

DPFIX 

009 1 

004E 

7FFF 

DPP 

DATA 

: 7FFF 

009 2 

004F 

007F 

M AN T 

DATA 

: 7F 

009 3 

0050 

003 0 

SBIT 

DATA 

: 8 0 

009 4 




END 


0000 

ERRORS 






168 A 


APPENDIX C 


0001 




NAM 

CMN D 


0002 




EXTR 

I ERj'BTTj STL, CRLF 

0003 




EXTR 

ERR,XEQ 


0004 

0100 



ABS 

; 100 


0005 



*PRINT 

N AM E 



0006 

0100 

0110 


ZAR 



0007 

0101 

F9 00 
0000 


J ST 

STL 


0003 

0102 

012B 


DATA 

NAME 


0009 

0010 

0103 

F201 

* 

JMP 

S+2 


0011 



* COMMAND INTERPRETER 


0012 



* 


• 


0013 

0104 

08 00 

CMM D 

ENT 



0014 

0105 

OAOO 


EIN 



0015 

0106 

4005 


Cl E 


ENABLE PANIC BUTT3N 

CD 1 6 

0107 

F9 00 
0000 

. 

J ST 

CRLF 

PRINT PRSMPT CHARACTER 

0017 

0103 

C687 


LAP 

; 8 7 


00 IS 

0109 

F9 00 
0000 


j ST 

0TT 


0019 

010A 

C6AA 


LAP 

•* • 


0020 

0103 

F9 00 
0000 


J ST 

0TT 


0021 

010C 

0108 


ZXR 



0022 

010D 

F9 00 
0000 


J ST 

IER 

INPUT COMMAND 

0023 

010E 

C0D4 


CAI 

t ■J* » 


0024 

010F 

E213 


LIK 

TRANS 


0025 

0110 

C0D2 


CAI 

»R, 


0026 

0111 

E212 


LDX 

READ 


0027 

0112 

C0C3 


CAI 

•c 


0028 

0113 

E21 1 


LDX 

CLEAR 


0029 

0114 

38 OA 


JXN 

SKF 


0030 

0115 

CODO 


CAI 

*P* 


0031 

0116 

E20F 


LEOC 

PUT CH 


0032 

0117 

C0C7 


CAI 

• 6 ' 


0033 

01 13 

E20F 


L IX 

FPC0T 


0034 

0119 

C0D6 


CAI • 

•V* 


0035 

01 1A 

E20E 


LDX 

CRT 


0036 

011B 

C0C4 


CAI 

* D f 


0037 

01 1C 

E20A 


LDX 

Datum 


0033 

01 ID 

33 01 


JXN 

SKF 


0039 

01 IE 

F9 00 

BAD 

J ST 

ERR 




0000 



• 


0040 

0 1 IF 

EAOA 

0KF 

STX 

FUMC 

SAVE CALL ADDRESS 

0041 

0120 

F9 00 
0000 


J ST 

XEQ 

WAIT F0R G0 

0042 

0121 

FB03 


J ST 

* FUN C 

CALL FUNCTI0M 

0043 

0122 

F61D 


JMP 

CMN EH- 1 


0044 

0123 


TRAMS 

REF 



0045 

0124 


READ 

REF 



0046 

0125 


CLEAR 

REF 




0047 

0126 


PUM CH 

0043 

0127 


DATUM 

0049 

0123 


FPL0T 

0050 

0129 


CRT 

0051 

012A 

0000 

FUI\1C 

0052 

0 12B 

B1B5 

NAME 


0 1 2C 

AAB2 



0 12D 

B5B5 



012E 

AOC 6 



012F 

D5CC 



0130 

CCAO 



0131 

C3D4 



0132 

D3A0 


0053 

0133 

SD3A 



0134 

0000 


0054 

0000 

ERR0RS 




REF 

REF 

REF 

REF 

DATA 0 

TEXT '154=255 FULL. HITS' 


DATA : 8 D3 A» 0 

EM D 


0001 


0002 


0003 


OOOA 

0000 

0005 

0000 

0006 

0001 

0007 

0002 

0008 

0003 

0009 

OOOA 

0010 

0005 

001 1 

0006 

0012 

0007 

0013 

0003 

001 A 

0009 

0015 

OOOA 

0016 

OOOB 

0017 

OOOC 

00 io 

GOOD 

0019 

OOOE 

0020 

OOOF 

0021 

0010 

0022 

0011 

0023 

0012 

002A 

0013 

0025 

00 1 A 

0026 

0015 

0027 

0016 

0023 

0017 

0029 

0013 

0030 

0019 

0031 

00 1 A 

0032 

00 IB 

0033 

001C 

00 3 A 

00 1 D 

0035 

00 1 E 

0Q36 

00 IF 

0037 

0020 

0038 

0021 

0039 

0022 


03 00 TRANS 
C6BF 
F9 00 
8 127 
013E 
F9 00 
8 128 
1200 
COCD 
F203 
CODA 
F20A 
F9 00 
8 12A 

F900 MSN 
8 129 
0103 
F20F 

F9 00 TAP 
8 129 
C63 A 
F9 00 
8 126 
F9 00 
8 125 
2 1 A 1 

F9 00 T 1 
8 126 
F9 00 
8 125 
C09 2 
F20i 
F60A 

F9 00 T2 

8 125 
COFF 
F20 1 
F603 

CA01 TA 
E9 00 BEG 
012B 
B 1 00 • 

0122 

9 AF3 
F9 00 
5 129 
C68F 
F9 00 
8 127 


NAM 

TRANS 

NAM 

RTTR 

DCTR 

DI SP 

REL 

ENT 

0 

LAP 

i •> * 

J ST 

*0TL 

DATA 

M 0 T A 

J ST 

*IKB 

R0V 

CAI 

*M ’ 

JMP 

MSN 

CAI 

t f I 

JMP 

tap 

J ST 

*ERR 

J ST 

*CRLF 

ZXR 

JMP 

BEG 

J ST 

=S= CRL F 

LAP 

: 8 A 

J ST 

*0TT 

J ST 

*BIPT 

JAZ 

S- 1 

J ST 

*0TT 

J ST 

*BI PT 

CAI 

: 9 2 

JMP 

T2 

JMP 

T I 

J ST 

*EIPT 

CAI 

: FF 

JMP 

ta 

JMP 

T2 

LXP 

1 

STX 

TETR 

L DA 

MASKO 

STA 

MASKi 

J ST 

* CRL F 

LAP 

’? ' 

J ST 

’ *0TL 


171 


0040 

0023 

0041 

0024 

0042 

0025 

0043 

0026 

0044 

0027 

0045 

0023 

0046 

0029 

0047 

002A 

0043 

002B 

0049 

002C 

0050 

002D 

0051 

002E 

0052 

002F 

0053 

0030 

0054 

0031 

0055 

0032 

0056 

0033 

0057 

0034 

0053 

0035 

0059 

0036 

0060 

0037 

0061 

0033 

0062 

0039 

0063 

003A 

0064 

003B 

0065 

003C 

0066 

. 003D 

0067 

003E 

0063 

003F 

0069 

0040 

007 0 

0041 

007 1 

0042 

0072 

0043 

007 3 

0044 

0074 

0045 

0075 

0046 

0076 

0047 

0077 

0043 

0073 

0 049 

0079 

004 A 

003 0 

004B 

J03 1 

004C 

003 2 

004D 

003 3 

004E 

003 4 

004F 

003 5 

0050 

00# 6 

0051 

003 7 

0052 

OOSS 

0053 

0039 

0054 

009 0 

0055 

009 1 

0056 


0130 
F9 00 
6 123 
C0C6 
F205 
C0C2 
F203 
COCF 
F20A 
F3FE 

C71E F0RD 

9AE2 

0103 

0403 

F207 

C7 1E BACK 
9 ADD 
0103 
F203 

C70F BOTH 
9 AD9 
0353 

EAF5 TKF 
E2F1 
23 02 
C70F 
9 AD 3 
E2D4 
9 AD2 
E2S3 
0110 

9 COO CLR 
0123 
DACD 
F603 
EAD1 
E2E4 
23 0 3 
E2E5 
33 01 
0210 

E2D3 N EXT 
0123 
0210 
CGOO 
C2FF 
EA3S 
9 ABA 
C7F5 
9 ACC 
C7FF 
9 A37 
B2D4 


DATA 

DI RE 

J ST 

*IKB 

CAI 

«F> 

JMP 

FORD 

CAI 

' B ’ 

JMP 

BACK 

CAI 

•0 * 

JMP 

B0TH 

J ST 

* ERR 

LAM 

30 

STA 

ZXR 

CNT2 

CXR 


JMP 

TKF 

LAM 

30 

STA 

CMT2 

ZXR 


JMP 

TKF 

LAM 

15 

STA 

CNT2 

AX P 


STX 

DI Rl 

LDX 

TETR 

JXZ 

£+3 

LAM 

15 

STA 

CMT2 

L DA 

ZCT 

STA 

CMT3 

LDX 

I3C 

Z.AR 


STA 

60 

IX R 


IMS 

CNT3 

JMP 

CLR 

STX 

IBP 

LDX 

TETR 

JXZ 

NEXT 

LDX 

DI Rl 

JXN 

NEXT 

CAR 


LDX 

MASK 

IX R 


CAR 


CAI 

0 

AX I 

255 

STX 

MSK1 

STA 

DIR' 

LAM 

245 

STA 

DCT2 

L All 

.255 

STA 

DCT 

L DA 

TETR 


si 


> 


172 


009 2 

0057 

2107 


JAZ 

TUR 1 

009 3 

0058 

FBCC 

T3 

J ST 

*BI PT 

009 4 

0059 

1357 


LLA 

8 

009 5 

005A 

9 AD4 


STA 

TEMT 

009 6 

005B 

FBC9 


J ST 

*BI PT 

009 7 

005C 

A2D2 


IBR 

TEMT 

0098 

005D 

9 ABO 


STA 

VAL 

0099 

005E 

F209 


JMP 

PRBC 

0100 

005F 

FAS 9 

TUR1 

J ST 

TURN0N 

0101 

0060 

FA9 3 

1LP 

J ST 

MBN S 

0102 

0061 

9 AAC 


STA 

VAL 

0103 

0062 

32CB 


LDA 

DIR1 

0104 

0063 

COOl 


CAI 

r 

0105 

0064 

F203 


JMP 

PRBC 

0106 

0065 

D2A6 


CMS 

DIR 

0107 

0066 

F235 


JMP 

FIN0NE 

0108 

0067 

F234 


JMP 

FIM0NE 

0109 

0068 

C7 OF 

PH0C 

LAM 

15 

0110 

0069 

9AA5 


STA 

CMT1 

01 1 1 

006A 

B2A3 


LDA 

MASK 1 

0112 

006B 

9 AA8 


STA 

MASK 2 

0113 

006C 

B2B6 


LDA 

IBC 

01 14 

006D 

9 AAS 


STA 

13 

01 15 

006E 

C7FF 

SPA 

LAM 

255 

0116 

006F 

9AA5 


STA 

SCT 

01 17 

0070 

3298 


LDA 

M SK 1 

01 18 

007 1 

9 A98 


STA 

M SK 2 

01 19 

0072 

B2A3 


LDA 

IB 

0120 

0073 

9 AA3 


STA 

IBP 

0121 

0074 

01 10 


ZAR 


0122 

"075 

D39 E 


CMS 

*MASK2 

0123 

0076 

F201 


JMP 

S+2 

0124 

0077 

F2 1 3 


JMP 

TLP2 

0125 

0078 

E29 5 

TL.P1 

L DM 

VAL 

0126 

0079 

0110 


ZAR 


0127 

007A 

D3SF 


CMS 

*MSK2 

0128 

007B 

F201 


JMP 

S+2 

0129 

007 C 

0508 


MXR 


0130 

007D 

0030 


TXA 


0131 

007E 

8 B9S 


ADD 

* I BP 

0132 

007F 

9B9 7 


STA 

*IBP 

0133 

008 0 

DA39 


IMS 

M$K2 

0134 

008 1 

DA9 5 


IMS 

IBP 

0135 

008 2 

DA9 2 


IMS 

SCT 

0136 

003 3 

F60B 


JMP 

TLP1 

0137 

003 4 

DAS F 


IMS 

MASK 2 

0138 

003 5 

C6FF 


LAP 

255 

0139 

003 6 

8 AS F 


ADD 

IB 

0140 

003 7 

9A8E 


STA 

IB 

0141 

0038 

DAS 6 


IMS 

CMT 1 

0142 

0089 

F61B 


JMP 

SPA 

C143 

008 A 

F21 1 


JMP 

* FIN0NE 

0144 

008B 

E28 2 

TLP2 

LEK 

VAL 


0145 

008 C 

0146 

008 D 

0147 

008 E 

0143 

003 F 

0149 

009 0 

0150 

009 1 

0151 

009 2 

0152 

009 3 

0153 

009 4 

0154 

009 5 

0155 

009 6 

0156 

009 7 

0157 

0098 

0158 

0099 

0159 

009 A 

0160 

009 B 

0161 

009 C 

0162 

009 D 

0163 

009 E 

0i64 

009 F 

0165 

00 AO 

0166 

OCA 1 

0167 

00A2 

0168 

OOA3 

0169 

00A4 

0170 

00A5 

017 1 

00A6 

0172 

00A7 

0173 

00 AS 

0174 

00A9 

0175 

OCAA 

0176 

OOAB 

0177 

00 AC 

0178 

OOAD 

0179 

OGAE 

013 0 

OOAF 

013 1 

0050 

018 2 

0CB1 

013 3 

0032 

018 4 

00B3 

0185 

0034 

0186 

00E5 

0187 

00E6 

0188 

0CB7 

0189 

00B8 

0190 

00B9 

019 1 

OOBA 

019 2 

OOEB 

01*9 3 

OOBC 

0194 

OOED 

0195 

OOBE 

0196 

OOBF 

0197 

ooco 


Oi 10 
D37C 
0508 
0030 
8B86 
9 38 5 
DA7 7 
DA3 3 
DAS 0 
r oOA 
DA7D 
C6FF 
8A7D 
9A7C 
DA74 
F62D 
E26C FIN8NE 
0128 
326D 
COOO 
C302 
EA67 
DA7D 
F20 1 
40C7 
DA67 
F20 1 
F204 

B23 2 TUR2 
2101 
F652 
F64B 

C7FF TUR3 
9 A6E 
E263 
0110 

9 COO CLR 1 

0123 

DA 69 

F603 

C7F 1 

9A63 

C7 OF 

9 A62 

9A62 

B269 

9 A63 

C70E 

9A60 

B260 MX BIN 
9 A60 

B35F NXR0W 
1 0D3 


ZAR 


CMS 

*MSK2 

MXR 


TX A 


ADD 

*IBP 

STA 

*IBP 

IMS 

M SK 2 

IMS 

IBP 

IMS 

SCT 

JMP 

TLP2 

IMS 

MASK 2 

LAP 

255 

ADD 

IB 

STA 

IB 

IMS 

CMT1 

JMP 

SPA' 

LEX 

M SK 1 

IX R 


LDA 

DIR 

CAI 

0 

SXI 

2 

SIX 

M SK 1 

IMS 

DCT2 

JMP 

S+ 2 

SEL 

24j 7 

IMS 

DCT 

JMP 

TUR2 

JMP 

TUR3 

LDA 

TETR 

JAZ 

S+2 

JMP 

T 3 

JMP 

ILP 

LAM 

255 

STA 

CM T 7 

LEX 

IBP 

ZAR 


STA 

©0 

IX R 


IMS 

CNT7 

JMP 

CL R 1 

LAM 

241 ' 

STA 

CMT4 

LAM 

15 

STA 

CMT5 

STA 

CMT6 

LDA 

IBC 

STA 

IBD 

LAM 

14 

STA 

HM 14 

LDA 

IBD 

STA 

I BE 

LDA 

’*IBE 

ARA 

4 


n[ 


! ... i I 


0198 

00C1 

8B55 


ADD 

I BP 

0199 

00C2 

9B54 


ST A 

*IBP 

0200 

00C3 

B254 


L DA 

H256 

0201 

00C4 

8 A5A 


ADD 

IBE 

0202 

00C5 

9A59 


STA 

I BE 

0203 

00C6 

DA53 


IMS 

CNT5 

0204 

00C7 

F603 


JMP 

NXR0U 

0205 

00C8 

B252 


LDA 

CNT6 

0206 

00C9 

9 A50 


STA 

CNT5 

0207 

OOCA 

DA53 


IMS 

IBL 

02 OB 

OOCB 

DA4B 


IMS 

IBP 

0209 

oocc 

DA4C 


IMS 

CNT4 

0210 

OOCD 

F610 


JMP 

NX BIN 

021 1 

OOCE 

B24B 


LDA 

CMT5 

0212 

OOCF 

0150 


IAR 


0213 

OODO 

9 A4A 


STA 

CNT 6 

0214 

00D1 

9A48 


STA 

CNT5 

0215 

00D2 

C701 


LAM 

1 

0216 

OOD3 

9A45 


STA 

CNT 4 

0217 

00D4 

DA48 


IMS 

HM 14 

0218 

OOD5 

F61S 


JMP 

NX BIN 

0219 

00D6 

0110 


2 AH 


0220 

00D7 

9A49 


STA 

OFFSET 

0221 

00D3 

F100 


JMP 

DI SP 



0000 




022 2 

00D9 

B251 

RTTR 

LDA 

TETR 

0223 

CODA 

2104 


JA2 

T5 

0224 

OODB 

B252 


LDA 

DIR1 

0225 

OODC 

313 4 


JAG 

T 6 ' 

0226 

OODD 

0210 


CAR 


0227 

OODE 

F206 


JMP 

T7 

0228 

OODF 

B24E 

T5 

LDA 

DIR1 

0229 

OOEO 

218 2 


JAL 

£+ 3 

0230 

OOE1 

322 A 

T6 

LDA 

DIR 

0231 

00E2 

F202 


JMP 

T7 

0232 

00E3 

B22S 


LDA 

DIR 

0233 

00E4 

3101 


JAM 

S+2 

0234 

00E5 

DA2D 

T7 

IMS 

MASK 1 

0235 

00E6 

DA29 


IMS 

CNT2 

0236 

00E7 

F69C 


JMP 

NEXT 

0237 

DOES 

F7E3 


RTN 

TRANS 

02 3B 

00E9 

03 0 0 

TURN0N 

ENT 


0239 

OOEA 

4006 


Cl D 


0240 

OOEB 

48 Cl 


SSM 

: Cl 

0241 

OOEC 

F60 1 


JMP 

5- 1 

0242 

OOED 

OEOO 


SBM 

• * 

0243 

OOEE 

49 Cl 


SEN 

: Cl 

0244 

OOEF 

F601 


JMP 

1 

0245 

OOFO 

OFOO 


SUM 

* 

0246 

OOF 1 

40C4 


SEL 

: C4 

0247 

00F2 

4005 


Cl E 


0248 

00F3 

F70A 


RTN 

• TURN0N 

0249 

00F4 

03 00 

MON 5 

EnJT 



175 


0250 

OOFF 

49 C6 


SEN 

: C6 : 

0251 

00F6 

F601 


JMP 

S- 1 • 

0252 

00F7 

SAC 6 


I NX 

: C6 

0253 

00F3 

28 43 


JXZ 

S- 3 

0254 

00F9 

C704 


LAM 

4 ‘ 

0255 

00 FA 

9 A 1 0 


STA 

CNT 

0256 

OOFB 

0110 


EAR 


0257 

OOFC 

9 AOD 


STA 

M SK 2 

0253 

OOFD 

F206 


JMP 

M2 

0259 

OOFE 

1350 

M 1 

LLA 

1 

0260 

OOFF 

9 AO A 


STA 

M SK 2 

0261 

0100 

1351 


LLA 

2 

0262 

0101 

SA03 


ADD 

M SK 2 

0263 

0102 

9 A07 


STA 

MSK2 

0264 

0103 

0110 


EAR 


0265 

0104 

1B03 

M2 

LLL 

4 

0266 

0105 

8 A04 


ADD 

M SK 2 

0267 

0106 

DA04 


IMS 

CNT 

02 63 

0107 

F609 


JMP 

M 1 

0269 

0103 

F7 1 4 


RTN 

MSN S 

027 0 

0109 

0000 

M SK 1 

DATA 

0 

027 1 

010A 

0000 

MSK2 

DATA 

0 

027 2 

01 OB 

0000 

CNT 

DATA 

0 

027 3 

OlOC 

0000 

DIR 

DATA 

0 

0274 

010D 

0000 

DOT 

DATA 

0 

027 5 

010E 

0000 

VAL 

DATA 

0 

0276 

01 OF 

0000 

CM T 1 

DATA 

0 

0277 

0110 

0000 

CM T 2 

DATA 

0 

0278 

Oil 1 

0000 

CM 7 3 

data 

0 

0279 

0112 

F010 

Z CT 

DATA 

-403 0 

028 0 

0113 

0000 

MASK 1 

DATA 

0 

028 1 

0114 

0000 

MASK 2 

DATA 

0 

028 2 

0115 

0000 

SCT 

DATA 

0 

023 3 

0116 

0000 

IB 

DATA 

0 

023 4 

0117 

0000 

IBP 

DATA 

0 

023 5 

01 IS 

0100 

ri 256 

DATA 

256 

023 6 

0119 

0000 

CNT 4 

DATA 

0 

028 7 

011A 

0000 

CNT 5 

DATA 

0 

0288 

01 IB 

0000 

CNT 6 

DATA 

0 

0289 

01 1C 

0000 

CNT7 

DATA 

0 

029 0 

01 ID 

FFF2 

HM 1 4 

DATA 

- 14 

029 1 

01 IE 

0000 

IBD 

DATA 

0 

029 2 

01 IF 

0000 

I BE 

DATA 

0 

029 3 

0120 

0000 

DCT2 

DATA 

0 

029 4 

0121 

0000 

OFFSET 

DATA 

0 

029 5 

0122 


MASKO 

REF 


029 6 

0123 


IBC 

REF 


029 7 

0124 


MASK 

REF 


0298 

0125 


BIPT 

REF 


0299 

0126 


0TT 

REF 


0300 

0127 


0TL 

REF 


0301 

0128 


1KB 

REF 

• 

0302 

0129 


CRLF 

REF 



4 


176 


0303 

012A- 


ERR 

REF 

0304 

0 1 2B 

0000 

TETR 

DATA 

0305 

012C 

0000 

TEDA 

DATA 

0306 

0 12D 

OFFO 

TCT1 

DATA 

0307 

012E 

0000 

DIRl 

DATA 

03 OB 

012F 

0000 

TE'lT 

data 

0309 

0130 

8D3A 

DIRE 

DATA 

0310 

0131 

C6CF 


TEXT 


0132 

D2D7 




0133 

C1D2 




0134 

C4AC 




0135 

C2C1 




0136 

C3CS 




0137 

D7C 1 




0133 

D2C4 




0139 

ACCF 




013A 

D2A0 




013B 

C2CF 




013C 

D4CS 




0 1 3D 

BFAO 



0311 

013E 

8D3A 

M0TA 

DATA 

0312 

013F 

C6D2 


TEXT 


0140 

CFCD 




0141 

AOGD 




0142 

CFCE 




0143 

D3C1 




0144 

CED4 




0145 

CFAO 




0146 

CFD2 




0147 

A0C6 




0148 

D2CF 




0149 

CDAO 




014A 

D4C1 




014B 

D0C5 




0 14C 

BFAO 



0313 




END 

0000 

ERRORS 





i 


o 

0 

403 0 
0 
0 

:SD3A 

1 FORWARD, BACKWARD, 0R BOTH? * 


: S D3 A 

* FROM MONSANTO OR FROM TAPE?' 



0001 




NAM 

DATUM/ RTDA 

0002 




EX TR 

01 SD 

0003 

0000 

03 00 

DATUM 

EMT 


0004 

0001 

C6BF 


LAP 

■ 

0005 

0002 

FB63 


J ST 

*0TL 

0006 

0003 

0067 


DATA 

DIRE 

0007 

0004 

FB7 1 


J ST 

*IKB 

0008 

0005 

C0C6 


CAI 

•F* 

0009 

0006 

F205 


JMP 

F0RD 

0010 

0007 

C0C2 


CAI 

*B» 

00 1 1 

0008 

F208 


JMP 

BACK 

0012 

0009 

COCF 


CAI 

.0 • 

0013 

000A 

F20A 


JMP 

B0TH 

0014 

0003 

FB6C 


J ST 

* ERR 

0015 

OOOC 

07 IE 

F0RD 

LAM 

30 

0016 

OOOD 

9 A7 3 


STA 

ENT2 

0017 

OOOE 

0108 


ZXR 


0018 

OOOF 

0408 


CXR 


0019 

0010 

F207 


JMP 

PKF 

0020 

0011 

C71E 

BACK 

LAM 

30 

0021 

0012 

9 A6E 


STA 

EMT2 

0022 

0013 

0103 


ZXR 


0023 

0014 

F203 


JMP 

PKF 

0024 

0015 

C70F 

B0TH 

LAM 

15 

0025 

0016 

9A6A 


STA 

DMT 2 

0026 

0017 

0353 


AXP 


0027 

0018 

EA60 

FKF 

STX 

DIRl 

0023 

0019 

B260 


L DA 

ZCT1 

0029 

00 1 A 

9A60 


STA 

CT 1 

0030 

00 IE 

E25B 


LEX 

IBC 

0C3 1 

001C 

EA5F 


STX 

IB 

0032 

00 1 D 

0110 


ZAR 


0033 

0C1E 

9 COO 

CLR 

STA 

to 

0034 

00 1 r 

0128 


IX R 


0035 

0020 

DA5A 


IMS 

CT 1 

0036 

0021 

F 603 


JMP 

CLR 

0037 

0022 

9A61 


STA 

R0N 

0038 

0023 

0210 

NEXT 

CAR 


0039 

0024 

9 A53 


STA 

. DIRO 

0040 

0025 

C7F5 


LAM 

245 

0041 

0026 

9 A53 


STA 

DDT 2 

0042 

0027 

C7FF 


LA-1 

255 

004 3 

0028 

9A55 


STA 

DDT 

0044 

0029 

FA19 


J ST 

DURN0N 

0045 

002A 

FA23 

CLP 

J ST 

DSN S 

0046 

002B 

9A54 


STA 

DAT 

0047 

002C 

B24C 


L DA 

DIRl 

0048 

002D 

C001 


CAI 

1 

0049 

002E 

F203 


JMP 

ST0 

0050 

002F 

D24D 


CMS 

DIRO 

0051 

0030 

F232 


JMP 

DLP2 

0052 

0031 

F23 1 


JMP 

• DLP2 

0053 

0032 

B24D 

ST0 

L DA 

DAT 


i 


0054 

0033 

£248 


LDX 

IB - 

0055 

0034 

9 COO 


STA 

eo ’ 

0056 

0035 

DA46 


IMS 

IB 

0057 

0036 

DA4S 


IMS 

DDT 2 

0053 

0037 

F20 1 


JMP 

5+2 

0059 

0038 

40C7 


SHL 

24; 7 

0060 

0039 

DA44 


IMS 

DDT 

0061 

003A 

F610 


JMP 

DLP 

0062 

003B 

DA48 


IMS 

R0W 

006 3 

003C 

B247 

DLP 1 

L DA 

R0M 

0064 

003D 

2101 


JAZ 

S+2 

0065 

003E 

FlOO 

0000 


JMP 

DISD 

0066 

003F 

B23D 

RTDA 

LDA 

DIRO 

0067 

0040 

DA40 


IMS 

DWT2 

00 6B 

0 041 

F61E 


JMP 

NEXT 

0069 

0042 

F742 


RTM 

DATUM 

007 0 

0043 

03 00 

DUR'J 0N 

EMT 


0071 

0044 

4006 


Cl D 


0072 

0045 

43 Cl 


SSW 

: Cl 

0073 

0046 

F601 


JMP 

S- 1 

0074 

0047 

OEOO 


SBM 


0075 

0043 

49 Cl 


SEW 

: Cl 

0076 

0049 

F601 


JMF 

5- 1 

0077 

004A 

OFOO 


SUM 

* 

0078 

0 04B 

40C4 


SEL 

: C4 

0079 

004C 

4005 


Cl E 


003 0 

004D 

F70A 


RTW 

DURN0N 

003 1 

004E 

03 00 

D0WS 

EMT 


003 2 

004F 

49 C6 


SEW 

: C6 

003 3 

0050 

F601 


JMP 

S- 1 

003 4 

0051 

5AC6 


I MX 

: C6 

003 5 

0052 

26 43 


JXZ 

$-3 

003 6 

0053 

C704 


LAM 

4 ' 

003 7 

0054 

9 A2D 


STA 

CNT 

0088 

0055 

0110 


EAR 


0089 

0056 

9A2C 


STA 

MSK2 

009 0 

0057 

F206 


JMP 

M2 

009 1 

0053 

1350 

M 1 

LLA 

1 

009 2 

0059 

9A29 


STA 

MSK2 

009 3 

005A 

1351 


LLA 

2 

009 4 

35B 

8 A27 


ADD 

MSK2 

009 5 

005C 

9 A26 


STA 

MSK2 

009 6 

005D 

0110 


EAR 


0097 

005E 

1B03 

M2 

LLL 

4 

0098 

005F 

8 A23 


ADD 

MSK2 

0099 

0060 

DA21 


IMS 

CNT 

0100 

0061 

F609 


JMP 

Ml 

01*01 

0062 

F714 


RTN 

D0NS 

0102 

0063 

DAI A 

DLP2 

IMS 

DDT 

0103 

0064 

F63A 


JMP 

DLP 

0104 

0065 

F629 


JMP 

DLP 1 

0105 

0066 


0TL 

REF 



0106 

0067 

8 IB A 

DI RE 

DATA 

0107 

0068 

C6CF 


TEXT 


0069 

D2D? 




006A 

C1D2 




006B 

C4AC 




006C 

C2C1 




• 006D 

G3CB 




006E 

D7C1 




006F 

D2C4 




0070 

AOCF 




0071 

D2A0 




0072 

C2CF 




0073 

D4CS 




0074 

BFAO 



0108 

0075 


CBLF 

REF 

0109 

0076 


IKB 

REF 

01 10 

0077 


IBC 

REF 

0111 

0078 


ERR 

REF 

01 12 

0079 

0000 

DIR1 

DATA 

0113 

007A 

F003 

ZCT1 

DATA 

01 14 

0073 

0000 

GT1 

DATA 

0115 

007C 

0000 

IB 

DATA 

0116 

007D 

0000 

DIRO 

DATA 

01 17 

007E 

0000 

DDT 

DATA 

01 IS 

007F 

0000 

DDT 2 

DATA 

01 19 

003 0 

0000 

DAT 

DATA 

0120 

008 1 

0000 

DMT 2 

DATA 

0121 

008 2 

0000 

CMT 

DATA 

0122 

003 3 

0000 

M SK2 

DATA 

0123 

008 4 

0000 

ROM 

DATA 

0124 




END 

0000 

ERRORS 





:8D3A 

' FORWARD, BACKWARD 0R B0TH 


0 

-408 5 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0001 




MAM 

IKBi IER»BIPTjXEQ 


0002 




HAM 

0TT.> £RR> CRLF 


0003 




NAM 

QTLa 0FPA.0DEC 


0004 

0000 



REL 

0 



0005 



* 





0006 



* PAM I C SWITCH 



0007 



* 





0003 

0000 

FBOO 


J ST 

* CMND 



0009 

0001 


CMND 

REF 




0010 



* 





0011 



*KEY BOARD INPUT 



0012 



* 





0013 

0002 

08 00 

I KB 

ENT 




0014 

0003 

4033 


SEL 

7* 0 

AUT0- ECH0 


0015 

0004 

4039 


SEL 

7* 1 

KBD M0DE 


0016 

0005 

59 39 


RDA 

7» 1 

READ 0N FLAG 


0017 

0006 

403C 


SEL 

7* 4 

RESET 


0018 

0007 

F705 


RTN 

IKB 



0019 

0003 

08 00 

IER 

ENT 




0020 

0009 

FE07 


J ST 

IKB 



0021 

000A 

CODF 


CAI 

: DF 

BACK ARR0W? 


0022 

OOOB 

FFOA 


J-ST 

*CMND 

YES 


0023 

OOOC 

COS A 


CAI 

S 8 A 

LINE FEED? 


0024 

GOOD 

FFOC 


J ST 

* CMND 

YES 


0025 

OOOE 

F706 


RTN 

IER 

N El THER> RETURN 


0026 



* 





0027 



* PAPER 

TAPE 

INPUT 



0028 



* DSO = 

= 0 F0R TTY* 



0029 



* 1 FOR HSR 




0030 








0031 

OOOF 

03 00 

BI PT 

ENT 




0032 

0010 

53 01 

BIP2 

ISA 


READ SWITCHES 


0033 

0011 

13D0 


LRA 

1 

DSO UP F0R TTY 


0034 

0012 

220A 


JOS 

HSR 

DOWN F0R HSR 


0035 

0013 

49 3B 


SEN 

7» 3 

TTY BUSY? 


0036 

0014 

F604 


JMP 

BI P2 

YES 


0037 

0015 

403A 


SEL 

7j 2 

N0, STEP READER 


0038 

0016 

48 39 

WT 

SSN 

1* 1 

FLAG? 


0039 

0017 

F203 


JMP 

IT 

YES 


0040 

0018 

0150 


IAR 


NO, BUMP COUNT 


0041 

0019 

2149 


jAZ 

BIP2 

RESTART IF TIME 

UP 

0042 

001 A 

F604 


JMP 

WT 

ELSE CHECK FLAG 

AGAIN 

0043 

00 IB 

53 33 

IT 

IN A 

7a 0 

INPUT FROM TTY 


0044 

001C 

F70D 


RTN 

BI PT 

a RETURN 


0045 

00 1 D 

49 33 

HSR 

SEN 

6; 3 

HSR BUSY? 


0046 

00 1 E 

F60E 


JMP 

BIP2 

yes 


0047 

00 IF 

4032 


SEL 

6j 2 

N0j STEP READER 


0048 

0020 

48 35 

m 

SSN 

6j 5 

FLAG? 


0049 

0021 

F203 


JMP 

IH 

yes 


0050 

0022 

0150 


IAR 


N0j BUMP COUNT 

, 

0051 

0023 

2153 


JAZ 

BIP2 

RESTART IF TIME 

UP 

0052 

0024 

F604 


JMP 

• WH 

ELSE CHECK FLAG 

AGAIN 

0053 

0025 

53 35 

1H 

INA 

6a 5 

INPUT FR0M HSR 



0054 

0055 

0056 

0057 
0053 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
0067 
0063 
0069 

007 0 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 
00S0 

008 1 
00S2 
008 3 
OOS 4 
008 5 
003 6 

008 7 

0038 

0039 

009 0 
009 1 
009 2 
009 3 
009 4 
009 5 
009 6 
0097 
0093 

0099 

0100 
0101 

0102 

0103 

0104 

0105 

0106 


0026 

F717 


RTN 

BIPt 

& RETURN 



*UAIT 

FOR EXECUTE SIGNAL 

0027 

08 00 

He 

XEQ 

ENT 



0028 

FE20 


J ST 

IER 

INPUT 

0029 

COS D 


CAI 

: 8 D 

CARRIAGE RETURN? 

002A 

F703 


RTN 

XEQ 

YES, RETURN 

002B 

F603 


JMP 

X CiQ+ 1 

N0, GET MORE 



n* 

* OUTPUT T0 

TTY 


002C 

08 00 

0TT 

E‘JT 



002D 

403C 


SEL 

7, 4 

RESET INTERFACE 

002E 

6D3B 


WRA 

If 3 

WRITE m NOT BUSY 

002F 

49 3B 


SEN 

3 

DONE? 

0030 

F601 


JMP 

$“ 1 

NS 

0031 

F705 


RTN 

0TT 

YES 


* COMMAND ERROR EXIT 




0032 

08 00 

ERR 

ent 



0033 

C6DF 


LAP 

: DF 

PRINT ARROW 

0034 

FE03 


J ST 

0TT 


0035 

FF34 


JST 

He Ci'IN D 

RESTART COMMAND 



He CARRI AGE- RETURN, LINE 

FEED 

0036 

08 00 

CRLF 

ENT 



0037 

C63D 


LAP 

: 8 D 

CR 

0033 

FEOC 


JST 

0TT 


0039 

C6SA 


LAP 

: 8 A 

LF 

003A 

FEOE 


JST 

OTT 


003B 

F705 


RTN 

CRLF 




*T* 

♦OUTPUT TEXT 

Jj 

FROM BUFFER 

003C 

03 00 

0TL 

ENT 



003D 

8 AOF 


ADD 

CAI 

MAKE COMPARE INSTRUCTION 

003E 

9A06 


STA 

0T2 

&SAVE IT 

003F 

9A09 


STA 

0T3 


0040 

E704 


LDX 

HeSTL 

GET TEXT POINTER 

0041 

DE05 


IMS 

3TL 

SET RETURN ADDRESS 

0042 

B400 

9T1 

L DA 

00 

GET W0RD 

0043 

1 1D7 


RRA 

8 

PRINT FIRST BYTE 

0044 

FE13 


JST 

0TT 


0045 

COOO 

0T2 

CAI 

0 

LAST ONE • 

0046 

F70A 


RTN 

0TL 

YES, RETURN 

0047 

1 157 


RL A 

8 

PRINT SECOND BYTE 

0048 

FE1C 


JST 

0TT 


0049 

COOO 

0T3 

CAI 

0 

LAST ONE? 

004A 

F70E 


RTN 

' STL 

YES RETURN 

004B 

0128 


IXR 


BUMP POINTER 


0107 004C F60A 

0108 0C4D C000 

0109 

0110 
0111 

0112 0C4E 08 00 

0113 004F E701 

0114 0050 DE02 

0115 0051 EA01 

0116 0052 FBOE 

0117 0053 0000 

0113 0054 0059 

0119 0055 0110 

0120 0056 F£1 A 

0121 0057 0059 

0122 0058 F70A 

0123 0059 0000 

0124 0061 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 0062 0800 

0137 0063 EA19 

0138 0064 C4AB 

0139 0065 3032 

0140 0066 0310 

0141 0067 C202 

0142 006S 9 A15 

0143 0069 0030 

0144 006A FE3E 

0145 006B B215 

01 46 006C 9A13 

0147 006D C705 

0148 006E 9A10 

0149 006F B20E 01 

0150 0070 C4AF 

0151 0071 9 30E 

0152 0072 0128 

0153 0073 30C2 

0J.54 0074 S BOB ' 

0155 0075 9 AOS 

0156 0076 0030 

0157 0077 FE4B 

0153 0078 DA07 

0159 0079 DAO 5 



JMP 

0T1 • 

CAI 

* 

CAI 

0 ' 

*0UTPUT FL0ATING P0INT 

* 

0FPA 

ENT 



LBX 

*0FPA 


IMS 

0FPA 


STIC 

0PT 


JST 

* FAS 

0PT 

DATA 

0 


DATA 

ZAR 

BUF 


JST 

0TL 


DATA 

BUF 


RTN 

0 FPA 

BUF 

RES 

8, 0 

FAS 

REF 



SIX 

S 

LXP 

’+ * 

JAP 

£+3 

MAR 


AXI 

2 

STA 

V 

TXA 


JST 

0TT 

LDA 

STRT 

STA 

PTR 

LAM 

5 

STA 

T 

LDA 

V 

LXP 

: AF 

SUB 

*PTR 

IX R 


JAP 

S-2 

ADD 

*PTR 

STA 

V 

TXA 


JST 

0TT 

IMS 

' PTR 

IMS 

T 


LB0P 


NUMBER 


GET P0 INTER 
SET RETURN ADDRESS 
SAVE POINTER 
CBN VERT T0 ASCII 


SET END FLAG 
PRINT NUMBER 


SAVE XR 


MAKE VALUE + 

MAKE SIGN' - 
SAVE VALUE * 

PRINT SIGN 

INITIALIZE TEL PTR 

SET F0R 5 DIGITS 

ZERO T0 -1 


PRINT DIGIT 


* 

* 0DEC OUTPUT DECIMAL C+/-DDDDD) 

** 0DEC CONVERTS THE BINARY VALUE IN THE 
*A REG AND PRINTS I T AS A SIGNED 5 

* DIGIT DECIMAL NUMBER 0N THE TELETYPE* 

* AR AND OV ARE DESTROYED. 

* 

* LDA VAL AR * VALUE 

* SIM MUST BE IN W0RD MODE 

* JST *0 DEC CALL ROUTINE 

* *** RETURN XR UNCHANGED 
0DEC ENT 


0160 

007A 

F60B 


JMP 

0 1. 


0161 

007B 

E201 


LEX 

S 

RESTORE XR 

0166 

007C 

F71A 


RTN 

0DEC 

RETURN 

0163 

007D 

08 do 

S 

HLT 


TEMP F0R XR 

01 64 

007E 

08 00 

V 

HLT 


VALUE 

0165 

007F 

08 00 

T 

HLT 


C0UNT 

0166 

003 0 

0800 

PTR 

HLT 


POINTER 

0167 

003 1 

003 2 

STRT 

DATA 

TBL 

TABLE ADDR 

0163 

005 2 

2710 

TBL 

DATA 

loooo. 

1000, 100 


008 3 

03E8 






008 4 

0064 





0169 

008 5 

OOOA 


DATA 

10> 1 



008 6 

0001 





0170. 

0000 

ERRORS 



END 




0001 

0002 

0003 

0004 

0000 

0000 

03 00 

READ 

MAM 

EXTR 

REL 

EMT 

0005 

0001 

F9 00 


J ST 

0006 

0002 

0000 
F9 00 


J ST 

0007 

0003 

0000 

COSO 


CAI 

0008 

0004 

F203 


JMP 

0009 

0005 

C0C7 


CAI 

0010 

0006 

F204 


JMP 

oo n 

0007 

F9 00 


J ST 

0012 

0008 

0000 

E21F 

REZE 

LDK 

'0013 

0009 

B21F 


LDA 

0014 

000 A 

F204 


JMP 

0015 

00 OB 

B21C 

RESM 

LDA 

0016 

000C 

8 A1E 


ADD 

0017 

OOOD 

0043 


TAX 

oo is 

OOOE 

B21B 


LDA 

0019 

000F 

9 A1C 

RELP 

STA 

0020 

0010 

C6S A 


LAP 

0021 

001 I 

F9 00 


J ST 

0022 

0012 

0000 
F9 00 


J ST 

0023 

0013 

0000 

2141 


JAZ 

0024 

0014 

F9 00 

R 1 

J ST 

0025 

0015 

0000 
F9 00 


J ST 

0026 

0016 

0000 
C09 2 


CAI 

0027 

0017 

F20 1 


JMP 

0023 

0018 

F 604 


JMP 

0029 

0019 

F9 00 

R2 

J ST 

0030 

00 1 A 

0000 

COFF 


CAI 

0031 

00 1 B 

F20 1 


JMP 

0032 

001C 

F603 


JMP 

0033 

00 1 D 

F9 00 

R3 

J ST 

0034 

00 IE 

0000 

1357 


LLA 

0035 

00 1 F 

9 AOE 


STA 

C036 

0020 

F9 00 


J ST 

0037 

0021 

0000 

A20C 


I0R 

0033 

0022 

6 COO 


ADD 

0639 

0023 

9 COO 


STA 

0040 

0024 

0123 


IX R 

0041 

0025 

DAO 6 


IMS 

0042 

0026 

F609 


JMP 

0043 

0027 

F727 


RTW 


READ# PUNCH* CLEAR 
I KB# 0 TT# BI PT# 0TL# CRLF# I ER* ERR 
0 

CRLF 

IER 

* 0 * 

REZE 
'■6 * 

RESM 

ERR 

IBC 
Z CTO 
RELP 
IBC 
SMO 

ZCT1 

CMT 

:SA 

0TT 

BIPT 

S-l SKIP LEADER 

0TT 

BI PT READ TAPE 

I 9 2 E\1D 0 F LABEL 

R2 YES 

R1 N0 

BIPT READ TAPE 

IFF FILE MARK 

R3 YES 

R2 M0 

BIPT 

S 

TEMP 
BI PT 

TEMP 

§0 

@0 

i 

CMT 
• R‘3 
READ 


185 


0044 0028 

0045 0029 

0046 002A 

0047 0025 

0043 0020 

0049 002D 

0050 002E 

0051 

0052 

0053 002F 

0054 0030 

0055 0031 

0056 0032 

0057 0033 
C05B 0034 

0059 0035 

0060 0036 

0061 0037 

0062 0038 

0063 0039 

0064 003A 

0065 0035 

0066 0030 

0067 003D 

0068 003E 

0069 003F 

0070 0040 

0071 0041 

0072 0042 

0073 0043 

0074 0044 

0075 0045 

0076 0046 

0077 0047 

0078 0043 

0079 0049 

008 0 004A 

003 1 0 0 4 5 

008 2 0040 

008 3 

008 4 0 04D 

003 5 004E 

008 6 004F 

0037 0050 

0038 0051 

0089 0052 

009 0 0053 


1 50 

F010 Z CTO 
FF01 ZCT1 
OEF 1 SNO 
0000 CNT 
0000 PNT 
0000 TEMP 
* 

* PUNCH 
03 00 PUMCH 
F9 00 

0000 
F9 GO 
0000 
COSO 
F203 
CGC7 
F204 
F9 00 
0000 

E60F PNZE 

B60F 

F204 

B612 FMSN 

8 El 0 
0043 
5613 

9 El 2 PMLP 
FAOD 

F9 00 PI 
0000 
C09 2 
F201 
F603 

C6FF P2 
F9 00 
0000 

B400 P3 

FAOD 

0123 

DE1D 

F 604 

FAO 1 

F71D 

* PUNCH 
08 00 LEAD 
C732 

9E22 

0110 

F9 0 0 L 2 
0000 
DE25 
F602 


REF 

DATA -403 0 
DATA “ 255 

DATA 3325 
DATA 0 
DATA 0 
DATA 0 


ENT 

J ST CRL F 

J ST I ER 

CAI ’ 0 • 

JMP PNZE 

CAI ’ G ’ 

JMP PM SN 

JS1T ERR 

LDX IBC 

L DA Z CTO 

JMP PiML P 

L DA IBC 

ADD SN 0 

TAX 

L DA Z CT 1 

STA CNT 

J ST L EAD 

J ST IKB 

CAI : 9 2 

JMP P2 

JMF PI 

LAP : FF 

J ST 0TT 

L DA 60 

J ST 0TU 

IXR 

IMS CNT 

JMP P3 

J ST L EAD 

HTM PUN CH 

5" 0F LEADER 
EM T * 

LAM 50 

STA PNT 

Z AR 

J ST 0TT 

IMS 
JMP 


ECH0 LABEL 

CTRL /TAPE 

YES 

N0 

PUNCH FILE MARK 


PUMCH BITS 



PNT 
L 2 


! 


009 1 0054 

009 2 

009 3 0055 

0094 0056 

009 5 0057 

009 6 0 058 

009 7 0059 

0098 005A 

0099 

0100 005B 

0101 005C 

0102 0050 

0103 005E 

0104 005F 

0105 0060 

0106 0061 

0107 0062 

0103 0063 

0109 0064 

0110 0065 

0111 0066 

0112 0067 

0113 0068 

0114 0069 

0115 006A 

0116 0C6B 

0117 006C 

01 18 006D 

0119 006E 

0120 006F 

0121 0070 
0122 

0000 ERRORS 


F707 

03 00 0TW 
1 1D7 
F9 00 
0000 
1 157 
F9 00 
0000 
F705 

* CL EAR 
08 00 CLEAR 
F9 00 
0000 
F9 00 
0000 
COBO 
F203 , 

C0C7 
F204 • 

F9 00 
0000 

E63B CLZE 

B63B 

F204 

B63E CL SN 
SE3C 
0048 
B63F 

9E3E CLLP 
0110 

9C00 CR 
0128 
DE42 
F603 
F715 


RTN 

LEAD 

EM T 

, 

RRA 

a 

J ST 

0TT 

RL A 

8 

J ST 

0TT 

RTN 

0TW 

PROGRAM 

EN T 


J ST 

CRLF 

J ST 

IER 

CAI 

• 0‘ 

JMP 

CLZE 

CAI 

* G * 

JMP 

CLSN 

J ST 

ERR 

L DM 

IBC 

LDA 

Z CTO 

JMP 

CLLP 

LDA 

IBC 

ADD 

SNO 

TAX 


LDA 

ZCT1 

STA 

CMT 

ZAR 


STA 

£0 

IXR 


IMS 

CMT 

JMP 

CR 

RTN 

CLEAR 

EM D 



187 


ODO 1 



MAM 

FPL0-T 

0002 

0000 


REL 

0 

0003 

0000 

OS 0 0 

FPL.0T ENT 


0004 

0001 

C70F 

LAM 

15 

0005 

0002 

9 ABB 

STA 

NRS'W 

0006 

0003 

B2D0 

L DA 

IBC 

0007 

0004 

9 ABA 

STA 

ib ; 

00 OS 

0005 

0108 

ZXR 


0009 

0006 

FBCF 

J ST 

*CRLF 

0010 

0007 

FBD1. 

' J ST 

41 ER . 

0011 

0008 

COBO 

CAI 

•O’ 

0012 

0009 

F221 

JMP 

?ZR0 

0013 

OOOA 

GOBI 

CAI 

* 1 * 

0014 

OOQB 

F222 

JMP 

F0ME 

0015 

OOOC 

C0B2 

CAI 

•2* 

0016 

OOOD 

F222 

JMP 

FT 00 

.0017 

00.0 E 

COB 3 

CAI 

* 3 * 

00 IS 

00 OF 

F.222 

JMP 

FTHRE 

0019 

0010 

C0B4 

CAI 

’4* 

0020 

0011 

F222 

JMP 

FF0UR 

0021 

0012 

C0B5 

CAI 

1 5 * 

0022 

0013 

F222 

JMP 

FFfVE 

0023 

0014 

CQB6 

CAI 

7 6* 

0024 

0015 

F222 

JMP 

Fsrx 

0025 

0016 

C0B7 

CAI 

*7* 

0026 

0017 

F222 

JMP 

FSEVE 

0027 

0018 

COBS 

CAI 

*8* 

0028 

0019 

F222 

JMP 

FEFGH 

0029 

001 A 

C0B9 

CAI 

*9* 

0030 

00 IB 

F222 

JMP 

FMTNE 

0031 

001C 

C.0C1 

CAI 

•A* 

0032 

00 1 D 

F222 

JMP 

FTEM 

0033 

00 IE 

C0C2 

CAI 

•B* 

0034 

00 IF 

7222 

JMP 

FECE 

0035 

0020 

C0C3 

CAI 

‘C' 

0036 

0021 

F222 

JMP 

FTOVE 

0037 

0022 

C0C4 

CAI 

* D* 

0038 

0023 

F222 

JMP 

FTHD 

0039 

0024 

C0C5 

CAI 

E* 

0040 

0025 

F222 

JMP 

FFRTM 

0041 

0026 

C0C6 

CAI 

’F* 

0042 

0027 

F222 

JMP 

FFOTN 

0043 

0028 

C0C7 

CAI 

•G’ 

0044 

0029 

F222 

v JMP 

FSXTM 

0045 

002A 

FBAF 

■- J ST : 

%ERR. 

0046 

002B 

0110 

FZ R0 . Z AR 


0047 

002C 

9A9C 

STA 

CALL 

00.48 

00 2D 

F ei2A 

JMP 

- >00 P I 

0049 

002E 

0601 

Fg»E LAP 

r 

0050 

002F 

721E 

JMP 

. FST0 

0051 

0030 

C602 

FTW0 LAP 

.2 

0052 

0031 

F21C 

JMP 

FST0 

0053 

0032 

C603 

FTHRE LAP 

3 


0054 

0033 

F21A 


JMP 

FST0 

0055 

0034 

C604 

FF0UR 

LAP 

4 

0056 

0035 

F218 


JMP 

FST0 

0057 ' 

0036 

C605 

FFIVE 

LAP 

5 

0053 

0037 

F216 


JMP 

FST0 

0059 

0038 

C606 

FSIX 

LAP 

6 

0060 

0039 

F214 


JMP 

FST0 

0061 

003A 

C607 

FSEVE 

LAP 

7 

0062 

003B 

F212 


JMP 

FST0 

0063 

003C 

C608 

FEI GH 

LAP 

8 

0064 

003D 

F210 


JMP 

FST0 

0065 

003E 

C609 

FNINE 

LAP 

9 

0066 

003F 

F20E 


JMP 

FST0 

0067 

0040 

C60A 

FTElv 

LAP 

10 

0063 

0041 

F20C 


JMP 

FST0 

0069 

0042 

;60B 

FEL& 

LAP 

11 

0070 

0043 

F20A 


JMP 

FST0 

0071 

0044 

C60C 

FT WE 

LAP 

12 

0072 

0045 

F208 


JMP 

FST0 

. 0073 

0046 

C60D 

FTHD 

LAP 

13 

0074 

0047 

F206 


JMP 

FST0 

0075 

0048 

C60E 

FFRTW 

LAP 

14 

0076 

0049 

F204 


JMP 

FST0 

0077 

004A 

C60F 

FFVTW 

LAP 

15 

0073 

004B 

F202 


JMP 

FST0 

0079 

004C 

C610 

FSXTN 

LAP 

16 

003 0 

004D 

F200 


JMP 

FST0 

003 1 

004E 

9A7A 

FST0 

STA 

CALL 

003 2 

004F 

0310 


NAR 


008 3 

0050 

9A79 


STA 

FR0W 

003 4 

0051 

B26D 


LDA 

IB 

003 5 

0052 

DA77 

FR0N0 

. IMS 

FR017 

003 6 

0053 

F201 


JMP 

S.+ 2 

003 7 

0054 

F202 


JMP 

5+3 

0088 

0055 

8 A75 


ADD 

H255 

0089 

0056 

F604 


JMP 

FR0N0 

009 0 

0057 

9A67 


STA 

IB 

009 1 

0058 

B270 

F00P1 

LDA 

CALL 

009 2 

0059 

2101 


JAZ 

S+2 

009 3 

005A 

F204 


JMP 

5+5 

009 4 

005B 

B26C 


1, DA 

ZCT 

009 5 

dose 

0150 


IAR 


009 6 

005D 

9A63 


STA 

CMT 

009 7 

005E 

F203 


JMP 

2+4 

0098 

005F 

C7FF 


LAM 

255 

0099 

0060 

0150 


IAR 


0100 

0061 

9A5F 


STA 

CMT 

0101 

0062 

B25C 


LDA 

IB 

01*02 

0063 

9A5E 


STA 

IPB1 

0103 

0064 

B35D 


LDA 

*IPB1 

0x04 

0065 

9A5E 


STA 

MAX 

0105 

0066 

DA5B 

L00P2 

IMS 

I PB1 

0106 

0067 

B35A 


LDA 

’*IPS1 



0107 

0C68 

D25B 

CMS 

MAX • 

0108 

0069 

F20 1 

JMP 

S+2 ' 

0109 

006A 

9A59 

STA 

MAX 

0110 

006B 

DA55 

IMS 

CNT 

0111 

006C 

F606 

JMP 

L00P2 

01 12 

006D 

0350 

ABP 


0113 

006E 

9A56 

STA 

SCALE 

01 14 

006F 

9A56 

STA 

R0UND 

0115 

0070 

B256 

LDA 

H50 

01 16 

0071 

D252 

CMS 

MAX 

0117 

0072 

F202 

JMP 

S+ 3 

0118 

0073 

F20C 

JMP 

L00P 

0119 

0074 

F20B 

JMP 

L00P 

0120 

0075 

B24E 

LDA 

MAX 

0121. 

0076 

0108 

ZXR 


0122 

0077 

9 24F 

SUB 

H50 

0123 

0078 

0128 

IX R 


0124 

0079 

31C2 

JAG 

$-2 

0125 

007A 

EA4A 

SIX 

SCALE 

0126 

007B 

1200 

R0V 


0127 

007C 

1 IAS 

RRX 

1 

0128 

007D 

3201 

J0R 

$-5*2 

0129 

007E 

0128 

IXR 


0130 

007F 

EA46 

STX 

R0UMD 

0131 

008 0 

F355 L00P 

J ST 

*CBLF 

0132 

008 1 

FB54 

J ST 

*CRLF 

0133 

008 2 

C6ED 

LAP 

T = f 

0134 

OOS 3 

FB54 

J ST 

£0tL 

0135 

008 4' 

OOCD 

DATA 

SF 

0136 

008 5 

C6A0 

LAP 

f I 

0137 

008 6 

FB50 

J ST 

£0TT 

0138 

003 7 

B23D 

LDA 

SCALE 

0139 

0038 

FB4C 

J ST 

^ 0 DEC 

0140 

0089 

C60 1 

LAP 

1 

0141 

OOS A 

9 A4l 

STA 

R0N0 

0142 

003 B 

FB4A 

J ST 

*CHLF 

0143 

008 C 

FB49. 

J ST 

*CRLF 

0144 

008 D 

B23E 
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APPENDIX D 


DEF FiJRCX) = IMTC 100*X+0» 5)/ 100 
DIM 0C 255) > F( 255), XC 255) 

PRINT "SPECTRUM CORRECTED F0R NEGETIVE DIP" 

PRINT 
PRINT 
CALL C 20) 

CALL C 5, 0( 0), 256.> 2) 

LET FC 0) = 2C 0) 

FOR N= 1 TB 255 

L ET I 1=N- 1 

IF I 1> G. THEN 3 5 

LET I 1=255 

L ET I 2=N+ 1 

IF 1 2< 256 THEN IOC 

LET 12=12-255 

LET I3=W+24 

IF I 3< 256 THEM 115 

LET 13=13-255 

LET I 4=M+ 25 

IF I 4< 256 THEM 130 

LET 14=14-255 

LET C=C 0CI 3) + 0CI4)-0CI 1)-0CI2))/2O 
LET FCIM) = 0CN) + C 
NEXT N > 

PBIMT "A: 1 FOR DISPLAY, 2 FOR GRAPH, 3 F0R TAPE" 

PRINT "B: 0 FBR ORIGINAL, 1 F0R FINAL" 

PRINT 

PRINT "A=",* 

INPUT A 
PRINT "B="J 
INPUT B 

PRINT "RESOLUTION I>"J 
INPUT D 

PRI NT "IN I TI AL WAVE L EM GTH L 0= "l 
INPUT LO 

PRINT . 

PRINT ’ 

IF A=3 THEM 730 
LET 2=0 
L ET M= 1E10 
FOR N= 1 TO 255 


201 


202 


350 IF 3=1 THEJ 400 

360 LET XCN)=0CNO 

370 G0T0 4 1 0 

400 LET X<IJ> = FCN ) 

410 IF XCN)>=M THEN 430 

420 LET M=X(N) 

430 I F X (N ) <=Z THEN 450 

440 LET Z=XCN> 

450 NEXT N 

460 IF M > = 0 THEN 510 

470 LET SO 255/CZ-M) 

43 0 LET S 1= 43/CZ-N) 

490 LET YO-M*S 

500 G0T0 540 

510 LET 50= 255/Z 

520 LET Sl=43/Z 

530 LET Y 0= 0 

54 0 PHI N T "M AX= ",* Z ; "H IM= "J M 

550 IF A= 2 THEN 630 

560 PRINT "SCALE FACT0R=";SO 

570 F0R N= 1 T0 255 

530 LET E=INT< S0*XCN) + 0. 5) 

59 0 CALL C3,N,Y0, 2, E> 

600 NEXT N 

620 G0T0 160 

630 PRINT "SCALE FACT0R="; SI 

635 PRINT 

636 PRINT 

640 F0R N= 1 T0 255 

650 LET E=INT( Sl*XCM)+0. 5) 

660 LET L=L 0+ C N- I ) * D 

670 PRINT FNRCLJJ TABCS),*XCN> 

715 NEXT N 

730 IF B= 0 THEn) 760 

740 IF E= 1 THEN 73 0 

750 G0T0 200 

760 CPLLt 6, 0C 0), 256, 2) 

770 G0T0 79 0 

76 0 CALLC 6, FC 0), 256, 2) 

79 0 CALL C 6, 0, 0, 3) 

800 ST0P 
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